[R] cluster, mona error
Sean Davis
sdavis2 at mail.nih.gov
Thu Jan 27 15:16:11 CET 2005
On Jan 27, 2005, at 9:06 AM, Morten Mattingsdal wrote:
> Hi
>
> I have a problem using the package cluster on my binary data. I want
> to try mona at first. But i get the an error.
>
> hc<-read.table("all.txt", header=TRUE, sep="\t", row.names=1)
> srt(hc)
> `data.frame': 51 obs. of 59 variables:
> $ G1p : int 2 1 1 1 1 1 1 1 1 1 ...
> $ G1q : int 1 1 1 1 1 1 1 1 1 1 ...
> $ G2p : int 1 1 1 1 1 1 1 1 1 1 ...
> $ G2q : int 1 1 1 1 1 1 1 1 1 1 ...
> $ G3p : int 1 1 1 1 1 1 1 1 1 1 ...
>
> m<-mona(hc)
> Error in mona(hc) : All variables must be binary (factor with 2
> levels).
>
You have to be careful that the data are indeed each factors with 2
levels (numeric variables with values 1 and 2 will not do). A summary
of the data will tell you that.
Sean
> I find this strange when the cluster dataset "animals" have the same
> structure as my data.
>
> srt(animals)
> `data.frame': 20 obs. of 6 variables:
> $ war: int 1 1 2 1 2 2 2 2 2 1 ...
> $ fly: int 1 2 1 1 1 1 2 2 1 2 ...
> $ ver: int 1 1 2 1 2 2 2 2 2 1 ...
> $ end: int 1 1 1 1 2 1 1 2 2 1 ...
> $ gro: int 2 2 1 1 2 2 2 1 2 1 ...
> $ hai: int 1 2 2 2 2 2 1 1 1 1 ...
>
> m<-mona(animals) #works fine
>
> what is this error trying to tell me?
> mvh
> morten
>
> ______________________________________________
> R-help at stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide!
> http://www.R-project.org/posting-guide.html
More information about the R-help
mailing list