Displaying 1 result from an estimated 1 matches for "log_abund".
Did you mean:
log2abun
2009 Nov 05
1
help with ols and contrast functions in Design library
...rast function in the same library to test various contrasts.
As a simple example, suppose I have three factors: feature (3
levels), group (2 levels), and patient (3 levels). Patient is coded
as a non-unique identifier and is therefore nested within group.
response <- rnorm(length(example$LOG_ABUNDANCE), mean = 12)
feature <- rep(c(1,2,3), 6)
group <- c(rep(c(1,2),each=9))
patient <- rep(rep(c(1,2,3), each=3),2)
myData <- data.frame(patient=factor(patient), group=factor(group),
feature=factor(feature), response=response)
I use the ols command to fit the linear model, but I recei...