Displaying 5 results from an estimated 5 matches for "trtmt".
2009 Sep 14
0
invalid type (list) for variable
...I cannot, however, figure
out how to add interactions between continuous*factor terms. I applied
a smoother to my continuous term because it has a non-linear trend
over time. Below is an example of my attempt, followed by the error
message:
MMF.Test <- gamm(MaxDepth ~
s(TripTime) + Trtmt +
s(TripTime):Trtmt,
random = list(Seal =~1), method = "REML",
data = EsealsAll)
"Error in model.frame.default(formula = MaxDepth ~ Trtmt +
s(TripTime):Trtmt + :
invalid type (list) for variable 's(TripTime)'
In addition: Warning messa...
2009 Nov 09
1
Getting Sphericity Tests for Within Subject Repeated Measure Anova (using "car" package)
...but
have no luck with Repeated-Measures Assuming Sphericity or Separate
Sphericity Tests.
I am trying to follow example from "car" package, but it seems that I am not
getting something right.
> Dataset$Sessn <- as.factor(Dataset$Sessn)
> LinearModel.1 <- lm(Response ~ Sessn*Trtmt, data=Dataset)
> summary(LinearModel.1)
All, good so far, but I have problem understanding "idata=" and "idesign="
functions pertaining to my example. Session is my repeated measure (Sessn 1
and Sessn 2 = two sessions, in reality I have more) and it is already
stacked. Any...
2009 Nov 09
1
Getting Sphericity Tests for Within Subject Repeated Measure Anova (using "car" package) (Adjusted Dataset)
...but
have no luck with Repeated-Measures Assuming Sphericity or Separate
Sphericity Tests.
I am trying to follow example from "car" package, but it seems that I am not
getting something right.
> Dataset$Sessn <- as.factor(Dataset$Sessn)
> LinearModel.1 <- lm(Response ~ Sessn*Trtmt, data=Dataset)
> summary(LinearModel.1)
All, good so far, but I have problem understanding "idata=" and "idesign="
functions pertaining to my example. Session is my repeated measure (Sessn 1
and Sessn 2 = two sessions, in reality I have more) and it is already
stacked. Any...
2006 Nov 11
2
Bayesian question (problem using adapt)
...<=0) {return("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:leng...
2009 May 01
3
adding zeros to dataframe
Greetings,
I am new to R and am hoping to get some tips from experienced R-programmers.
I have a dataset that I've read into R as a dataframe. There are 5 columns: Plot location,species name, a species number code (unique to each species name), abundance, and treatment. There are 272 plots in each treatment, but only the plots in which the species was recorded have an abundance value. For