Displaying 3 results from an estimated 3 matches for "deathday".
2006 Apr 16
0
[S] Problems with lme and 2 levels of nesting:Summary
...in each cup (the cups were not meant to differ in any way, it's just impossible to count 50 larvae in motion, but ok with 10). The day of death was recorded for each larva. Since there's only one data point per larva, I see no need of including larva ID in the model. So:
>
> Response: DeathDay
> Fixed factor: Treatment
> Random factors: Clucth, Cup
> Design: Cup nested in Clutch
>
> I'm primarlity interested in the effect of Treatment on DeathDay, but would also like to know about the variation within "Clutch" and "Cup". The unit of replication is &...
2006 Nov 11
2
Bayesian question (problem using adapt)
...uot;function undefined")}
val=c()
for (i in 1:length(x)) {
val=c(val, integrate(hazard2, lower=0, upper=x, theta=theta,
param=param)$value)}
return(val)
}
param0 = c(1.11, .833)
theta0 = c(1.11, .833)
yn1 <-read.table("yang1.txt", col.names = c("trtmt", "deathday", "reason"))
yn2 <-read.table("yang2.txt", col.names = c("trtmt", "deathday", "reason"))
for (i in 1:length(yn1$reason)) {
if (yn1$reason[i] == -1) (yn1$reason[i] = 0) else (yn1$reason[i] = 1) }
for (i in 1:length(yn2$reason)) {
if...
2006 Apr 08
1
dim(x) error message in lme (nlme package)
Hello
I am trying to analyse mortality data from fish larvae using lme from
the package nlme as well as using lmer in the package lme4
Response is DeathDay
Fixed factor is Treatment
Random factors are Clucth, Cup
Design: Cup nested in Clutch
If I do this in lme, I use the syntax:
model1 <- lme(DeathDay ~ Treatment, random=~ 1 | Clutch/Cup)
summary(model1)
I get the first part of the output, but where should be "Correlation of
Fixed Effects&...