Displaying 2 results from an estimated 2 matches for "areab".
Did you mean:
area
2007 Jul 07
2
No convergence using ADAPT
...ry short time. The program is below:
## START PROGRAM
trial <- function(input)
{
pmvnorm(lower = c(0,0), upper = c(2, 2), mean = input, sigma = matrix(c(.1, 0,
0, .1), nrow = 2, ncol = 2, byrow = FALSE))
}
require(mvtnorm)
require(adapt)
bottomB <- -5*sqrt(.1)
topB <- 2 + 5*sqrt(.1)
areaB <- (topB - bottomB)^2
unscaled.Po.in.a <- adapt(2, lo = c(bottomB, bottomB), up = c(topB, topB),
minpts = 1000, eps = 1e-4, functn = trial)
(1/areaB)*unscaled.Po.in.a$value
## FINISH PROGRAM
I tried to run the program again changing a.) sigma in the trial function, b.)
upper in the tria...
2011 Mar 20
2
R as a non-functional language
I'm reading Torgo (2010) *Data Mining with
R*<http://www.liaad.up.pt/~ltorgo/DataMiningWithR/code.html>in
preparation for a class I'll be teaching next quarter. Here's an
example
that is very non-functional.
> pH <- c(4.5,7,7.3,8.2,6.3)
> names(pH) <- c('area1','area2','mud','dam','middle')
> pH
area1 area2 mud dam