Bioanalyzer HS results in.
Each sample has a tiny peak at 68-69 bp. Are these the adapters? What happened to DNA???
Meeting with Nick
Robust regression
Concerned that simple lm
maybe too liberal in detecting patterns. After some reading, decided that I should use robust regression for analysis of gene expression data. No default F statistic to assess overall model fit, but this page explains that fitting a full model and a minimal model with no covariates Y ~ 1
and then doing anova(full.model, minimal.model)
is equivalent to an F test. I confirmed this for the simple lm
case.
Error message where model doesn’t fit - need to continue ddply
even with error. Using try
as suggested here.
Hof, C., Rahbek, C. & Araújo, M.B. (2010). Phylogenetic signals in the climatic niches of the world’s amphibians. Ecography, 33, 242–250.
Botkin, D., Saxe, H., Araujo, M.B., Betts, R., Bradshaw, R.H.W., Cedhagen, T., Chesson, P., Dawson, T.P., Etterson, J.R., Faith, D.P., Ferrier, S., Guisan, A., Hansen, A., Hilbert, D.W., Loehle, C., Margules, C., New, M., Sobel, M.J. & Stockwell, D.R.B. (2007). Forecasting the Effects of Global Warming on Biodiversity. BioScience, 57, 227–236.
This work is licensed under a Creative Commons Attribution 4.0 International License.