Displaying 1 result from an estimated 1 matches for "anlsx".
2012 Mar 26
1
Seeming failure of options(width=60)
...ests this should work.
options(width=60)
pastured <- data.frame(
time=c(9, 14, 21, 28, 42, 57, 63, 70, 79),
yield= c(8.93, 10.8, 18.59, 22.33, 39.35,
56.11, 61.73, 64.62, 67.08))
regmod<-"yield ~ t1 - t2*exp(-exp(t3+t4*log(time)))"
huetstart<-c(t1=70, t2=60, t3=0, t4=1)
anlsx<-try(nls(regmod, start=huetstart, trace=TRUE, data=pastured))
print(anlsx)
John Nash