[R] function gdist, dist and vegdist in mvpart

Sina Muster sina.muster at web.de
Wed Mar 8 19:12:11 CET 2006


Dear R community,

I am analyzing plant communities with the function mvpart, using a 
dissimilarit matrix as input. The matrix is calculated with the funtion 
gdist.

fit <- mvpart(gdist (ba12[,18:29], meth="maximum", full=TRUE, 
sq=F) ~ 			beers + slope_dem + elev_dem+ plc_dem + pr_curv+ 			 
+curv+max_depth+doc_rocks+ abandon+land_use+ca_old,
	data=ba12, xv="p")

This works fine. Now I would like to use other dissimilarity measures as 
can be found in the function dist (STATS) or vegdist (VEGAN).

De'Ath notes that gdist should be interchangeable with dist - but I 
receive following error message:

 > fit21 <- mvpart(dist (ba12[,18:29], meth="minkowski", diag=T, 
upper=T 		p=2) ~
	beers + slope_dem + elev_dem+ plc_dem + pr_curv+ 			 
curv+max_depth+doc_rocks+ abandon+land_use+ca_old,
	data=ba12, xv="p")
Error in model.frame(formula, rownames, variables, varnames, extras, 
extranames,  :
        Variable lengths differ


When computing the matrices directly...

ba12.dist1 <- gdist (ba12[,18:29], method="maximum", full=T, sq=F)
ba12.dist2 <- vegdist (ba12[,18:29], method="jaccard", diag=T, upper=T)
ba12.dist3 <- dist (ba12[,18:29], method="minkowski", p=1, diag=T, upper=T)

and looking at them, there appears to be no difference though only 
ba12.dist1 works with mvpart. With the other two I get the same error 
message as above.

What I am missing? Where is the problem?

Thanks for helping.
Sina

University of Potsdam, Germany
Department of Geoecology




More information about the R-help mailing list