Question: Analyzing sub sampling data CRD ANOVA in R
1
2.0 years ago by
Wox320
HUJI
Wox320 wrote:

My experiment is CRD: Two Factor Factorial Experiment, Fixed Effects. Factor 1: Levels: a,b (2) Factor 2: Levels: p,q,r (3) for each experimental units, there are four subsamples (Observational units) I want to add subsampling to this model and do the analysis, As I've referred this is only possible in GenStat S/W (http://www.stats.net.au/dropbox/ANOVA%20&%20REML%20manual.pdf Page 48). How can I modify the R codes to get the above ANOVA output and mean separation?

``````stem_lengths <- c(3.5, 2.5, 3, 5, 3.5, 4.5, 5, 5.5, 5.5, 8.5, 6.5, 7, 6, 6, 6.5, 7, 6, 11, 4, 4.5, 3, 5.5, 3.5, 4, 4.5, 6, 4.5, 6, 7, 7, 5.5, 8.5, 6.5, 9, 7, 7, 3, 5.5, 2.5, 4, 3, 4, 5, 5, 6.5, 9, 8, 7, 3.5, 4.5, 8.5, 8.5, 7, 9, 4.5, 5, 3, 3.5, 4, 5, 4.5, 5, 5.5, 8.5, 6.5, 7, 7, 7.5, 7.5, 8.5, 7, 8)

stems <- data.frame( temp = rep(rep(as.factor(c("high", "low")), each = 9), 4), daylight = rep(rep(as.factor(c(8, 12, 16)), each = 3), 8), pot = rep(as.factor(1:18), 4), len = stem_lengths )

aov.stem <- aov(len ~ temp * daylight, data = stems)  aov.stemsummary(aov.stem)coefficients(aov.stem)
``````
R • 1.2k views
written 2.0 years ago by Wox320

Hello BIOAWY!

We believe that this post does not fit the main topic of this site.

For this reason we have closed your question. This allows us to keep the site focused on the topics that the community can help with.

If you disagree please tell us why in a reply below, we'll be happy to talk about it.

Cheers!

Thanks Ram. no worries. I just thought to post this Q in here since it's related to Biological data analysis / Design experiments thinking that most of the experts in here are familiar with this type of Q.

Thanks. I'll post this CV.

1

No problem - if someone thinks this belongs on the forum, we can always reopen it. Cheers!