similar to: under what conditions would rtmvnorm (from package tmvtnorm) produce all NaNs

Hello! I was wondering if it's possible to see the actual code of checkTmvArgs function that is part of the code for rtmvnorm (which is below - I just typed "rtmvnorm" on the prompt). I get an error: Error in checkTmvArgs(mean, sigma, lower, upper) : sigma must be a symmetric matrix At the same time I am pretty sure that the matrix I am passing as sigma is a var-covar matrix

Dear R users, the tmvtnorm package, the package for the truncated multivariate normal and Student-t distribution, has been updated on CRAN. The major changes in version 1.0-2 (2010-03-04) are: * The package now provides methods for the truncated multivariate Student-t distribution, i.e. random number generation, density function, distribution functions like rtmvt(), dtmvt() und ptmvt() and

Dear R users, the tmvtnorm package, the package for the truncated multivariate normal and Student-t distribution, has been updated on CRAN. The major changes in version 1.0-2 (2010-03-04) are: * The package now provides methods for the truncated multivariate Student-t distribution, i.e. random number generation, density function, distribution functions like rtmvt(), dtmvt() und ptmvt() and

I want to use the rtmvt from the {tmvtnorm} package using the "gibbs" algorithm but how to i specify the nested function rtmvnorm to use gibbs as well? Right now I am using the code: for (i in 1:g){ for (j in 1:n){ sgamma[,,i,j] = rtmvt(n=50, mean=mu[i,j], sigma[i,j], df=nu[i], lower=rep(0,2),algorithm="gibbs") } } heres an example of one iteration: >

When sampling from a multi-variate truncated normal using rtmvnorm from the tmvnorm-package, I experience extreme performance differences between two of my computers. On my laptop computer, draws take ~5s, on my desktop ~30s. I need to run MCMCs with repeat calls to rtmvnorm on my desktop. The code in both cases is exactly the same. rtmvnorm seems not to have any machine dependent defaults

Hi, dear all, I wish to get one sample (2500-d vector) from the truncated multivariate normal distribution, so I choose use the R function rtmvnorm() to do this. But the error information shows that for this function, the dimension should be lower than 1000, So could you help me to find out if there's any solution could do such sampling from a truncated multivariate normal distribution

Dear R users, the tmvtnorm package, the package for the truncated multivariate normal distribution, has been updated on CRAN. The major changes in version 0.8 are: * Reimplemented the Gibbs' sampler for random number generation in Fortran for performance reasons. This compiled code is now even faster than conventional rejection sampling and allows for the generation of large amounts of

Dear R users, the tmvtnorm package, the package for the truncated multivariate normal distribution, has been updated on CRAN. The major changes in version 0.8 are: * Reimplemented the Gibbs' sampler for random number generation in Fortran for performance reasons. This compiled code is now even faster than conventional rejection sampling and allows for the generation of large amounts of

I am sampling from the truncated multivariate student t distribution "rtmvt" in the package {tmvtnorm}. My question is about the mean vector. Is it possible to define a mean vector outside of the truncated region? Thank you in advance for any help. -- View this message in context: http://r.789695.n4.nabble.com/truncated-distributions-tp3422245p3422245.html Sent from the R help mailing

I have been using the rtmvt function in the {tmvtnorm} package i'm getting the warning: "Acceptance rate is very low and rejection sampling becomes inefficient. Consider using Gibbs sampling." but i AM specifying the gibbs algorithm!!: rtmvt(M, mean=q[,,i,j], sigma=((u[i,j] + nu[i])/(p+nu[i]))*delta[,,i], df=ceiling(nu[i]+p), lower=c(0,0), algorithm="gibbs") Any

Dears Samba users I have a network configurated with a some samba servers e all the clients are Windows (95, 98, Me and NT). For some rasons my access to samba servers is with plaintextpassord and not for encryptedpassword. I need in a moment connect som machines with Windows XP. I read the documentation e follow the steps to configurations, but is not work. I can?t configure

>>>>> David Winsemius <dwinsemius at comcast.net> >>>>> on Tue, 9 May 2017 14:33:04 -0700 writes: >> On May 9, 2017, at 2:05 PM, Czarek Kowalski <czarek230800 at gmail.com> wrote: >> >> I have already posted that in attachement - pdf file. > I see that now. I failed to scroll to the 3rd page. from a late reader:

Dear Members, I am working with 6-dimensional Student-t distribution with 4 degrees of freedom truncated to [20; 60]. I have generated 100 000 samples from truncated multivariate Student-t distribution using rtmvt function from package ?tmvtnorm?. I have also calculated mean vector using equation (3) from attached pdf. The problem is, that after summing all elements in one column of rtmvt result

Hi, all I didn't get any response from swig for my question. see if I can get some help here Thanks ---------- Forwarded message ---------- From: charlie <charlie.xia.fdu@gmail.com> Date: Tue, Nov 4, 2008 at 1:55 PM Subject: SWIG with R and C++ STL To: swig-user@lists.sourceforge.net Hi all, I am new to SWIG. I encountered some problem when I try to SWIG to R some C++ modules. Here

Dear all, I have a large vector (lets call it myVector) and I want to plot its value with the logic below yaxis<-myVector[1] yaxis<-c(xaxis,mean(myvector[2:3]) yaxis<-c(xaxis,mean(myvector[4:8]) yaxis<c(xaxis,mean(myvector[9:16]) yaxis<c(xaxis,mean(myvector[17:32]) this has to stop when the new ..... yaxis<c(xaxis,mean(myvector[1024:2048]) will not find the correspondent number

I have been trying to write a new class that reimplements the vector class. As a test of my overloading I decided to try and and call t.test on two vectors which were objects of my class rather than the default class. The overloaded length I wrote seems to work correctly. I used setMethod("length", "myvector", ...) setMethod("mean", "myvector", ...)

Hello, I have an expression that I want to substitute for something else. myvector<-c("ajkhfkiuwe","Variable(kg/min)") if I use the following code to extract "variable(kg/min)" and substitute it for "va" gsub(myvector[grep("var", myvector, ignore=T)], "va", myvector) grep identifies the element of the vector that matches my

Hi Eric, Thanks for having a look into this. I think you have a small typo... B <- matrix(x, nrow=3, byrow = TRUE) should read B <- matrix(y, nrow=3, byrow = 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

>From: Johannes Graumann <johannes_graumann at web.de> >Date: 2007/12/18 Tue PM 04:40:37 CST >To: r-help at stat.math.ethz.ch >Subject: [R] All anchored series from a vector? lapply(1:length(myvector) function(.length) { c(myvector[1}:myvector[.length]) }) but test it because i didn't. >Hi all, > >What may be a smart, efficient way to get the following result:

Dear R users, I'm having trouble with calculating pvalues for my 2d dataset. First I performed clustering and I would like to get some info about the strength of cluster membership for each point. I've calculated (thanks to nice people help) the multivariate normal densities (mnd) using dmvnorm function: