source("weight_functions.s")
source("commands.s")
?nls
?nls.control
 ?nls
?nls.control
source("commands.s")
?nls.control
source("commands.s")
ls()
summary(ty.mwc)
summary(ty.mwc.wt)
ls()
summary(ty.ad.wt)
?lines
?lowess
source("commands.s")
source("commands.s")
source("commands.s")
source("commands.s")
source("commands.s")
source("commands.s")
source("commands.s")
source("commands.s")
source("commands.s")
source("commands.s")
source("commands.s")
source("commands.s")
help.search("smooth")
(.packages)
?library
(.packages())
?smooth
?smooth.spline
source("commands.s")
source("commands.s")
?smooth.spline
     y18 <- c(1:3,5,4,7:3,2*(2:5),rep(10,4))
     xx  <- seq(1,length(y18), len=201)
     (s2  <- smooth.spline(y18)) # GCV
     (s02 <- smooth.spline(y18, spar = 0.2))
     plot(y18, main=deparse(s2$call), col.main=2)
     lines(s2, col = "gray"); lines(predict(s2, xx), col = 2)
     lines(predict(s02, xx), col = 3); mtext(deparse(s02$call), col = 3)
source("commands.s")
source("commands.s")
source("commands.s")
source("commands.s")
source("commands.s")
source("commands.s")
source("commands.s")
source("commands.s")
source("commands.s")
source("commands.s")
source("commands.s")
source("commands.s")
source("commands.s")
?points
?points
source("commands.s")
source("commands.s")
source("commands.s")
source("commands.s")
ls()
fitted(ty.mwc)
fitted(ty.mwc.wt)
y.mwc.wt
y.mwc
ls()
fitted(ty.ad.wt)
ls(,2)
y*(1-y)
residuals(ty.mwc.wt)
residuals(ty.mwc.wt)*y*(1-y)
y
ls()
y.mwc
fitted(ty.mwc)
y-residuals(yy.mwc)
y-residuals(ty.mwc)
ls()
wt.adair
y*(1-y)
ls()
wt.mwc
source("commands.s")
sum(residuals(ty.mwc.wt)*y*(1-y))
sum(residuals(ty.ad.wt)*y*(1-y))
sum(residuals(ty.mwc))
sum(residuals(ty.mwc)^2)
sum(residuals(ty.ad.wt)*y*(1-y)^2)
sum(residuals(ty.mwc.wt)*y*(1-y)^2)
sum((residuals(ty.mwc.wt)*y*(1-y))^2)
sum((residuals(ty.ad.wt)*y*(1-y))^2)
sum(residuals(ty.mwc)^2)
source("commands.s")
?nls
source("commands.s")
source("commands.s")
summary(ty.ad.wt)
wt.adair
var(log(y/(1-y)))
var(y)k
var(log(y/(1-y)))/varf(y)
var(log(y/(1-y)))/var(y)
ls()
sum((y-y.mwc)^2)
sum((y-y.mwc.wt)^2)
sum((y-y.mwc.wt)^2)*59.9
sum((y-y.ad.wt)^2)*59.9
sum((y-y.ad.wt)^2)
sum((y-fitted(ty.mwc.wt))^2)*59.9
sum((y-fitted(ty.mwc.wt))^2)
sum((y-fitted(ty.mwc))^2)
nlsModel
nls
history()
c(y,fitted(ty.mwc),fitted(ty.mwc.wt),fitted(ty.ad.wt))
cbind(y,fitted(ty.mwc),fitted(ty.mwc.wt),fitted(ty.ad.wt))
cbind(y,residuals(ty.mwc),residuals(ty.mwc.wt),residuals(ty.ad.wt))
history()
cbind(y,residuals(ty.mwc.wt),fitted(ty.mwc.wt),residuals(ty.ad.wt),fitted(ty.ad.wt))
sum(residuals(ty.mwc.wt))
sum(residuals(ty.mwc))
sum(residuals(ty.ad.wt))
sum(residuals(ty.mwc)^2)
sum(residuals(ty.ad.mwc)^2)
sum(residuals(ty.mwc.wt)^2)
sum((residuals(ty.mwc.wt)*y*(1-y))^2)
sum((y - residuals(ty.mwc.wt)*y*(1-y))^2)
cbind(y,y.mwc,y.mwc.wt,y.ad.wt)
cbind(y,fitted.values(ty.mwc),(y - residuals(ty.mwc.wt)*y*(1-y)),(y - residuals(ty.ad.wt)*y*(1-y)))
l
sum((y - residuals(ty.ad.wt)*y*(1-y))^2)
((y - residuals(ty.ad.wt)*y*(1-y))^2)
cbind((y - residuals(ty.ad.wt)*y*(1-y))^2)
cbind(y,(y - residuals(ty.ad.wt)*y*(1-y))^2)
cbind(y,(y - residuals(ty.ad.wt)*y*(1-y)))
sum((2y - residuals(ty.ad.wt)*y*(1-y))^2)
sum((2*y - residuals(ty.ad.wt)*y*(1-y))^2)
sum((residuals(ty.ad.wt)*y*(1-y))^2)
sum((residuals(ty.mwc.wt)*y*(1-y))^2)
sum((residuals(ty.mwc))^2)
cbind(y,(residuals(ty.ad.wt)*y*(1-y)))
cbind(y,residuals(ty.mwc),(residuals(ty.ad.wt)*y*(1-y)))
cbind(y,residuals(ty.mwc),y.ad.wt,(residuals(ty.ad.wt)*y*(1-y)))
cbind(y,y.mwc,residuals(ty.mwc),y.ad.wt,(residuals(ty.ad.wt)*y*(1-y)))
history()
c(sum((residuals(ty.mwc))^2), sum((residuals(ty.mwc.wt)*y*(1-y))^2), sum((residuals(ty.ad.wt)*y*(1-y))^2))
c(summary(ty.mwc)$sigma)
(summary(ty.mwc)
)
(summary(ty.mwc)$sigma)^2
(summary(ty.mwc)$sigma)^2*17
(summary(ty.mwc.wt)$sigma)^2*17
(summary(ty.mwc.wt)
)
(summary(ty.ad.wt))
(summary(ty.ad.wt)$sigma)^2*16
(summary(ty.ad.wt)$sigma)^2
(summary(ty.mwc.wt)$sigma)^2
(summary(ty.mwc.wt)$sigma)^2*17
(summary(ty.mwc)$sigma)^2*17
ty.mwc$sigma
ty.mwc$m$sigma
ty.mwc
list(ty.mwc)
names(ty.mwc)
list(ty.mwc$m)
list(ty.mwc$data)
list(ty.mwc$call)
names(ty.mwc)
list(ty.mwc$dataClasses)
list(ty.mwc$m)
list(ty.mwc$m$resid)
list(ty.mwc$m$resid())
list(ty.mwc$m$Rmat())
list(ty.mwc$m$conv())
list(ty.mwc$m$sigma())
list(ty.mwc$m$deviance())
list(ty.mwc.wt$m$deviance())
(ty.mwc.wt$m$deviance())
(summary(ty.mwc)$sigma)
(summary(ty.mwc.wt)$sigma)
(summary(ty.mwc.wt)$sigma)^2*17
(summary(ty.mwc.wt)
)
vcov((ty.mwc.wt)
)
6.554661e+10*0.0653
6.554661e+10*0.0653^2*17
0.0653^2*17
sqrt(6.554661e+10*0.0653^2*17)
summary(ty.mwc)
vcov(ty.mwc)
vcov(ty.mwc)[1,1]
sqrt(vcov(ty.mwc)[1,1])
sqrt(vcov(ty.mwc.wt)[1,1])
sqrt(diag(vcov(ty.mwc.wt)))
summary(ty.ad.wt)
nls(rep(0,length(y)) ~ wt.adair(y, p, K1,K2,K3,K4),tyuma.df,list(K1=.024,K2=.0743,K3=.0858,K4=7.37),control=nls.control(maxiter=100,minFactor=.001),trace=T)
vi(wt.adair)
wt.adair.mod<-vi(wt.adair)
nls(rep(0,length(y)) ~ wt.adair.mod(y, p, K1,K2,K3,K4),tyuma.df,list(K1=.024,K2=.0743,K3=.0858,K4=7.37),control=nls.control(maxiter=100,minFactor=.001),trace=T)
nls(rep(0,length(y)) ~ wt.adair.mod(y, p, K1,K2,K3,K4),tyuma.df,list(K1=.24,K2=.0743,K3=.0858,K4=7.37),control=nls.control(maxiter=100,minFactor=.001),trace=T)
nls(rep(0,length(y)) ~ wt.adair(y, p, K1,K2,K3,K4),tyuma.df,list(K1=.24,K2=.0743,K3=.0858,K4=7.37),control=nls.control(maxiter=100,minFactor=.001),trace=T)
nls(rep(0,length(y)) ~ wt.adair.mod(y, p, K1,K2,K3,K4),tyuma.df,list(K1=.024,K2=.0743,K3=.0858,K4=7.37),control=nls.control(maxiter=100,minFactor=.001),trace=T)
nls(rep(0,length(y)) ~ wt.adair.mod(y, p, K1,K2,K3,K4),tyuma.df,list(K1=.02,K2=.07,K3=.08,K4=7.),control=nls.control(maxiter=100,minFactor=.001),trace=T)
nls(rep(0,length(y)) ~ wt.adair(y, p, K1,K2,K3,K4),tyuma.df,list(K1=.02,K2=.07,K3=.08,K4=7.),control=nls.control(maxiter=100,minFactor=.001),trace=T)
nls(rep(0,length(y)) ~ wt.adair(y, p, K1,K2,K3,K4),tyuma.df,list(K1=.04,K2=.07,K3=.08,K4=7.),control=nls.control(maxiter=100,minFactor=.001),trace=T)
nls(rep(0,length(y)) ~ wt.adair.mod(y, p, K1,K2,K3,K4),tyuma.df,list(K1=.04,K2=.07,K3=.08,K4=7.),control=nls.control(maxiter=100,minFactor=.001),trace=T)
nls(rep(0,length(y)) ~ wt.adair.mod(y, p, K1,K2,K3,K4),tyuma.df,list(K1=.04,K2=.07,K3=.08,K4=7.),control=nls.control(maxiter=100,minFactor=.0001),trace=T)
nls(rep(0,length(y)) ~ wt.adair.mod(y, p, K1,K2,K3,K4),tyuma.df,list(K1=.04,K2=.07,K3=.08,K4=7.),control=nls.control(maxiter=200,minFactor=.001),trace=T)
nls(rep(0,length(y)) ~ wt.adair.mod(y, p, K1,K2,K3,K4),tyuma.df,list(K1=.04,K2=.07,K3=.08,K4=7.),control=nls.control(maxiter=200,minFactor=.0001),trace=T)
nls(rep(0,length(y)) ~ wt.adair.mod(y, p, K1,K2,K3,K4),tyuma.df,list(K1=.01,K2=.07,K3=.08,K4=7.),control=nls.control(maxiter=200,minFactor=.0001),trace=T)
nls(rep(0,length(y)) ~ wt.adair.mod(y, p, K1,K2,K3,K4),tyuma.df,list(K1=.02,K2=.07,K3=.08,K4=7.),control=nls.control(maxiter=200,minFactor=.0001),trace=T)
nls(rep(0,length(y)) ~ wt.adair.mod(y, p, K1,K2,K3,K4),tyuma.df,list(K1=.02,K2=.7,K3=.08,K4=7.),control=nls.control(maxiter=200,minFactor=.0001),trace=T)
nls(rep(0,length(y)) ~ wt.adair.mod(y, p, K1,K2,K3,K4),tyuma.df,list(K1=.2,K2=.7,K3=.08,K4=7.),control=nls.control(maxiter=200,minFactor=.0001),trace=T)
nls(rep(0,length(y)) ~ wt.adair.mod(y, p, K1,K2,K3,K4),tyuma.df,list(K1=.1,K2=.1,K3=.2,K4=7.),control=nls.control(maxiter=200,minFactor=.0001),trace=T)
nls(rep(0,length(y)) ~ wt.adair.mod(y, p, K1,K2,K3,K4),tyuma.df,list(K1=.1,K2=.1,K3=.2,K4=1000.),control=nls.control(maxiter=200,minFactor=.0001),trace=T)
nls(rep(0,length(y)) ~ wt.adair.mod(y, p, K1,K2,K3,K4),tyuma.df,list(K1=.01,K2=.1,K3=.001,K4=100.),control=nls.control(maxiter=200,minFactor=.0001),trace=T)
nls(rep(0,length(y)) ~ wt.adair.mod(y, p, K1,K2,K3,K4),tyuma.df,list(K1=.01,K2=.1,K3=.001,K4=1000.),control=nls.control(maxiter=200,minFactor=.0001),trace=T)
nls(rep(0,length(y)) ~ wt.adair.mod(y, p, K1,K2,K3,K4),tyuma.df,list(K1=.01,K2=.1,K3=.1,K4=1000.),control=nls.control(maxiter=200,minFactor=.0001),trace=T)
nls(rep(0,length(y)) ~ wt.adair.mod(y, p, K1,K2,K3,K4),tyuma.df,list(K1=.01,K2=.1,K3=.1,K4=10.),control=nls.control(maxiter=200,minFactor=.0001),trace=T)
nls(rep(0,length(y)) ~ wt.adair.mod(y, p, K1,K2,K3,K4),tyuma.df,list(K1=.01,K2=.1,K3=.1,K4=40.),control=nls.control(maxiter=200,minFactor=.0001),trace=T)
nls(rep(0,length(y)) ~ wt.adair(y, p, K1,K2,K3,K4),tyuma.df,list(K1=.01,K2=.1,K3=.1,K4=40.),control=nls.control(maxiter=200,minFactor=.0001),trace=T)
nls(rep(0,length(y)) ~ wt.mwc(y, p, K1,K2,K3,K4),tyuma.df,list(K1=.01,K2=.1,K3=.1,K4=40.),control=nls.control(maxiter=200,minFactor=.0001),trace=T)
sqrt(0.04652726/16)
history()
q()
ls()
sumdelsq.mwc<-sum(y - y.mwc)^2)
sumdelsq.mwc.wt<-sum(y - y.mwc.wt)^2)
sumdelsq.ad.wt<-sum(y - y.ad.wt)^2)
c(sumdelsq.mwc,sumdelsq.mwc.wt,sumdelsq.ad.wt)
sumdelsq.mwc<-sum((y - y.mwc)^2)
sumdelsq.mwc.wt<-sum((y - y.mwc.wt)^2)
sumdelsq.ad.wt<-sum((y - y.ad.wt)^2)
c(sumdelsq.mwc,sumdelsq.mwc.wt,sumdelsq.ad.wt)
attach(tyuma.df)
sumdelsq.mwc<-sum((y - y.mwc)^2)
sumdelsq.mwc.wt<-sum((y - y.mwc.wt)^2)
sumdelsq.ad.wt<-sum((y - y.ad.wt)^2)
c(sumdelsq.mwc,sumdelsq.mwc.wt,sumdelsq.ad.wt)
c(sumdelsq.mwc,sumdelsq.mwc.wt,sumdelsq.ad.wt)/c(17,17,16)
sqrt(c(sumdelsq.mwc,sumdelsq.mwc.wt,sumdelsq.ad.wt)/c(17,17,16))
c(sumdelsq.mwc=sumdelsq.mwc,sumdelsq.mwc.wt,sumdelsq.ad.wt)
c(sumdelsq.mwc=sumdelsq.mwc,sumdelsq.mwc.wt=sumdelsq.mwc.wt,sumdelsq.ad.wt=sumdelsq.ad.wt)
sqrt(c(sumdelsq.mwc,sumdelsq.mwc.wt,sumdelsq.ad.wt)/c(17,17,16))
sumdelsq.vec<-c(sumdelsq.mwc=sumdelsq.mwc,sumdelsq.mwc.wt=sumdelsq.mwc.wt,sumdelsq.ad.wt=sumdelsq.ad.wt)
sqrt(sumdelsq.vec/c(17,17,16))
?help.search("cat")
help.search("cat")
help.search("output")
?capture.output
?print
?cat
source("commands.s")
source("commands.s")
source("commands.s")
sink("scr.txt")
postscript(file="scr.ps",horizontal=FALSE)
source("commands.s")
sink()
graphics.off()
q()
source("commands.s")
sink("scr.txt")
source("commands.s")
sink()
q()
sum((y - y.mwc)^2/(y*(1-y))/17
)
sum((y - y.mwc)^2/(y*(1-y)/17)
)
attach("tyuma.df")
ls()
attach(tyuma.df)
sum(((y - y.mwc)/(y*(1-y))^2)/17)
sum(((y - y.mwc.wt)/(y*(1-y))^2)/17)
sum(((y - y.ad.wt)/(y*(1-y))^2)/16)
sumdelsq.wt.vec<-c(sum(((y - y.mwc)/(y*(1-y))^2)), sum(((y - y.mwc.wt)/(y*(1-y))^2)), sum(((y - y.ad.wt)/(y*(1-y))^2)))/c(17,17,16)
sumdelsq.wt.vec
sqrt(sumdelsq.wt.vec)
sumdelsq.wt.vec<-c(sum(((y - y.mwc)/(y*(1-y)))^2), sum(((y - y.mwc.wt)/(y*(1-y)))^2), sum(((y - y.ad.wt)/(y*(1-y)))^2))/c(17,17,16)
sumdelsq.wt.vec
sqrt(sumdelsq.wt.vec)
q()
ls()
sumdelsq.wt.vec
sumdelsq.vec
sqrt(sumdelsq.wt.vec)
sqrt(sumdelsq.vec)
sumdelsq.wt.vec<-c(sum(((y - y.mwc)/(y*(1-y)))^2), sum(((y - y.mwc.wt)/(y*(1-y)))^2), sum(((y - y.ad.wt)/(y*(1-y)))^2))
attach(tyuma.df)
sumdelsq.wt.vec<-c(sum(((y - y.mwc)/(y*(1-y)))^2), sum(((y - y.mwc.wt)/(y*(1-y)))^2), sum(((y - y.ad.wt)/(y*(1-y)))^2))
sumdelsq.vec
sumdelsq.wt.vec
sumdelsq.wt.vec/c(17,17,16)
sumdelsq.vec/c(17,17,16)
sqrt(sumdelsq.wt.vec/c(17,17,16))
sqrt(sumdelsq.vec/c(17,17,16))
history()
ls()
sumdelsq.wt.vec<-c(mwc=sum(((y - y.mwc)/(y*(1-y)))^2), mwc.wt=sum(((y - y.mwc.wt)/(y*(1-y)))^2), ad.wt=sum(((y - y.ad.wt)/(y*(1-y)))^2))
sumdelsq.wt.vec
sqrt(sumdelsq.wt.vec/c(17,17,16))
source("commands.s")
source("commands.s")
sink("scr,txt")
source("commands.s")
sink()
ls()
qqnorm(ty.mwc.wt)
?qqnorm
qqnorm(y)
qqnorm(log(y/(1-y)))
plot(residuals(ty.mwc))
plot(residuals(ty.mwc.wt))
plot(residuals(ty.mwc))
lines(residuals(ty.mwc.wt))
ls()
plot(residuals(ty.mwc))
lines((y-y.mwc)
)
lines((y-y.mwc.wt))
plot((y-y.mwc.wt))
lines((y-y.mwc.wt))
lines((y-y.mwc))
lines((y-y.ad.wt))
qqnorm((y-y.ad.wt))
qqnorm((y-y.mwc.wt))
qqnorm((y-y.mwc))
qqline((y-y.mwc))
qqnorm((y-y.mwc.wt))
qqline(())
qqnorm(residuals(ty.mwc))
qqnorm(residuals(ty.mwc.wt))
qqnorm(residuals(ty.ad.wt))
qqline(residuals(ty.ad.wt))
qqnorm(residuals(ty.mwc)); qqline(residuals(ty.mwc))
readline()
readline()
qqnorm(residuals(ty.mwc)); qqline(residuals(ty.mwc))
qqnorm(residuals(ty.mwc.wt)); qqline(residuals(ty.mwc.wt))
qqnorm(residuals(ty.ad.wt)); qqline(residuals(ty.ad.wt))
> qqnorm((y-y.mwc))
qqnorm((y-y.mwc))
qqline((y-y.mwc))
qqnorm((y-y.mwc)); qqline((y-y.mwc))
qqnorm((y-y.mwc.wt)); qqline((y-y.mwc.wt))
qqnorm((y-y.ad.wt)); qqline((y-y.ad.wt))
qqnorm((y-y.mwc)); qqline((y-y.mwc))
qqnorm(residuals(ty.mwc)); qqline(residuals(ty.mwc))
qqnorm(residuals(ty.mwc.wt)); qqline(residuals(ty.mwc.wt))
q()
