Displaying 3 results from an estimated 3 matches for "fm1soy".
Did you mean:
fm1co2
2009 May 04
1
how to change nlme() contrast parametrization?
How to set the nlme() function to return the answer without the intercept parametrization?
#=========================================================================================
library(nlme)
Soybean[1:3, ]
(fm1Soy.lis <- nlsList(weight ~ SSlogis(Time, Asym, xmid, scal),
data = Soybean))
(fm1Soy.nlme <- nlme(fm1Soy.lis))
fm2Soy.nlme <- update(fm1Soy.nlme, weights = varPower())
plot(ranef(fm2Soy.nlme, augFrame = TRUE),
form = ~ Year * Variety, layout = c(3,1))
soyFix <...
2001 Jan 17
1
Pinheiro/Bates Soybean nlme failure
Dear Mixed Effect Friends,
Somehow, R(1021, Windows) seem to run differently from S Plus:
The soybean example from Pinheiro/Bates on page 290 fails
in R. (Soybean1 is Soybean with the NA and "critical" case
removed. Same procedure with full Soybean).
> fm1Soy.lis<-nlsList(weight~SSlogis(Time,Asym,xmid,scal),data=Soybean1)
> fm1Soy.nlme<-nlme(fm1Soy.lis)
Error: Singularity in backsolve at level 0, block 1
In addition: There were 50 or more warnings (use warnings() to see the first 50)
> warnings()
Warning messages:
1: Singular precision matri...
2004 Feb 20
1
nlme and multiple comparisons
This is only partly a question about R, as I am not quite sure about the
underlying statistical theory either.
I have fitted a non-linear mixed-effects model with nlme. In the fixed
part of the model I have a factor with three levels as explanatory
variable. I would like to use Tukey HSD or a similar test to test for
differences between these three levels.
I have two grouping factors: