Displaying 1 result from an estimated 1 matches for "species10".
Did you mean:
species1
2011 Jul 29
1
Using perm.t.test() upon Matrix/Dataframe columns parted by factor instead of t.test()
...uot;, 5))->Faktor
matrix(rnorm(100, mean=20, sd=4), nrow=10, ncol=10)->M
colnames(M) <- c("species1","species2",
"species3","species4","species5","species6","species7","species8","species9",
"species10")
###Conventional T-Test to test for differences of each species per factor
lapply(
M, function(x)
t.test(x~ Faktor)
)
###Trying it for perm.t.test without the helpful formular expression.
Does not work
lapply(
M, function(x)
perm.t.test(subset(M, Faktor=="A"),...