Displaying 1 result from an estimated 1 matches for "hfdata".
Did you mean:
fdata
2008 Dec 12
0
Help with a permutation test
...permute the MNNUM column once
rand<- function(DF){
new.DF<-DF
new.DF$MNNUM<-sample(new.DF$MNNUM)
new.DF
}
#this function does one model I am interested in.
modeltree<-function(DF){
MLM.plot <- multinom(MN_fact ~ Canpy + mean_dbh + num_beechoak +
num_class5 + prop_hard , data=hfdata, trace=FALSE)
MLM.plot
}
# this replicates the 'rand' function and applies a model
resamp.funct<-function(DF,funct, n){
list<-replicate(n,rand(DF), simplify = FALSE)
sapply(list, funct, simplify = FALSE)
}
#So if I paste below:
l<-resamp.funct(hfdata, modeltree, 3)
#...