Re: [R] pvmnorm, error message

From: Prof Brian Ripley <ripley_at_stats.ox.ac.uk>
Date: Wed, 09 May 2007 15:05:09 +0100 (BST)

But it should not give that error (once pmvnorm is spelled correctly I get no error).

I believe the message indicates that the installed package is corrupt.

On Wed, 9 May 2007, Dimitris Rizopoulos wrote:

> A is not correlation matrix; try this instead:
>
> A <- diag(rep(0.5, 3))
> A[1, 2] <- 0.5
> A[1, 3] <- 0.25
> A[2, 3] <- 0.5
> A <- A + t(A)
> pmvnorm(lower = rep(-Inf, 3), upper = rep(2, 3), corr = A)
>
>
> Best,
> Dimitris
>
> ----
> Dimitris Rizopoulos
> Ph.D. Student
> Biostatistical Centre
> School of Public Health
> Catholic University of Leuven
>
> Address: Kapucijnenvoer 35, Leuven, Belgium
> Tel: +32/(0)16/336899
> Fax: +32/(0)16/337015
> Web: http://med.kuleuven.be/biostat/
> http://www.student.kuleuven.be/~m0390867/dimitris.htm
>
>
> ----- Original Message -----
> From: "Andreas Faller" <afaller__at_web.de>
> To: <r-help_at_stat.math.ethz.ch>
> Sent: Wednesday, May 09, 2007 3:16 PM
> Subject: [R] pvmnorm, error message
>
>
> Hello there!
>
> My operating system is Windows XP, my version of R is the latest
> (R-2.5.0). Recently I have downloaded the package "mvtnorm" and a
> problem with the command "pmvnorm" occured. Trying to enter the lines
> ...
>
> A <- diag(3)
> A[1,2] <-0.5
> A[1,3] <- 0.25
> A[2,3] <- 0.5
> pvmnorm(lower=c(-Inf,-Inf,-Inf), upper=c(2,2,2),mean = c(0,0,0),
> corr=A)
>
> I got the following error message:
>
> .Fortran("mvtdst", N = as.integer(n), NU=as.integer(df), lower =
> as.double(lower), :
> Fortran Symbolname "mvtdst" nicht in der DLL für Paket "mvtnorm"
>
> Can anyone advise what to do now to get rid of this problem? Thank you
> very much indeed.
>
> Regards, Andreas Faller
>
> ______________________________________________

> R-help_at_stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>
>
> Disclaimer: http://www.kuleuven.be/cwis/email_disclaimer.htm
>
> ______________________________________________
> R-help_at_stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>

-- 
Brian D. Ripley,                  ripley_at_stats.ox.ac.uk
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272866 (PA)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595


______________________________________________ R-help_at_stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

Received on Wed 09 May 2007 - 14:22:11 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Wed 09 May 2007 - 15:01:48 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.