Entering edit mode
3.5 years ago
Claire
•
0
I have 30 row samples. The first column is the names of my samples. My code then ran like this:
ENV = read.table("E:--\lastmd.txt",header=TRUE,sep="\t",row.names=1,dec=",")
aov.res <- aov(activity ~ genotype*gender)
> Df Sum Sq Mean Sq F value Pr(>F)
genotype 2 87842 43921 241.125 1.011e-13 ***
gender 2 689096 344548 1891.565 < 2.2e-16 ***
genotype :gender 4 15402 3850 21.139 1.331e-06 ***
Residuals 18 3279 182
LSD_A = LSD.test(y = activity,
trt = genotype,
DFerror = aov.res$df.residual,
MSerror = deviance(aov.res)/aov.res$df.residual,
alpha = 0.05,
p.adj = "bonferroni",
group = TRUE,
console = TRUE)
Here i get this error:
> Error in (1 - h) * qs[i] : non numeric argument to binary operator
> mean.default(junto[, 1]) :argument not numeric