Analyzing sub sampling data CRD ANOVA in R

Question

[WUSCHEL]

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)

[Ram]

Hello BIOAWY!

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

Please ask Cross Validated

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!

[WUSCHEL]

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.

[Ram]

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