Jeff,
The dataset are 3D landmarks taken from a skull (Homo being the dataset of 3
species of the genus Homo) I am trying to perform Generalised Procrustes
Distance and superimposition and as you can see, I am told x is not True.
If I also perform is.numeric(Homo) I am told that it is FALSE. I hope that
is enough context for you. I also hope you know what's going wrong in my
dataset.
> Homo=read.table("D:/Homo.txt", sep="")
> Homo
X1x X1y X1z X2X X2Y X2Z X3X
X3Y
SK 125.2288 -207.2243 95.13449 94.08087 -184.7764 98.85405 159.1837
-189.3104
GIB 121.4254 -228.1540 33.61633 84.07050 -209.2446 52.82806 159.5556
-205.1629
KBW 128.5904 -232.1195 47.11266 92.96498 -215.2787 48.39832 164.3043
-213.7220
X3Z X4X X4Y X4Z X5X X5Y X5Z
X6x
SK 94.84497 129.3257 -192.9469 120.64075 102.09383 -167.0360 140.4830
161.4940
GIB 48.10354 123.3033 -235.6632 68.96723 82.34183 -220.0837 104.7756
164.1860
KBW 51.10052 126.2948 -223.1627 81.66055 75.98707 -201.2843 108.0087
177.3304
X6y X6z X7x X7y X7z X8x X8y
X8z
SK -169.4289 136.2563 56.07269 -156.3827 89.49998 195.9156 -165.8967
79.33151
GIB -221.0951 101.2779 57.06448 -194.0515 61.04704 188.6181 -184.5234
57.21327
KBW -199.6966 110.2110 55.50871 -172.5727 55.57302 201.8666 -173.5339
57.47801
X9x X9y X9z X10x X10y X10z
SK 71.01711 -158.7208 126.92698 192.5974 -158.2639 116.59634
GIB 57.25282 -198.6531 93.25098 186.9386 -195.8407 88.64017
KBW 55.54756 -180.8983 92.43064 197.7266 -180.1343 96.22781
> Homo=arrayspecs(Homo, 10, 3, byLand=F)
> Homo
, , SK
[,1] [,2] [,3]
[1,] 125.2288 -207.2243 95.13449
[2,] 94.08087 -184.7764 98.85405
[3,] 159.1837 -189.3104 94.84497
[4,] 129.3257 -192.9469 120.6407
[5,] 102.0938 -167.036 140.483
[6,] 161.494 -169.4289 136.2563
[7,] 56.07269 -156.3827 89.49998
[8,] 195.9156 -165.8967 79.33151
[9,] 71.01711 -158.7208 126.927
[10,] 192.5974 -158.2639 116.5963
, , GIB
[,1] [,2] [,3]
[1,] 121.4254 -228.154 33.61633
[2,] 84.0705 -209.2446 52.82806
[3,] 159.5556 -205.1629 48.10354
[4,] 123.3033 -235.6632 68.96723
[5,] 82.34183 -220.0837 104.7756
[6,] 164.186 -221.0951 101.2779
[7,] 57.06448 -194.0515 61.04704
[8,] 188.6181 -184.5234 57.21327
[9,] 57.25282 -198.6531 93.25098
[10,] 186.9386 -195.8407 88.64017
, , KBW
[,1] [,2] [,3]
[1,] 128.5904 -232.1195 47.11266
[2,] 92.96498 -215.2787 48.39832
[3,] 164.3043 -213.722 51.10052
[4,] 126.2948 -223.1627 81.66055
[5,] 75.98707 -201.2843 108.0087
[6,] 177.3304 -199.6966 110.211
[7,] 55.50871 -172.5727 55.57302
[8,] 201.8666 -173.5339 57.47801
[9,] 55.54756 -180.8983 92.43064
[10,] 197.7266 -180.1343 96.22781
> Homo=gpagen(Homo, Proj=T, ProcD=T,ShowPlot=T,curves=NULL, surfaces=NULL,
> pointscale=1)
Error: is.atomic(x) is not TRUE
>
--
View this message in context:
http://r.789695.n4.nabble.com/is-numeric-FALSE-tp4666530p4666564.html
Sent from the R help mailing list archive at Nabble.com.
______________________________________________
R-help@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.
