Hi,
I was wondering if anyone could help me address this issue? I am trying to use the treeWAS package on R to perform SNP analysis. I have prepared my two data-sets (one is a matrix (geney) the other is factor (pheny)) as required by the package, and then ran the following treeWAS command. I keep getting the "Error in axis" comment below:
> out <- treeWAS(snps = geney,
+ phen = pheny,
+ na.rm=TRUE,
+ tree = "BIONJ",
+ seed = 1)
Error in axis(side = side, at = x, labels = lab, ...) :
no locations are finite
When I use the two example data-sets, provided with the package, the command works fine. I have checked the two data-sets I have created against the two examples in the package and the class, attributes and formatting match. So I'm not sure what the error means. The following is a short description of my data:
> class(geney)
[1] "matrix"
> str(geney)
num [1:27, 1:62] NA NA NA NA NA 0 NA NA NA NA ...
- attr(*, "dimnames")=List of 2
..$ : chr [1:27] "FR3883.vcf" "FR3977.vcf" "FR4887.vcf" "FR4928.vcf" ...
..$ : chr [1:62] "X1000756.G" "X1000756.t" "X1001126.g" "X1001126.T" ...
> class(pheny)
[1] "factor"
> str(pheny)
Factor w/ 2 levels "A","B": 1 1 1 1 1 1 1 1 1 1 ...
- attr(*, "names")= chr [1:27] "FR4887.vcf" "FR4928.vcf" "FR4930.vcf" "FR5128.vcf" ...
This is the description of the exemplary data-sets:
> class(snps)
[1] "matrix"
> str(snps)
num [1:100, 1:20003] 0 1 1 1 1 1 1 1 1 1 ...
- attr(*, "dimnames")=List of 2
..$ : chr [1:100] "A_1090" "B_1078" "A_1070" "B_1083" ...
..$ : chr [1:20003] "1.g" "1.a" "2.t" "2.g" ...
> class(phen)
[1] "factor"
> str(phen)
Factor w/ 2 levels "A","B": 2 2 2 2 1 2 1 1 2 2 ...
- attr(*, "names")= chr [1:100] "B_1132" "B_1105" "B_1079" "B_1133" ...
I'm not sure if this is an R issue or a treeWAS issue. Any feedback would be greatly appreciated!
Cheers =)