Re: [MORPHMET] Question on Morpho
Dear Pablo, Yes, I wrote a function to read all the landmark configurations stored in a folder and sampled in Amira... Try : require(devtools) install.github("evoshape/Arothron", local=FALSE) data=read.amira.dir("pathofyourfolder","auto") Best, Antonio Il 14 lug 2017 16:51, "Pablo Fisichella" ha scritto: Dear Antonio Excuse my delayed response. Thank you, your suggestion work quite fine, I made the GPA. However, as Paolo mentioned the coords looks so strange. I made .csv files from .pts files (in landmark). Here I attach one of my .csv and the corresponding .pts file. There is a way to read multiple .txt files from a folder similar to read.csv.folder?. which is the right format of the data frame? I have 3D coordinates derived from Amira and microscribe. Thank you 2017-07-13 5:30 GMT-03:00 Antonio Profico : > Dear Pablo, > > try: proc <- ProcGPA(a$arr). > > The object "a" is a list, the array is an element (arr) of the list. > > Let me know if works, > > Antonio > > 2017-07-13 4:59 GMT+01:00 Pablo Fisichella : > >> Dear All >> >> I am using the package Morpho in order to make some GM analyses but I >> obtained some errors. Now I am trying to use the function read.csv.folder >> >> I ran this command >> >> a<-read.csv.folder("E:/Documents/Zan2", x = 2:44, y = 2:4, rownames = >> NULL, header = TRUE, dec = ",", sep = ";", pattern = "csv", addSpec = NULL, >> back = TRUE) >> >> In the r console appeared my 5 files .csv (see below) with 43 >> landmarks/semilandmark in 3 dimensions for 5 individuals (0093, 0137, 0377, >> 1865, 1880) >> >> Then I tried perform a Procrustes Analysis using this function: proc <- >> ProcGPA(a) however this error appear: Error in addo(a3) : please provide >> 3D numeric array >> >> I feel that this error is related to my data array but I dont know how I >> can solve this in order to run a GPA >> >> thanks in advance, I appreciate any help you can give me. >> >> cheers, >> >> Pablo >> >> $arr >> , , 0093 >> >> X YZ >> 1 3.30e+08 -1.08e+08 5.22e+08 >> 2 2.94e+08 -5.35e+07 5.12e+08 >> 3 2.76e+08 -7.79e+07 5.59e+08 >> 4 3.13e+08 -1.37e+08 5.48e+08 >> 5 3.21e+08 -1.35e+08 5.45e+08 >> 6 3.27e+08 -1.33e+08 5.42e+08 >> 7 3.32e+08 -1.30e+08 5.39e+08 >> 8 3.36e+08 -1.28e+08 5.36e+08 >> 9 3.38e+08 -1.25e+08 5.34e+08 >> 10 3.39e+08 -1.23e+08 5.32e+08 >> 11 3.39e+08 -1.20e+08 5.30e+08 >> 12 3.37e+08 -1.17e+08 5.29e+08 >> 13 3.33e+08 -1.14e+08 5.27e+08 >> 14 3.29e+08 -9.43e+07 5.21e+08 >> 15 3.30e+08 -9.20e+07 5.19e+08 >> 16 3.31e+08 -8.92e+07 5.17e+08 >> 17 3.31e+08 -8.59e+07 5.16e+08 >> 18 3.30e+08 -8.20e+07 5.15e+08 >> 19 3.28e+08 -7.77e+07 5.14e+08 >> 20 3.24e+08 -7.33e+07 5.13e+08 >> 21 3.19e+08 -6.86e+07 5.11e+08 >> 22 3.12e+08 -6.36e+07 5.10e+08 >> 23 3.04e+08 -5.84e+07 5.08e+08 >> 24 2.84e+08 -4.42e+07 5.20e+08 >> 25 2.81e+08 -4.40e+07 5.27e+08 >> 26 2.79e+08 -4.37e+07 5.32e+08 >> 27 2.78e+08 -4.28e+07 5.38e+08 >> 28 2.77e+08 -4.35e+07 5.43e+08 >> 29 2.76e+08 -4.60e+07 5.47e+08 >> 30 2.76e+08 -4.93e+07 5.50e+08 >> 31 2.75e+08 -5.34e+07 5.53e+08 >> 32 2.76e+08 -5.83e+07 5.55e+08 >> 33 2.76e+08 -6.43e+07 5.56e+08 >> 34 2.83e+08 -8.41e+07 5.67e+08 >> 35 2.85e+08 -8.59e+07 5.69e+08 >> 36 2.87e+08 -8.77e+07 5.70e+08 >> 37 2.89e+08 -8.99e+07 5.72e+08 >> 38 2.91e+08 -9.27e+07 5.73e+08 >> 39 2.93e+08 -9.67e+07 5.74e+08 >> 40 2.96e+08 -1.02e+08 5.74e+08 >> 41 2.97e+08 -1.09e+08 5.73e+08 >> 42 2.97e+08 -1.17e+08 5.69e+08 >> 43 2.96e+08 -1.26e+08 5.65e+08 >> >> , , 0137 >> >> XYZ >> 1 5.31e+08 5660 1.45e+08 >> 2 5.64e+08 -330 1.47e+08 >> 3 5.85e+08 2190 1.03e+08 >> 4 5.35e+08 7310 9.86e+07 >> 5 5.35e+08 7580 1.07e+08 >> 6 5.34e+08 7780 1.14e+08 >> 7 5.35e+08 7940 1.20e+08 >> 8 5.34e+08 7990 1.26e+08 >> 9 5.34e+08 7910 1.31e+08 >> 10 5.33e+08 7750 1.35e+08 >> 11 5.33e+08 7490 1.38e+08 >> 12 5.31e+08 7090 1.40e+08 >> 13 5.30e+08 6620 1.41e+08 >> 14 5.39e+08 5470 1.50e+08 >> 15 5.42e+08 5320 1.53e+08 >> 16 5.44e+08 5160 1.55e+08 >> 17 5.47e+08 4940 1.57e+08 >> 18 5.49e+08 4610 1.59e+08 >> 19 5.52e+08 4120 1.61e+08 >> 20 5.54e+08 3400 1.61e+08 >> 21 5.56e+08 2490 1.59e+08 >> 22 5.57e+08 1450 1.57e+08 >> 23 5.58e+08 291 1.54e+08 >> 24 5.77e+08 -654 1.40e+08 >> 25 5.82e+08 -335 1.34e+08 >> 26 5.88e+08 632000 1.31e+08 >> 27 5.92e+08 417 1.27e+08 >> 28 5.94e+08 687 1.23e+08 >> 29 5.95e+08 931 1.19e+08 >> 30 5.94e+08 1150 1.16e+08 >> 31 5.93e+08 1340 1.13e+08 >> 32 5.92e+08 1530 1.11e+08 >> 33 5.90e+08 1730 1.09e+08 >> 34 5.85e+08 3280 9.65e+07 >> 35 5.85e+08 3610 9.61e+07 >> 36 5.84e+08 3930 9.58e+07 >> 37 5.84e+08 4270 9.54e+07 >> 38 5.83e+08 4650 9.45e+07 >> 39 5.81e+08 5040 9.30e+07 >> 40 5.76e+08 5420 9.10e+07 >> 41 5.70e+08 5830 8.85e+07 >> 42 5.62e+08 6080 8.60e+07 >> 43 5.52e+08 6310 8.31e+07 >> >> , , 0377 >> >> XYZ >> 1 5.
RE: [MORPHMET] Question on Morpho
Your values in csv are non-sensical. I am not aware of a single function that will read the pts files, but you can loop over the list of files and use read.table to read any ascii format. So something like this would work for your case pts.files=dir(patt='pts') temp=read.table(file=pts.files[1], skip=2, sep=' ', header = F) LM=array(dim=c(nrow(temp)), 3, length(pts.files)) for (i in 1:length) LM[,,i]=read.table(file=pts.files[i], skip=2, sep=" ", header=F)[,2:4] I didn’t check it, so you may need to edit a few things. M From: Pablo Fisichella [mailto:fisichellapa...@gmail.com] Sent: Friday, July 14, 2017 7:51 AM To: Antonio Profico Cc: MORPHMET ; morphmet_modera...@morphometrics.org Subject: Re: [MORPHMET] Question on Morpho Dear Antonio Excuse my delayed response. Thank you, your suggestion work quite fine, I made the GPA. However, as Paolo mentioned the coords looks so strange. I made .csv files from .pts files (in landmark). Here I attach one of my .csv and the corresponding .pts file. There is a way to read multiple .txt files from a folder similar to read.csv.folder?. which is the right format of the data frame? I have 3D coordinates derived from Amira and microscribe. Thank you 2017-07-13 5:30 GMT-03:00 Antonio Profico mailto:antonio.prof...@uniroma1.it>>: Dear Pablo, try: proc <- ProcGPA(a$arr). The object "a" is a list, the array is an element (arr) of the list. Let me know if works, Antonio 2017-07-13 4:59 GMT+01:00 Pablo Fisichella mailto:fisichellapa...@gmail.com>>: Dear All I am using the package Morpho in order to make some GM analyses but I obtained some errors. Now I am trying to use the function read.csv.folder I ran this command a<-read.csv.folder("E:/Documents/Zan2", x = 2:44, y = 2:4, rownames = NULL, header = TRUE, dec = ",", sep = ";", pattern = "csv", addSpec = NULL, back = TRUE) In the r console appeared my 5 files .csv (see below) with 43 landmarks/semilandmark in 3 dimensions for 5 individuals (0093, 0137, 0377, 1865, 1880) Then I tried perform a Procrustes Analysis using this function: proc <- ProcGPA(a) however this error appear: Error in addo(a3) : please provide 3D numeric array I feel that this error is related to my data array but I dont know how I can solve this in order to run a GPA thanks in advance, I appreciate any help you can give me. cheers, Pablo $arr , , 0093 X YZ 1 3.30e+08 -1.08e+08 5.22e+08 2 2.94e+08 -5.35e+07 5.12e+08 3 2.76e+08 -7.79e+07 5.59e+08 4 3.13e+08 -1.37e+08 5.48e+08 5 3.21e+08 -1.35e+08 5.45e+08 6 3.27e+08 -1.33e+08 5.42e+08 7 3.32e+08 -1.30e+08 5.39e+08 8 3.36e+08 -1.28e+08 5.36e+08 9 3.38e+08 -1.25e+08 5.34e+08 10 3.39e+08 -1.23e+08 5.32e+08 11 3.39e+08 -1.20e+08 5.30e+08 12 3.37e+08 -1.17e+08 5.29e+08 13 3.33e+08 -1.14e+08 5.27e+08 14 3.29e+08 -9.43e+07 5.21e+08 15 3.30e+08 -9.20e+07 5.19e+08 16 3.31e+08 -8.92e+07 5.17e+08 17 3.31e+08 -8.59e+07 5.16e+08 18 3.30e+08 -8.20e+07 5.15e+08 19 3.28e+08 -7.77e+07 5.14e+08 20 3.24e+08 -7.33e+07 5.13e+08 21 3.19e+08 -6.86e+07 5.11e+08 22 3.12e+08 -6.36e+07 5.10e+08 23 3.04e+08 -5.84e+07 5.08e+08 24 2.84e+08 -4.42e+07 5.20e+08 25 2.81e+08 -4.40e+07 5.27e+08 26 2.79e+08 -4.37e+07 5.32e+08 27 2.78e+08 -4.28e+07 5.38e+08 28 2.77e+08 -4.35e+07 5.43e+08 29 2.76e+08 -4.60e+07 5.47e+08 30 2.76e+08 -4.93e+07 5.50e+08 31 2.75e+08 -5.34e+07 5.53e+08 32 2.76e+08 -5.83e+07 5.55e+08 33 2.76e+08 -6.43e+07 5.56e+08 34 2.83e+08 -8.41e+07 5.67e+08 35 2.85e+08 -8.59e+07 5.69e+08 36 2.87e+08 -8.77e+07 5.70e+08 37 2.89e+08 -8.99e+07 5.72e+08 38 2.91e+08 -9.27e+07 5.73e+08 39 2.93e+08 -9.67e+07 5.74e+08 40 2.96e+08 -1.02e+08 5.74e+08 41 2.97e+08 -1.09e+08 5.73e+08 42 2.97e+08 -1.17e+08 5.69e+08 43 2.96e+08 -1.26e+08 5.65e+08 , , 0137 XYZ 1 5.31e+08 5660 1.45e+08 2 5.64e+08 -330 1.47e+08 3 5.85e+08 2190 1.03e+08 4 5.35e+08 7310 9.86e+07 5 5.35e+08 7580 1.07e+08 6 5.34e+08 7780 1.14e+08 7 5.35e+08 7940 1.20e+08 8 5.34e+08 7990 1.26e+08 9 5.34e+08 7910 1.31e+08 10 5.33e+08 7750 1.35e+08 11 5.33e+08 7490 1.38e+08 12 5.31e+08 7090 1.40e+08 13 5.30e+08 6620 1.41e+08 14 5.39e+08 5470 1.50e+08 15 5.42e+08 5320 1.53e+08 16 5.44e+08 5160 1.55e+08 17 5.47e+08 4940 1.57e+08 18 5.49e+08 4610 1.59e+08 19 5.52e+08 4120 1.61e+08 20 5.54e+08 3400 1.61e+08 21 5.56e+08 2490 1.59e+08 22 5.57e+08 1450 1.57e+08 23 5.58e+08 291 1.54e+08 24 5.77e+08 -654 1.40e+08 25 5.82e+08 -335 1.34e+08 26 5.88e+08 632000 1.31e+08 27 5.92e+08 417 1.27e+08 28 5.94e+08 687 1.23e+08 29 5.95e+08 931 1.19e+08 30 5.94e+08 1150 1.16e+08 31 5.93e+08 1340 1.13e+08 32 5.92e+08 1530 1.11e+08 33 5.90e+08 1730 1.09e+08 34 5.85e+08 3280 9.65e+07 35 5.85e+08 3610 9.61e+07 36 5.84e+08 3930 9.58e+07 37 5.84e+08 4270
Re: [MORPHMET] Question on Morpho
Dear Antonio Excuse my delayed response. Thank you, your suggestion work quite fine, I made the GPA. However, as Paolo mentioned the coords looks so strange. I made .csv files from .pts files (in landmark). Here I attach one of my .csv and the corresponding .pts file. There is a way to read multiple .txt files from a folder similar to read.csv.folder?. which is the right format of the data frame? I have 3D coordinates derived from Amira and microscribe. Thank you 2017-07-13 5:30 GMT-03:00 Antonio Profico : > Dear Pablo, > > try: proc <- ProcGPA(a$arr). > > The object "a" is a list, the array is an element (arr) of the list. > > Let me know if works, > > Antonio > > 2017-07-13 4:59 GMT+01:00 Pablo Fisichella : > >> Dear All >> >> I am using the package Morpho in order to make some GM analyses but I >> obtained some errors. Now I am trying to use the function read.csv.folder >> >> I ran this command >> >> a<-read.csv.folder("E:/Documents/Zan2", x = 2:44, y = 2:4, rownames = >> NULL, header = TRUE, dec = ",", sep = ";", pattern = "csv", addSpec = NULL, >> back = TRUE) >> >> In the r console appeared my 5 files .csv (see below) with 43 >> landmarks/semilandmark in 3 dimensions for 5 individuals (0093, 0137, 0377, >> 1865, 1880) >> >> Then I tried perform a Procrustes Analysis using this function: proc <- >> ProcGPA(a) however this error appear: Error in addo(a3) : please provide >> 3D numeric array >> >> I feel that this error is related to my data array but I dont know how I >> can solve this in order to run a GPA >> >> thanks in advance, I appreciate any help you can give me. >> >> cheers, >> >> Pablo >> >> $arr >> , , 0093 >> >> X YZ >> 1 3.30e+08 -1.08e+08 5.22e+08 >> 2 2.94e+08 -5.35e+07 5.12e+08 >> 3 2.76e+08 -7.79e+07 5.59e+08 >> 4 3.13e+08 -1.37e+08 5.48e+08 >> 5 3.21e+08 -1.35e+08 5.45e+08 >> 6 3.27e+08 -1.33e+08 5.42e+08 >> 7 3.32e+08 -1.30e+08 5.39e+08 >> 8 3.36e+08 -1.28e+08 5.36e+08 >> 9 3.38e+08 -1.25e+08 5.34e+08 >> 10 3.39e+08 -1.23e+08 5.32e+08 >> 11 3.39e+08 -1.20e+08 5.30e+08 >> 12 3.37e+08 -1.17e+08 5.29e+08 >> 13 3.33e+08 -1.14e+08 5.27e+08 >> 14 3.29e+08 -9.43e+07 5.21e+08 >> 15 3.30e+08 -9.20e+07 5.19e+08 >> 16 3.31e+08 -8.92e+07 5.17e+08 >> 17 3.31e+08 -8.59e+07 5.16e+08 >> 18 3.30e+08 -8.20e+07 5.15e+08 >> 19 3.28e+08 -7.77e+07 5.14e+08 >> 20 3.24e+08 -7.33e+07 5.13e+08 >> 21 3.19e+08 -6.86e+07 5.11e+08 >> 22 3.12e+08 -6.36e+07 5.10e+08 >> 23 3.04e+08 -5.84e+07 5.08e+08 >> 24 2.84e+08 -4.42e+07 5.20e+08 >> 25 2.81e+08 -4.40e+07 5.27e+08 >> 26 2.79e+08 -4.37e+07 5.32e+08 >> 27 2.78e+08 -4.28e+07 5.38e+08 >> 28 2.77e+08 -4.35e+07 5.43e+08 >> 29 2.76e+08 -4.60e+07 5.47e+08 >> 30 2.76e+08 -4.93e+07 5.50e+08 >> 31 2.75e+08 -5.34e+07 5.53e+08 >> 32 2.76e+08 -5.83e+07 5.55e+08 >> 33 2.76e+08 -6.43e+07 5.56e+08 >> 34 2.83e+08 -8.41e+07 5.67e+08 >> 35 2.85e+08 -8.59e+07 5.69e+08 >> 36 2.87e+08 -8.77e+07 5.70e+08 >> 37 2.89e+08 -8.99e+07 5.72e+08 >> 38 2.91e+08 -9.27e+07 5.73e+08 >> 39 2.93e+08 -9.67e+07 5.74e+08 >> 40 2.96e+08 -1.02e+08 5.74e+08 >> 41 2.97e+08 -1.09e+08 5.73e+08 >> 42 2.97e+08 -1.17e+08 5.69e+08 >> 43 2.96e+08 -1.26e+08 5.65e+08 >> >> , , 0137 >> >> XYZ >> 1 5.31e+08 5660 1.45e+08 >> 2 5.64e+08 -330 1.47e+08 >> 3 5.85e+08 2190 1.03e+08 >> 4 5.35e+08 7310 9.86e+07 >> 5 5.35e+08 7580 1.07e+08 >> 6 5.34e+08 7780 1.14e+08 >> 7 5.35e+08 7940 1.20e+08 >> 8 5.34e+08 7990 1.26e+08 >> 9 5.34e+08 7910 1.31e+08 >> 10 5.33e+08 7750 1.35e+08 >> 11 5.33e+08 7490 1.38e+08 >> 12 5.31e+08 7090 1.40e+08 >> 13 5.30e+08 6620 1.41e+08 >> 14 5.39e+08 5470 1.50e+08 >> 15 5.42e+08 5320 1.53e+08 >> 16 5.44e+08 5160 1.55e+08 >> 17 5.47e+08 4940 1.57e+08 >> 18 5.49e+08 4610 1.59e+08 >> 19 5.52e+08 4120 1.61e+08 >> 20 5.54e+08 3400 1.61e+08 >> 21 5.56e+08 2490 1.59e+08 >> 22 5.57e+08 1450 1.57e+08 >> 23 5.58e+08 291 1.54e+08 >> 24 5.77e+08 -654 1.40e+08 >> 25 5.82e+08 -335 1.34e+08 >> 26 5.88e+08 632000 1.31e+08 >> 27 5.92e+08 417 1.27e+08 >> 28 5.94e+08 687 1.23e+08 >> 29 5.95e+08 931 1.19e+08 >> 30 5.94e+08 1150 1.16e+08 >> 31 5.93e+08 1340 1.13e+08 >> 32 5.92e+08 1530 1.11e+08 >> 33 5.90e+08 1730 1.09e+08 >> 34 5.85e+08 3280 9.65e+07 >> 35 5.85e+08 3610 9.61e+07 >> 36 5.84e+08 3930 9.58e+07 >> 37 5.84e+08 4270 9.54e+07 >> 38 5.83e+08 4650 9.45e+07 >> 39 5.81e+08 5040 9.30e+07 >> 40 5.76e+08 5420 9.10e+07 >> 41 5.70e+08 5830 8.85e+07 >> 42 5.62e+08 6080 8.60e+07 >> 43 5.52e+08 6310 8.31e+07 >> >> , , 0377 >> >> XYZ >> 1 5.48e+08 7410 1.08e+08 >> 2 5.79e+08 1400 1.01e+08 >> 3 6.04e+08 4710 5.61e+07 >> 4 5.56e+08 9770 5.76e+07 >> 5 5.57e+08 9850 6.97e+07 >> 6 5.56e+08 9740 8.03e+07 >> 7 5.55e+08 9630 8.91e+07 >> 8 5.54e+08 9440 9.59e+07 >> 9 5.52e+08 9190 1.01e+08 >> 10 5.50e+08 8850 1.0
Re: [MORPHMET] Question on Morpho
I think the problem is in the input: 3.30e+08 -1.08e+08 5.22e+08 etc. are very strange coordinates; I suspect some trouble in the csv file with comma symbol (, or .) ... maybe you prepared it with excel? There could be the problem. best paolio 2017-07-13 10:30 GMT+02:00 Antonio Profico : > Dear Pablo, > > try: proc <- ProcGPA(a$arr). > > The object "a" is a list, the array is an element (arr) of the list. > > Let me know if works, > > Antonio > > 2017-07-13 4:59 GMT+01:00 Pablo Fisichella : > >> Dear All >> >> I am using the package Morpho in order to make some GM analyses but I >> obtained some errors. Now I am trying to use the function read.csv.folder >> >> I ran this command >> >> a<-read.csv.folder("E:/Documents/Zan2", x = 2:44, y = 2:4, rownames = >> NULL, header = TRUE, dec = ",", sep = ";", pattern = "csv", addSpec = NULL, >> back = TRUE) >> >> In the r console appeared my 5 files .csv (see below) with 43 >> landmarks/semilandmark in 3 dimensions for 5 individuals (0093, 0137, 0377, >> 1865, 1880) >> >> Then I tried perform a Procrustes Analysis using this function: proc <- >> ProcGPA(a) however this error appear: Error in addo(a3) : please provide >> 3D numeric array >> >> I feel that this error is related to my data array but I dont know how I >> can solve this in order to run a GPA >> >> thanks in advance, I appreciate any help you can give me. >> >> cheers, >> >> Pablo >> >> $arr >> , , 0093 >> >> X YZ >> 1 3.30e+08 -1.08e+08 5.22e+08 >> 2 2.94e+08 -5.35e+07 5.12e+08 >> 3 2.76e+08 -7.79e+07 5.59e+08 >> 4 3.13e+08 -1.37e+08 5.48e+08 >> 5 3.21e+08 -1.35e+08 5.45e+08 >> 6 3.27e+08 -1.33e+08 5.42e+08 >> 7 3.32e+08 -1.30e+08 5.39e+08 >> 8 3.36e+08 -1.28e+08 5.36e+08 >> 9 3.38e+08 -1.25e+08 5.34e+08 >> 10 3.39e+08 -1.23e+08 5.32e+08 >> 11 3.39e+08 -1.20e+08 5.30e+08 >> 12 3.37e+08 -1.17e+08 5.29e+08 >> 13 3.33e+08 -1.14e+08 5.27e+08 >> 14 3.29e+08 -9.43e+07 5.21e+08 >> 15 3.30e+08 -9.20e+07 5.19e+08 >> 16 3.31e+08 -8.92e+07 5.17e+08 >> 17 3.31e+08 -8.59e+07 5.16e+08 >> 18 3.30e+08 -8.20e+07 5.15e+08 >> 19 3.28e+08 -7.77e+07 5.14e+08 >> 20 3.24e+08 -7.33e+07 5.13e+08 >> 21 3.19e+08 -6.86e+07 5.11e+08 >> 22 3.12e+08 -6.36e+07 5.10e+08 >> 23 3.04e+08 -5.84e+07 5.08e+08 >> 24 2.84e+08 -4.42e+07 5.20e+08 >> 25 2.81e+08 -4.40e+07 5.27e+08 >> 26 2.79e+08 -4.37e+07 5.32e+08 >> 27 2.78e+08 -4.28e+07 5.38e+08 >> 28 2.77e+08 -4.35e+07 5.43e+08 >> 29 2.76e+08 -4.60e+07 5.47e+08 >> 30 2.76e+08 -4.93e+07 5.50e+08 >> 31 2.75e+08 -5.34e+07 5.53e+08 >> 32 2.76e+08 -5.83e+07 5.55e+08 >> 33 2.76e+08 -6.43e+07 5.56e+08 >> 34 2.83e+08 -8.41e+07 5.67e+08 >> 35 2.85e+08 -8.59e+07 5.69e+08 >> 36 2.87e+08 -8.77e+07 5.70e+08 >> 37 2.89e+08 -8.99e+07 5.72e+08 >> 38 2.91e+08 -9.27e+07 5.73e+08 >> 39 2.93e+08 -9.67e+07 5.74e+08 >> 40 2.96e+08 -1.02e+08 5.74e+08 >> 41 2.97e+08 -1.09e+08 5.73e+08 >> 42 2.97e+08 -1.17e+08 5.69e+08 >> 43 2.96e+08 -1.26e+08 5.65e+08 >> >> , , 0137 >> >> XYZ >> 1 5.31e+08 5660 1.45e+08 >> 2 5.64e+08 -330 1.47e+08 >> 3 5.85e+08 2190 1.03e+08 >> 4 5.35e+08 7310 9.86e+07 >> 5 5.35e+08 7580 1.07e+08 >> 6 5.34e+08 7780 1.14e+08 >> 7 5.35e+08 7940 1.20e+08 >> 8 5.34e+08 7990 1.26e+08 >> 9 5.34e+08 7910 1.31e+08 >> 10 5.33e+08 7750 1.35e+08 >> 11 5.33e+08 7490 1.38e+08 >> 12 5.31e+08 7090 1.40e+08 >> 13 5.30e+08 6620 1.41e+08 >> 14 5.39e+08 5470 1.50e+08 >> 15 5.42e+08 5320 1.53e+08 >> 16 5.44e+08 5160 1.55e+08 >> 17 5.47e+08 4940 1.57e+08 >> 18 5.49e+08 4610 1.59e+08 >> 19 5.52e+08 4120 1.61e+08 >> 20 5.54e+08 3400 1.61e+08 >> 21 5.56e+08 2490 1.59e+08 >> 22 5.57e+08 1450 1.57e+08 >> 23 5.58e+08 291 1.54e+08 >> 24 5.77e+08 -654 1.40e+08 >> 25 5.82e+08 -335 1.34e+08 >> 26 5.88e+08 632000 1.31e+08 >> 27 5.92e+08 417 1.27e+08 >> 28 5.94e+08 687 1.23e+08 >> 29 5.95e+08 931 1.19e+08 >> 30 5.94e+08 1150 1.16e+08 >> 31 5.93e+08 1340 1.13e+08 >> 32 5.92e+08 1530 1.11e+08 >> 33 5.90e+08 1730 1.09e+08 >> 34 5.85e+08 3280 9.65e+07 >> 35 5.85e+08 3610 9.61e+07 >> 36 5.84e+08 3930 9.58e+07 >> 37 5.84e+08 4270 9.54e+07 >> 38 5.83e+08 4650 9.45e+07 >> 39 5.81e+08 5040 9.30e+07 >> 40 5.76e+08 5420 9.10e+07 >> 41 5.70e+08 5830 8.85e+07 >> 42 5.62e+08 6080 8.60e+07 >> 43 5.52e+08 6310 8.31e+07 >> >> , , 0377 >> >> XYZ >> 1 5.48e+08 7410 1.08e+08 >> 2 5.79e+08 1400 1.01e+08 >> 3 6.04e+08 4710 5.61e+07 >> 4 5.56e+08 9770 5.76e+07 >> 5 5.57e+08 9850 6.97e+07 >> 6 5.56e+08 9740 8.03e+07 >> 7 5.55e+08 9630 8.91e+07 >> 8 5.54e+08 9440 9.59e+07 >> 9 5.52e+08 9190 1.01e+08 >> 10 5.50e+08 8850 1.03e+08 >> 11 5.49e+08 8580 1.04e+08 >> 12 5.49e+08 8350 1.04e+08 >> 13 5.49e+08 8060 1.04e+08 >> 14 5.56e+08 6830 1.14e+08 >> 15 5.58e+08 6600 1.16e+08 >> 16 5.60e+08 6360 1.17e+08 >
Re: [MORPHMET] Question on Morpho
Dear Pablo, try: proc <- ProcGPA(a$arr). The object "a" is a list, the array is an element (arr) of the list. Let me know if works, Antonio 2017-07-13 4:59 GMT+01:00 Pablo Fisichella : > Dear All > > I am using the package Morpho in order to make some GM analyses but I > obtained some errors. Now I am trying to use the function read.csv.folder > > I ran this command > > a<-read.csv.folder("E:/Documents/Zan2", x = 2:44, y = 2:4, rownames = > NULL, header = TRUE, dec = ",", sep = ";", pattern = "csv", addSpec = NULL, > back = TRUE) > > In the r console appeared my 5 files .csv (see below) with 43 > landmarks/semilandmark in 3 dimensions for 5 individuals (0093, 0137, 0377, > 1865, 1880) > > Then I tried perform a Procrustes Analysis using this function: proc <- > ProcGPA(a) however this error appear: Error in addo(a3) : please provide > 3D numeric array > > I feel that this error is related to my data array but I dont know how I > can solve this in order to run a GPA > > thanks in advance, I appreciate any help you can give me. > > cheers, > > Pablo > > $arr > , , 0093 > > X YZ > 1 3.30e+08 -1.08e+08 5.22e+08 > 2 2.94e+08 -5.35e+07 5.12e+08 > 3 2.76e+08 -7.79e+07 5.59e+08 > 4 3.13e+08 -1.37e+08 5.48e+08 > 5 3.21e+08 -1.35e+08 5.45e+08 > 6 3.27e+08 -1.33e+08 5.42e+08 > 7 3.32e+08 -1.30e+08 5.39e+08 > 8 3.36e+08 -1.28e+08 5.36e+08 > 9 3.38e+08 -1.25e+08 5.34e+08 > 10 3.39e+08 -1.23e+08 5.32e+08 > 11 3.39e+08 -1.20e+08 5.30e+08 > 12 3.37e+08 -1.17e+08 5.29e+08 > 13 3.33e+08 -1.14e+08 5.27e+08 > 14 3.29e+08 -9.43e+07 5.21e+08 > 15 3.30e+08 -9.20e+07 5.19e+08 > 16 3.31e+08 -8.92e+07 5.17e+08 > 17 3.31e+08 -8.59e+07 5.16e+08 > 18 3.30e+08 -8.20e+07 5.15e+08 > 19 3.28e+08 -7.77e+07 5.14e+08 > 20 3.24e+08 -7.33e+07 5.13e+08 > 21 3.19e+08 -6.86e+07 5.11e+08 > 22 3.12e+08 -6.36e+07 5.10e+08 > 23 3.04e+08 -5.84e+07 5.08e+08 > 24 2.84e+08 -4.42e+07 5.20e+08 > 25 2.81e+08 -4.40e+07 5.27e+08 > 26 2.79e+08 -4.37e+07 5.32e+08 > 27 2.78e+08 -4.28e+07 5.38e+08 > 28 2.77e+08 -4.35e+07 5.43e+08 > 29 2.76e+08 -4.60e+07 5.47e+08 > 30 2.76e+08 -4.93e+07 5.50e+08 > 31 2.75e+08 -5.34e+07 5.53e+08 > 32 2.76e+08 -5.83e+07 5.55e+08 > 33 2.76e+08 -6.43e+07 5.56e+08 > 34 2.83e+08 -8.41e+07 5.67e+08 > 35 2.85e+08 -8.59e+07 5.69e+08 > 36 2.87e+08 -8.77e+07 5.70e+08 > 37 2.89e+08 -8.99e+07 5.72e+08 > 38 2.91e+08 -9.27e+07 5.73e+08 > 39 2.93e+08 -9.67e+07 5.74e+08 > 40 2.96e+08 -1.02e+08 5.74e+08 > 41 2.97e+08 -1.09e+08 5.73e+08 > 42 2.97e+08 -1.17e+08 5.69e+08 > 43 2.96e+08 -1.26e+08 5.65e+08 > > , , 0137 > > XYZ > 1 5.31e+08 5660 1.45e+08 > 2 5.64e+08 -330 1.47e+08 > 3 5.85e+08 2190 1.03e+08 > 4 5.35e+08 7310 9.86e+07 > 5 5.35e+08 7580 1.07e+08 > 6 5.34e+08 7780 1.14e+08 > 7 5.35e+08 7940 1.20e+08 > 8 5.34e+08 7990 1.26e+08 > 9 5.34e+08 7910 1.31e+08 > 10 5.33e+08 7750 1.35e+08 > 11 5.33e+08 7490 1.38e+08 > 12 5.31e+08 7090 1.40e+08 > 13 5.30e+08 6620 1.41e+08 > 14 5.39e+08 5470 1.50e+08 > 15 5.42e+08 5320 1.53e+08 > 16 5.44e+08 5160 1.55e+08 > 17 5.47e+08 4940 1.57e+08 > 18 5.49e+08 4610 1.59e+08 > 19 5.52e+08 4120 1.61e+08 > 20 5.54e+08 3400 1.61e+08 > 21 5.56e+08 2490 1.59e+08 > 22 5.57e+08 1450 1.57e+08 > 23 5.58e+08 291 1.54e+08 > 24 5.77e+08 -654 1.40e+08 > 25 5.82e+08 -335 1.34e+08 > 26 5.88e+08 632000 1.31e+08 > 27 5.92e+08 417 1.27e+08 > 28 5.94e+08 687 1.23e+08 > 29 5.95e+08 931 1.19e+08 > 30 5.94e+08 1150 1.16e+08 > 31 5.93e+08 1340 1.13e+08 > 32 5.92e+08 1530 1.11e+08 > 33 5.90e+08 1730 1.09e+08 > 34 5.85e+08 3280 9.65e+07 > 35 5.85e+08 3610 9.61e+07 > 36 5.84e+08 3930 9.58e+07 > 37 5.84e+08 4270 9.54e+07 > 38 5.83e+08 4650 9.45e+07 > 39 5.81e+08 5040 9.30e+07 > 40 5.76e+08 5420 9.10e+07 > 41 5.70e+08 5830 8.85e+07 > 42 5.62e+08 6080 8.60e+07 > 43 5.52e+08 6310 8.31e+07 > > , , 0377 > > XYZ > 1 5.48e+08 7410 1.08e+08 > 2 5.79e+08 1400 1.01e+08 > 3 6.04e+08 4710 5.61e+07 > 4 5.56e+08 9770 5.76e+07 > 5 5.57e+08 9850 6.97e+07 > 6 5.56e+08 9740 8.03e+07 > 7 5.55e+08 9630 8.91e+07 > 8 5.54e+08 9440 9.59e+07 > 9 5.52e+08 9190 1.01e+08 > 10 5.50e+08 8850 1.03e+08 > 11 5.49e+08 8580 1.04e+08 > 12 5.49e+08 8350 1.04e+08 > 13 5.49e+08 8060 1.04e+08 > 14 5.56e+08 6830 1.14e+08 > 15 5.58e+08 6600 1.16e+08 > 16 5.60e+08 6360 1.17e+08 > 17 5.63e+08 6120 1.19e+08 > 18 5.66e+08 5820 1.21e+08 > 19 5.69e+08 5410 1.22e+08 > 20 5.71e+08 4800 1.22e+08 > 21 5.72e+08 3960 1.20e+08 > 22 5.73e+08 2960 1.19e+08 > 23 5.71e+08 2000 1.13e+08 > 24 5.93e+08 1110 9.49e+07 > 25 5.99e+08 1380 9.20e+07 > 26 6.04e+08 1650 8.87e+07 > 27 6.08e+08 1930 8.54e+07 > 28 6.11e+08 2230 8.19e+07 > 29 6.13e+08 2530 7.83e+07 > 30 6.14e+08 2830 7.44e+07 > 31 6.13e+08 3140 7