[R] ? ~ ADONIS
Tom Wright
tom at maladmin.com
Thu Mar 27 21:56:07 CET 2014
I'm no expert on the adonis functions by the error you are getting
usually means the dataset is not quite what you expected. Are you sure
the event dataset is loading as you expect?
On Tue, 2014-03-25 at 16:36 -0400, Maggie Wisniewska wrote:
> Hello,
> I am an R noivce, so excuse my simple questions .
> My data have 5 respone variales (each 109 rows) and three independent
> variables, each with two levels [Age Class(A,SA), Repro State(L,NL) and Sex
> of Calf(F,M. *some females do not have calves so NA)]. My data are not
> normally distributed with unequal number of samples (within each
> independent variale).
> example:
> AClass RState CSex Chemo/min Aggress/min Nonagress/min
> A L F 0.5 0.5
> 0.3
> SA L M 0.4 0.2
> 1
> A NL F 0 1
> 0.1
> A NL NA 2 0
> 0.3
> SA NL F 0.4 0.4
> 0
>
> I am trying to run ADONIS (as multivariate nonparametric anova) to test the
> differences and interactions.
> my code is the following:
> event<-read.csv("C:\\Users\\Owner\\Google
> Drive\\EventBehaviorResults\\Event.csv")
> adonisresults=adonis(Chemo.min~AClass+CSex+AClass*CSex,data=event,permutations=10,method="horn",strata=NULL,binary=FALSE)
> (and, separatelly for AClass and RState becaue I don't know how ot deal
> with NA values in )
> adonisresults=adonis(Chemo.min~AClass+RState+AClass*RState,data=event,permutations=10,method="horn",strata=NULL,binary=FALSE)
>
> I get the following error message and am not able to proceed
> Error in rowSums(x, na.rm = TRUE) :
> 'x' must be an array of at least two dimensions
>
> I would appreciate any input as to how to continue (and am happy to
> elaborate if more info would be of help).
> Thank you for readign.
> Maggie
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help at r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
More information about the R-help
mailing list