Dear Ari,
you can try this code:
library(geomorph)
library(Morpho)
X=rnorm(100)
Y=rnorm(100)
Z=rnorm(100)
XYZ=cbind(X,Y,Z)
Array=(array(cbind(X,Y,Z),dim=c(30,3,10)))
#dimnames(Array)[[3]]=list(c("specimen1")) if n=1 the array k*d*n
dimnames(Array)[[3]]=c(paste("specimen",c(1:10),sep="_"))
writeland.tps(A=Array, "myarray.tps", scale = NULL)
write.pts(x=Array, filename = "myarray", rownames =
c(paste("specimen",c(1:10),sep="_")))
probably the error is in the imported data. Do you convert your Landmark sets
in a matrix or (if you have more specimens) in an array?
try this code:
is.matrix(yourmatrix)
is.array(youarrray)
is.numeric(yourmatrix)
is.numeric(yourarray)
______________________________________
Antonio Profico
PhD student
Department of Environmental Biology – Dipartimento di Biologia Ambientale
SAPIENZA Università di Roma
Lab. 06 4991 2690
Mob. 3293440766
Da: Ariadne Schulz
Data invio: venerdì 22 gennaio 2016 13:47
A: MORPHMET
Hello all,
I have another hopefully simple geomorph question. How does one write a .pts
file? I've found how to write a .tps file, and I have managed to create a .pts
file, but it was empty as far as I could tell. (I used writeland.tps and for
the filename simply appended .pts, but I'm pretty sure that's not how it's
done) Is there a simple way of doing this?
Best,
Ari
--
MORPHMET may be accessed via its webpage at http://www.morphometrics.org
---
You received this message because you are subscribed to the Google Groups
"MORPHMET" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
--
MORPHMET may be accessed via its webpage at http://www.morphometrics.org
---
You received this message because you are subscribed to the Google Groups
"MORPHMET" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
