xyplot(log(Length)*Conc | Strip,muscle,as.table=T)
X<-model.matrix(~Strip-1,muscle)
class(X)
fm<-nls(log(Length)~cbind(X,X*rho^Conc),muscle,start=c(rho=0.4),trace=T,algorithm="plinear")
fm<-nls(log(Length)~cbind(X,X*rho^Conc),muscle,start=c(rho=0.8),trace=T,algorithm="plinear")
coef(fm)
fm<-nls(log(Length)~cbind(X,X*rho^Conc),muscle,start=c(rho=0.8),trace=T,algorithm="plinear")
fm<-nls(log(Length)~cbind(X,X*rho^Conc),muscle,start=c(rho=0.5),trace=T,algorithm="plinear")
coef(fm)
fm<-nls(log(Length)~cbind(X,X*rho^Conc),muscle,start=c(rho=0.05),trace=T,algorithm="plinear")
coef(fm)
names(fm)
mdata<-with(muscle,range(Conc))
mdata
mdata<-with(muscle,expand.grid(Strip=levels(Strip),Conc=seq(0.25,4,legth=50))
mdata
names(mdata)
q()
hist(CD)
hist(Price)
with(CD,hist(Price))
m30.mu
CD1$m30.mu
fm<-loess(Price~Load,CD)
help(loess)
find("petrol")
fixef(petrol.lm3)
fixef(petrol.lme)
ranef(petrol.lme)
with(petrol,SG)
with(petrol,EP)
help(lme)
