similar to: mvtnorm package

Dear all, I'm using the ellipse package and I would like to verify if the confidence region that I build with this package can be equivalent to an Union Intersection Test. I used different value for the t-statistic but I can not find the right equivalence. Does someone know how to choose the right value? Thanks a lot. Antonio -- Antonio Lucadamo, Dipartimento di Scienze Economiche e

Dear all, I need to use the "qmvnorm" function in mtvnorm package. Here is the error message I got > qmvnorm(0.05, tail="upper", sigma=var)$quantile Error in uniroot(pfct, interval = interval) : f() values at end points not of opposite sign There is no problem for 50th quantile. > qmvnorm(0.5, tail="upper", sigma=var)$quantile [1] -0.5325527 Can anyone

Dear all, I wrote the following function previously. It worked fine with the old mvtnorm package. Somehow with the updated package, I got a error message when trying to use the function. I need some help. It is sort of urgent. Can anyone please take a look. The function is the following. Thank you very much! Hannah library(mvtnorm) cc_f <- function(m, rho,

Dear all, I have the following program for a multiple comparison procedure. There are two functions for the two steps. First step is to calculate the critical values, while the second step is the actual procedure [see below: program with two functions]. This work fine. However, However I want to put them into one function for the convenience of later use [see below: program with one

Dear all, I am trying to calculate certain critical values from bivariate normal distribution (please see the function below). m <- 10 rho <- 0.1 k <- 2 alpha <- 0.05 ## calculate critical constants cc_z <- numeric(m) var <- matrix(c(1,rho,rho,1), nrow=2, ncol=2, byrow=T) for (i in 1:m){ if (i <= k) {cc_z[i] <- qmvnorm((k*(k-1))/(m*(m-1))*alpha,

Hello I have read that under Solaris one can use DTrace to get I/O request size distribution on a global scale (also on a per process/pid basis). See for example http://prefetch.net/articles/observeiodtk.html Can anyone?recommend?an alternative to get similar information under CentOS? I looked into dtrace for linux but it seems still work in progress, even putting aside CDDL issues ...

I have these two dataframes in which 'id' is the key field > tabella id nome 1 1 PIEMONTE 2 2 VALLED'AOSTA 3 3 LOMBARDIA 4 4 TRENTINO 5 5 VENETO 6 6 FRIULI AND > tab id nome 1 1 PIEMONTE 2 2 VALLED'AOSTA 3 3 LOMBARDIA 4 4 TRENTINO 5 25 CAMPANIA 6 28 LAZIO Is there any

Dear list, I tried to use the nparcomp to run some post hoc non-parametric comparison and got and error. Error in uniroot(pfct, interval = interval) : f() values at end points not of opposite sign Appreciate any comments. the command line: >nparcomp(Ulceration~Group,data=test,type='Dunnett',control='Non-treated') Jun

Sorry to be pedantic but you are a little careless about version numbers .... > It was actually Debian "testing" that has R 1.4.0. Debian "stable" has R 0.9 or something really old > like that. 1.4.1 last I checked, and that is in the Debian mirrors (where testing is frozen, hence the older version). On CRAN it should be 1.5.1, built on my testing system. >

Hi all, I've been trying to install the "mvtnorm" package (in a Linux R version) without sucess. I write install.packages("mvtnorm",lib="/home/posmae/cnaber",repos="http://cran.uk.r-project.org/") and the following message arises ================================================== downloaded 160Kb * Installing *source* package 'mvtnorm'

Hi all, I know pnorm() is able to calculate a sequence of probabilities by providing a vector of means, and a single number for location and standard deviation. For example, pnorm(0.5,c(0.5,2),1) will produce [1] 0.5000000 0.0668072. However, I wonder if its multivariate counterpart mvtnorm can do the same thing? Namely, if I provide a sequence of (vector-)means and a single vector for

Jonathan, > I tried to re-intall R from Debian "stable", in which R's version is 1.4.0. And the installation > of "mvtnorm" works. Debian's own "stable" R is way older; this must have been a user-contributed one. And I suspect that it was build without readline support (see below). > I then re-installed R yet again from Debian

Hi all, I tried to install the mvtnorm library on linux (mandrake version) and i obtained the following result : /usr/bin/ld: cannot find -lreadline collect2: ld returned 1 exit status make: *** [mvtnorm.so] Erreur 1 ERROR: compilation failed for package 'mvtnorm' Could you help me ? Best regards, Olivier

Hi all! My apologies for posting such a naive question. I've been trying to run multiple comparison contrasts on several GLM models. For doing so I've tried to use the mulcomp CRAN package. However, before I can make the package to work I have to load the mvtnorm package which I cannot find in CRAN (or anywhere else). Any inputs? Thanks a lot, /J Jose A. Andres. Post-doc Associate.

Hi R-users, I just started learning R. I have a project on lot quality assurance sampling (LQAS). In this project I have to develop LQAS plans to make decision on stopping / continuing a programme. The LQAS plans is based on cluster sampling: selection of k clusters (villages) of m children each for a total sample of n. In this connection I would like to estimate the classification errors ( alpha

Hi R-users, I just started learning R. I have a project on lot quality assurance sampling (LQAS). In this project I have to develop LQAS plans to make decision on stopping / continuing a programme. The LQAS plans is based on cluster sampling: selection of k clusters (villages) of m children each for a total sample of n. In this connection I would like to estimate the classification errors ( alpha

Dear All, For function: rmvnorm{mvtnorm} in (library mvtnorm, not splus2R), if I generate 2 bivariate normal samples as follows: > rmvnorm(2,mean=rep(0,2),sigma=diag(2)) [,1] [,2] [1,] 2.0749459 1.4932752 [2,] -0.9886333 0.3832266 Where is the first sample, it is stored in the first row or the first column? Does this function store samples row-wise or column-wise? Thank

Hi Dirk, It was actually Debian "testing" that has R 1.4.0. Debian "stable" has R 0.9 or something really old like that. Anyways, I did what A.J. Rossini suggested: # apt-get install libreadline4-dev Then I installed R-1.5.1 from Debian "unstable". The installation of "mvtnorm" has no trouble after this. I think that you are probably right: the

Announcement: mvtnorm Multivariate Normal and T Distribution mvtnorm implements two R functions for the computation of the multivariate t and normal distribution: pmvt: Computes the the distribution function of the multivariate t distribution for arbitary limits, degrees of freedom and correlation matrices based on algorithms by Genz and Bretz. pmvnorm: Computes the distribution

Announcement: mvtnorm Multivariate Normal and T Distribution mvtnorm implements two R functions for the computation of the multivariate t and normal distribution: pmvt: Computes the the distribution function of the multivariate t distribution for arbitary limits, degrees of freedom and correlation matrices based on algorithms by Genz and Bretz. pmvnorm: Computes the distribution