marginals3d: marginals3d
In twDEMC: parallel DEMC

Description Usage Arguments Value Author(s) See Also Examples

integrate over all others but 3 columns

1 2	marginals3d(x, vars = c(2, 3, 4), intCol = 1, n = round(nrow(x)^(1/3)/12), dimnames = .colNames(x, vars), fAgg = mean, ...)

`x`	numeric matrix
`vars`	index (integer or string): column to define integration boxes
`intCol`	index (integer or string): column to average
`n`	numeric vector: number of groups per variable
`dimnames`	dimension names of the result array
`fAgg`	function applied over the subsets of x[,intCol]
`...`	further arguments to fAgg

numeric array with each item representing application of fAgg to x[,intCol] and dimnames equal to means of the classes

Thomas Wutzler

marginals1d twDEMCBlockInt

.tmp.f <- function(){
	res <- marginals3d(sample,vars=c("kY","kO","tLagLeaf"),n=8)
	#library(Rcmdr)
	ds <- melt(res)
	scatter3d(ds$kY, ds$tLagLeaf, ds$kO
		, surface=FALSE
		,bg="white", axis.scales=TRUE, grid=TRUE, ellipsoid=FALSE
		, xlab="kY",ylab="tLagLeaf", zlab="kO"
		, point.col=rev(heat.colors(100))[round(rescale(ds$value,to=c(1,100)))]
		,sphere.size=1.5
	)
}