----- Forwarded message from Jason Hemingway
<[email protected]> -----
Date: Fri, 15 Nov 2013
06:03:15 -0500
From: Jason Hemingway
<[email protected]>
Reply-To: Jason Hemingway
<[email protected]>
Subject: RE: 3D confidence
ellipsoids in R
To: "[email protected]"
<[email protected]>
Hi Andrea
You can try the ellipse3d function in the rgl package. I used a for loop to draw each ellipse.
http://cran.r-project.org/web/packages/rgl/
open3d(FOV=1)
mesh<-mshape(X)
rgl.lines(mesh[t(polylines),1], mesh[t(polylines),2], mesh[t(polylines),3], col="black", alpha=c(1), line_antialias = TRUE)
for (i in 1:dim(X[,,1])[1]) {
plot3d(ellipse3d(cov(t(X[i,,])), centre=apply(X[i,,], 1, mean)), col="green", alpha=0.95, specular="green", shininess=70.0, add = TRUE)
}
Regards,
Jason
Jason Hemingway
Lecturer (Morphological Anatomy) | School of Anatomical Sciences | University of the Witwatersrand | Johannesburg, South Africa
Tel: (+27)11 717 2114 | Fax: (+27)11 717 2422 | Email: [email protected]
-----Original Message-----
From: [email protected] [mailto:[email protected]]
Sent: 15 November 2013 12:08 PM
To: [email protected]
Subject: 3D confidence ellipsoids in R
----- Forwarded message from andrea cardini <[email protected]> -----
Date: Wed, 13 Nov 2013 08:14:14 -0500
From: andrea cardini <[email protected]>
Reply-To: andrea cardini <[email protected]>
Subject: 3D confidence ellipsoids in R
To: [email protected]
Dear All,
I'd greatly appreciate any suggestion about how to draw 3D confidence ellipsoids in R (packages? scripts? examples?)?
Thanks a lot in advance for your help.
Cheers
Andrea
Dr. Andrea Cardini
Researcher in Animal Biology, Dipartimento di Scienze Chimiche e Geologiche, Università di Modena e Reggio Emilia, l.go S. Eufemia 19, 41121 Modena, Italy Honorary Fellow, Centre for Anatomical and Human Sciences, University of Hull, Cottingham Road, Hull, HU6 7RX, UK & University of York, Heslington, York YO10 5DD, UK Adjunct Associate Professor, Centre for Forensic Science , The University of Western Australia, 35 Stirling Highway, Crawley WA 6009, Australia
E-mail address: [email protected], [email protected]
WEBPAGE: http://sites.google.com/site/hymsfme/drandreacardini
Summary of research interests at:
http://www.dscg.unimore.it/site/home/ricerca/aree-di-ricerca/evolution-taxonomy-and-forensics.html
FREE Yellow BOOK on Geometric Morphometrics:
http://www.italian-journal-of-mammalogy.it/issue/view/405
or full volume at:
http://www.italian-journal-of-mammalogy.it/public/journals/3/issue_241_complete_100.pdf
Editorial board for:
Zoomorphology:
http://www.springer.com/life+sciences/animal+sciences/journal/435
Journal of Zoological Systematics and
Evolutionary Research: http://www.wiley.com/bw/journal.asp?ref=0947-5745&site=1
Hystrix, the Italian Journal of
Mammalogy: http://www.italian-journal-of-mammalogy.it/
----- End forwarded message -----
<table width="100%" border="0" cellspacing="0" cellpadding="0" style="width:100%;">
<tr>
<td align="left" style="text-align:justify;"><font face="arial,sans-serif" size="1" color="#999999"><span style="font-size:11px;">This communication is intended for the addressee only. It is confidential. If you have received this communication in error, please notify us immediately and destroy the original message. You may not copy or disseminate this communication without the permission of the University. Only authorised signatories are competent to enter into agreements on behalf of the University and recipients are thus advised that the content of this message may not be legally binding on the University and may contain the personal views and opinions of the author, which are not necessarily the views and opinions of The University of the Witwatersrand, Johannesburg. All agreements between the University and outsiders are subject to South African Law unless the University agrees in writing to the contrary. </span></font></td>
</tr>
</table>
You can try the ellipse3d function in the rgl package. I used a for loop to draw each ellipse.
http://cran.r-project.org/web/packages/rgl/
open3d(FOV=1)
mesh<-mshape(X)
rgl.lines(mesh[t(polylines),1], mesh[t(polylines),2], mesh[t(polylines),3], col="black", alpha=c(1), line_antialias = TRUE)
for (i in 1:dim(X[,,1])[1]) {
plot3d(ellipse3d(cov(t(X[i,,])), centre=apply(X[i,,], 1, mean)), col="green", alpha=0.95, specular="green", shininess=70.0, add = TRUE)
}
Regards,
Jason
Jason Hemingway
Lecturer (Morphological Anatomy) | School of Anatomical Sciences | University of the Witwatersrand | Johannesburg, South Africa
Tel: (+27)11 717 2114 | Fax: (+27)11 717 2422 | Email: [email protected]
-----Original Message-----
From: [email protected] [mailto:[email protected]]
Sent: 15 November 2013 12:08 PM
To: [email protected]
Subject: 3D confidence ellipsoids in R
----- Forwarded message from andrea cardini <[email protected]> -----
Date: Wed, 13 Nov 2013 08:14:14 -0500
From: andrea cardini <[email protected]>
Reply-To: andrea cardini <[email protected]>
Subject: 3D confidence ellipsoids in R
To: [email protected]
Dear All,
I'd greatly appreciate any suggestion about how to draw 3D confidence ellipsoids in R (packages? scripts? examples?)?
Thanks a lot in advance for your help.
Cheers
Andrea
Dr. Andrea Cardini
Researcher in Animal Biology, Dipartimento di Scienze Chimiche e Geologiche, Università di Modena e Reggio Emilia, l.go S. Eufemia 19, 41121 Modena, Italy Honorary Fellow, Centre for Anatomical and Human Sciences, University of Hull, Cottingham Road, Hull, HU6 7RX, UK & University of York, Heslington, York YO10 5DD, UK Adjunct Associate Professor, Centre for Forensic Science , The University of Western Australia, 35 Stirling Highway, Crawley WA 6009, Australia
E-mail address: [email protected], [email protected]
WEBPAGE: http://sites.google.com/site/hymsfme/drandreacardini
Summary of research interests at:
http://www.dscg.unimore.it/site/home/ricerca/aree-di-ricerca/evolution-taxonomy-and-forensics.html
FREE Yellow BOOK on Geometric Morphometrics:
http://www.italian-journal-of-mammalogy.it/issue/view/405
or full volume at:
http://www.italian-journal-of-mammalogy.it/public/journals/3/issue_241_complete_100.pdf
Editorial board for:
Zoomorphology:
http://www.springer.com/life+sciences/animal+sciences/journal/435
Journal of Zoological Systematics and
Evolutionary Research: http://www.wiley.com/bw/journal.asp?ref=0947-5745&site=1
Hystrix, the Italian Journal of
Mammalogy: http://www.italian-journal-of-mammalogy.it/
----- End forwarded message -----
<table width="100%" border="0" cellspacing="0" cellpadding="0" style="width:100%;">
<tr>
<td align="left" style="text-align:justify;"><font face="arial,sans-serif" size="1" color="#999999"><span style="font-size:11px;">This communication is intended for the addressee only. It is confidential. If you have received this communication in error, please notify us immediately and destroy the original message. You may not copy or disseminate this communication without the permission of the University. Only authorised signatories are competent to enter into agreements on behalf of the University and recipients are thus advised that the content of this message may not be legally binding on the University and may contain the personal views and opinions of the author, which are not necessarily the views and opinions of The University of the Witwatersrand, Johannesburg. All agreements between the University and outsiders are subject to South African Law unless the University agrees in writing to the contrary. </span></font></td>
</tr>
</table>
----- End forwarded message
-----
