search for: fm5

Displaying 7 results from an estimated 7 matches for "fm5".

Did you mean: fc5
2012 Jan 17
1
MuMIn package, problem using model selection table from manually created list of models
...on object using dredge fm0 <- lm(y ~ ., data = Cement) print(dd <- dredge(fm0)) fm1 <- lm(formula = y ~ X1 + X2, data = Cement) fm2 <- lm(formula = y ~ X1 + X2 + X4, data = Cement) fm3 <- lm(formula = y ~ X1 + X2 + X3, data = Cement) fm4 <- lm(formula = y ~ X1 + X4, data = Cement) fm5 <- lm(formula = y ~ X1 + X3 + X4, data = Cement) # ranked with AICc by default # obviously this works model.avg(get.models(dd, delta < 4)) # option 2: the aim is to produce a model selection object comparable to that from get.models(dd, delta < 4) # but from a manually-specified list of...
2012 Dec 16
1
nls for sum of exponentials
...following model with a sum of exponentials - y ~ Ae^(-md) + B e^(-nd) + c the model has 5 parameters A, b, m, n, c I am using nls to fit the data and I am using DEoptim package to pick the most optimal start values - fm4 <- function(x) x[1] + x[2]*exp(x[3] * -dist) + x[4]*exp(x[5] * -dist) fm5 <- function(x) sum((wcorr-fm4(x))^2) fm6 <- DEoptim(fm5, lower=c(0,0.1,1,0.1,1), upper=c(10e7, 100, 300,100, 300 ), control=list(trace=FALSE)) par2 <- fm6$optim$bestmem names(par2) <- c("c", "A", "n1", "B", "n2") fit2 <- nls(wcorr...
2004 Apr 22
1
lme correlation structure error
...ing identical to that in the book. fm2<-update(fm1,correlation=corAR1()) #### The parameters of fm2 are different to that in the book and plot(ACF(fm2)) ##### signifies that serial correlation still exists in the residuals. ###### When I try and run this (which runs fine in the book) fm5<-update(fm1,corr=corARMA(p=1,q=1)) #### I get the following error message #Error in "coef<-.corARMA"(*tmp*, value = c(62.3428455941166, 62.3428517930051 : # Coefficient matrix not invertible I have tried running the example on R for windows versions 1.7.1 and 1.8.1 with...
2006 Sep 06
1
Help on estimated variance in lme4
Dear all, I get an error message when I run my model and I am not sure what to do about it. I try to determine what factors influence the survival of voles. I use a mixed-model because I have several voles per site (varying from 2 to 19 voles). Here is the model: ### fm5 <-lmer(data=cdrgsaou2, alive~factor(pacut)+factor(agecamp)+factor(sex)+ResCondCorp+(1|factor(cdrgsa ou2$ids)), family=binomial, method="Laplace", ) ### Description of variables Alive: 0 or 1; dead or alive pacut: 0 or 1; presence of parasites agecamp: a or j; adult or juvenile sex:...
2003 Oct 15
2
Example of cell means model
This is an example from chapter 11 of the 6th edition of Devore's engineering statistics text. It happens to be a balanced data set in two factors but the calculations will also work for unbalanced data. I create a factor called 'cell' from the text representation of the Variety level and the Density level using '/' as the separator character. The coefficients for the linear
2006 Oct 20
1
Translating lme code into lmer was: Mixed effect model in R
...th ~ year +(year|schoolid:childid) +(year|schoolid), egsingle, control=list(gradient = FALSE, niterEM = 0)) Both result in parameter estimates that are exactly the same. The newest release of lme4 allows for an even easier transition between lme and lmer. Model syntax in lmer can be specified as: fm5 <- lmer(math ~ year +(year|schoolid/childid), egsingle, control=list(gradient = FALSE, niterEM = 0)) Notice how the grouping structure is used in this case, which resembles the lme syntax for random terms. The reason that the control statement is used is because Doug Bates told me to. I hope h...
2009 Jan 20
0
"syslog is smaller" warning
Hi list, I've installed logcheck some time ago. I also added some ignore filters, but there's one thing I can't figure out. Every day I get the following mail from logcheck: System Events =-=-=-=-=-=-= *** WARNING ***: Log file /var/log/syslog is smaller than last time checked! *************** *************** This could indicate tampering. I guess it has something to do with the