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