similar to: How to specify Variance Covariance matrix of residuals?

Dear R team, I would like to do a multivariate meta-analysis in R using the nlme package. In meta-analysis I fix the residuals to known sampling errors. As I do a multivariate analysis, I have a variance-covariance matrix of sampling errors. Unfortunately, via varFixed I can only fix a vector of sampling errors and no matrix. In the R package metafor using the rma.mv function I can insert the

> Date: Mon, 20 Feb 2006 16:43:55 -0600 > From: Aldi Kraja <aldi at wustl.edu> > > Hi, > Using package gclus in R, I have created some graphs that show the > trends within subgroups of data and correlations among 9 variables (v1-v9). > Being interested for more details on these data I have produced also the > var-covar matrices. > Question: From a pair of two

I'm using the lme function in nmle to estimate the variance components of a fully nested two-level model: Y_ijk = mu + a_i + b_j(i) + e_k(j(i)) lme computes estimates of the variances for a, b, and e, call them v_a, v_b, and v_e, and I can use the intervals function to get confidence intervals. My understanding is that these intervals are probably not that robust plus I need intervals on the

Dear all, I am trying to simulate data sets from a model fitted with glmmPQL, in order to compute the distribution of a summary statistics. My data are binomial and I have a correlation term in my model. My model is structured in the following way m <- glmmPQL( fixed = cbind(sucess,failure) ~ x1 + x2 + ... , random = ~ 1 | bidon, correlation = corGaus(form=~

Dear R users, i have a matrix with 31 rows and 444 columns and i want to extract every 37th column of that matrix starting from 1. more precisely i want to select columns 1, 38,75, 112 and so on. then doing the same by starting from column number 2(2,39,76,113.......). i know that there is a manual way of doing it but i wanted to make it more quickly as i have fairly large data to dealth with.

I am pleased to announce the immediate availability of Maven in version 3.3 on CentOS Linux 7 x86_64, delivered via a Software Collection (SCL) built by the SCLo Special Interest Group (https://wiki.centos.org/SpecialInterestGroup/SCLo). QuickStart ---------- You can get started in three easy steps: $ sudo yum install centos-release-scl $ sudo yum install rh-maven33 $ scl

I am pleased to announce the immediate availability of Maven in version 3.3 on CentOS Linux 7 x86_64, delivered via a Software Collection (SCL) built by the SCLo Special Interest Group (https://wiki.centos.org/SpecialInterestGroup/SCLo). QuickStart ---------- You can get started in three easy steps: $ sudo yum install centos-release-scl $ sudo yum install rh-maven33 $ scl

I am pleased to announce the immediate availability of Maven in version 3.5 on CentOS Linux 7 x86_64, delivered via a Software Collection (SCL) built by the SCLo Special Interest Group (https://wiki.centos.org/SpecialInterestGroup/SCLo). QuickStart ---------- You can get started in three easy steps: $ sudo yum install centos-release-scl $ sudo yum install rh-maven35 $ scl enable

On 25.11.2014 15:02, net.rsync at io7m.com wrote: > 'Lo. > > I've run into a frustrating issue when trying to synchronize a > directory hierarchy over a reliable (but slow) connection to an > unreliable remote. Basically, I have the following: > > http://mvn-repository.io7m.com/com/io7m/ > > This is a set of nested directories containing binaries and

(until a better answer comes along) The killed rsync process should leave a temporary file .<file-name>.<random> If you rename the temporary to the real file name, rsync should continue from about where it left off. -----Original Message----- From: rsync-bounces at lists.samba.org [mailto:rsync-bounces at lists.samba.org] On Behalf Of net.rsync at io7m.com Sent: Tuesday, November 25,

Dear all, I have a multivariate dataset containing 100,000 or more points. I want find the p-value for the dataset of points coming from a particular multivariate normal distribution With mean vector u Covariance matrix s2 So H0: points ~ MVN( u, s2) H1: points not ~ MVN( u, s2) How do I find the p-value in R? To me this is a likelihood ratio test problem. In H0 the parameters are

Hi Folks, Given k RVs with MVN distribution N(mu,S) (S a kxk covariance matrix), let (w.l.o.g.) X1 denote the first r of them, and X2 the last (k-r). Likewise, let mu1 and mu2 denote their respective expectations. Then, of course, the expectation of X2 given X1=x1 is mu2 + S21*inv(S22)*(x1 - mu1) and the covariance matrix of X2 given X1=x2 is S22 - S21*inv(X11)*S12 where Sij is the

Hi All, I need to fit a mutlivariate normal model to a dataset in order to obtain the mean and covariance parameters. I see that the MVN function in the MCLUST package can do this, however when I try to run even the simplest example provided in the documentation, as below, I get the following error: n <- 1000 set.seed(0) x <- rnorm(n, mean = -1, sd = 2) mvn(modelName = "X", x)

'Lo. I've run into a frustrating issue when trying to synchronize a directory hierarchy over a reliable (but slow) connection to an unreliable remote. Basically, I have the following: http://mvn-repository.io7m.com/com/io7m/ This is a set of nested directories containing binaries and sources for projects I develop/maintain. Every time a new release is made, I deploy the binaries and

<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN"> <html> <head> <meta content="text/html;charset=UTF-8" http-equiv="Content-Type"> <title></title> </head> <body bgcolor="#ffffff" text="#000000"> <font face="Helvetica, Arial, sans-serif">Hi,<br> <br> I am

Hi, For some 64-bit immediates, ARM can generate a two instruction sequence. For example, for the following code: ; 734439407618 = 0x000000ab00000002 define i64 @f1(i64 %a) { %tmp = sub i64 %a, 734439407618 ret i64 %tmp } We should be able to generate: subs r0, r0, #1 sbc r1, r1, #171 @ 0xab But instead we get: $ llvm-as < %s | llc -march=thumb -mattr=+thumb2 mvn r2, #1

Hi R-experts, I have a problem using nlme. I use the following code to group my data: Parameterg <- groupedData( result ~ time | Batch, data = Batchdata, labels = list( x = "Time", y = "analysis") ) and then uses the nlme function to fit a nonlinear mixed model that includes Process as a fixed covariate: nlme.model001epr <- nlme(result ~ A0 * exp(- ( exp(A1)

Dear R users, I'm applying a correlation structure in a mixed model (nmle function) to control for spatial correlation between land parcels that are adjacent to each other. I generated X,Y coordinates in ArcGIS for each land parcel and used them in the correlation form like this: test.exp<-corExp(1, form = ~ X + Y) test.exp<- Initialize(test.exp,dataset) However, the correlation

Hi All: Not sure why my previous question never got posted. Here I am seeking some help on my code. I am using the following code to run MCMC simulation on the following data using the model below: # Data matrix<-NULL > csvs<-paste("MVN", 1:2,".csv",sep="") > for (i in 1:length(csvs)){ + matrix[[i]]<- read.csv(file=csvs[i],header=TRUE) +

Dear list Having just started to use the Windows version, I am very impressed with it's package handling as well as the gui. So I tried to see what was due for update and packages such as Hmisc, Matrix and others came up. But when I had updated them - which took a few goes as something hung between here and Bristol - I noticed that the default packages such as nmle, MASS had disappeared. I