search for: miwa

Currently doing some work on local maxima on a random field and have encountered an issue with the Miwa algorithm used with the pmvnorm function in the mvtnorm R package. Based on recommendations by Mi et al., we ran the mvtnorm package using the Miwa algorithm, since we have a maximum of 4 dimensions with non-singular matrices. However, running the estimation procedure in this way, we obtained inco...

...---------------- > *From:* R-help <r-help-bounces at r-project.org> on behalf of Eric Berger < > ericjberger at gmail.com> > *Sent:* 02 October 2017 16:16:13 > *To:* Bert Gunter > *Cc:* r-help at r-project.org; Hollie Johnson (PGR) > *Subject:* Re: [R] Issues with 'Miwa' algorithm in mvtnorm package > > Hi Hollie, > I tried to reproduce your example but could not. Here is what I did. Please > explain if you did something different. > > x <- c(9.358*10^-3, -8.165*10^-3, -1.689*10^-8, > -8.165*10^-3, 9.358*10^-3, 1.854*10^-8, &...

...o it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) On Mon, Oct 2, 2017 at 5:16 AM, Hollie Johnson (PGR) < h.a.johnson at newcastle.ac.uk> wrote: > Currently doing some work on local maxima on a random field and have > encountered an issue with the Miwa algorithm used with the pmvnorm function > in the mvtnorm R package. > > Based on recommendations by Mi et al., we ran the mvtnorm package using > the Miwa algorithm, since we have a maximum of 4 dimensions with > non-singular matrices. However, running the estimation procedure in th...

...TRUE) Regards, Hollie ________________________________ From: R-help <r-help-bounces at r-project.org> on behalf of Eric Berger <ericjberger at gmail.com> Sent: 02 October 2017 16:16:13 To: Bert Gunter Cc: r-help at r-project.org; Hollie Johnson (PGR) Subject: Re: [R] Issues with 'Miwa' algorithm in mvtnorm package Hi Hollie, I tried to reproduce your example but could not. Here is what I did. Please explain if you did something different. x <- c(9.358*10^-3, -8.165*10^-3, -1.689*10^-8, -8.165*10^-3, 9.358*10^-3, 1.854*10^-8, -1.689*10^-8, 1.854*10^-8...

Hello. So as my english is not so good, i'll try to write my story in bash commands :o) miwa@zork3$ uname -mrs FreeBSD 5.5-STABLE i386 miwa@zork3$ su root@zork3# cd /usr/ports/sysutils/nat root@zork3# make && make install root@zork3# cd ../nut-usb root@zork3# make && make install root@zork3# cd /usr/local/etc/nut root@zork3# cat ups.conf user=root [MGE-Ellipse] driver=newhi...

.../ash/libs ** R ** preparing package for lazy loading ** help *** installing help indices ** building package indices ... ** testing if installed package can be loaded * DONE (ash) * installing *source* package ‘mvtnorm’ ... ** libs gcc -I/usr/share/R/include -fpic -std=gnu99 -O3 -pipe -g -c miwa.c -o miwa.o gfortran -fpic -O3 -pipe -g -c mvt.f -o mvt.o gcc -I/usr/share/R/include -fpic -std=gnu99 -O3 -pipe -g -c randomF77.c -o randomF77.o gfortran -fpic -O3 -pipe -g -c tvpack.f -o tvpack.o gfortran -fpic -O3 -pipe -g -c tvpackAPI.f -o tvpackAPI.o gcc -I/usr/share/R/includ...