search for: dpomatrix

After loading the package "Matrix" (version 0.995-12), using the summary function with an lme (package nlme version 3.1-75) object results in an error message saying Fehler in dim(x) : kein Slot des Namens "Dim" f?r dieses Objekt der Klasse "correlation" (translated: 'Error in dim(x) : no slot of the name "Dim" for this object of class

...a bit of help extracting the residual error variance from the VarCorr structure from lmer. #Here's a 2-way random effects model lmer.1 <- lmer(rating ~ (1|person)+(1|rater), data = dat) #Get the structure vc.fit <- VarCorr(lmer.1) #results in..... $person 1 x 1 Matrix of class "dpoMatrix" (Intercept) (Intercept) 0.7755392 $rater 1 x 1 Matrix of class "dpoMatrix" (Intercept) (Intercept) 0.2054469 attr(,"sc") [1] 0.5051518 #I can pull out the person and rater variance components easy enough. For example... vc.person <- vc.fit...

Hi: I have some problems when I use the function "solve" function in a loop. In the following code, I have a diagonal martix "ttt" whose elements change in every iteration in a loop. I defined a "dpoMatrix"class before the loop so I do not need to define this class every time in the loop. The reason is to save some computing time. The code is below. The inverse of the matrix "ttt" should change according to the change of "ttt" in the loop. However, the values in "sigma...

Dear all, I try to run the example of lmer and get the following error message. > library(lme4) > example(lmer) lmer> (fm1 <- lmer(Reaction ~ Days + (Days | Subject), sleepstudy)) [[1]] Error in get(x, envir, mode, inherits) : variable "as.dpoMatrix" was not found This error message is similar to what I get with other models. It looks like the mer class has a slightly different structure. Anybody an idea how to solve this? I am using R 2.4.1 under linux and the latest releases of lme4 and Matrix lme4_0.9975-10 Matrix_0.9975-8 > ve...

Hi, I have been having a tedious issue with lmer models with lots of factors and lots of levels. In order to get the basic information at the beginning of the print out I also have to generate these enormous tables as well. Is there a method command to leave off all of the effects and correlations? Or, do I have to go to string commands?

hello all, i'm trying to use QPmat, from the popbio package. it appears to be based on solve.QP and is intended for making a population projection matrix. QPmat asks for: nout, A time series of population vectors and C, C constraint matrix, (with two more vectors, b and nonzero). i believe the relevant code from QPmat is: function (nout, C, b, nonzero) { if (!"quadprog" %in%

...ing, but what about the level-1 variance, viz., the 3.215072 term? (actually this term squared) Didn't see anything in the archives on this. Cheers, David > fm <- lmer( dv ~ time.num*drug + (1 | Patient.new), data=dat.new ) > VarCorr(fm) $Patient.new 1 x 1 Matrix of class "dpoMatrix" (Intercept) (Intercept) 8.519916 attr(,"sc") scale 3.215072 > VarCorr(fm)[[1]][1] [1] 8.519916 > VarCorr(fm)[[2]][1] Error in VarCorr(fm)[[2]] : subscript out of bounds ########################################################## set.seed(500) n.time...

...ne until the last step. make[1]: Leaving directory `/tmp/R.INSTALL.TW3399/Matrix/src/AMD' pgCC -L/usr/lib64 -L/usr/X11R6/lib64 -pgf90libs -o Matrix.so CHMfactor.o Csparse.o TMatrix_as.o Tsparse.o init.o Mutils.o chm_common.o cs.o cs_utils.o dense.o dgCMatrix.o dgTMatrix.o dgeMatrix.o dpoMatrix.o dppMatrix.o dsCMatrix.o dsyMatrix.o dspMatrix.o dtCMatrix.o dtTMatrix.o dtrMatrix.o dtpMatrix.o factorizations.o ldense.o lgCMatrix.o sparseQR.o SPQR.a CHOLMOD.a COLAMD.a AMD.a -L/usr/local/apps/R/R-2.8.1/lib64/R/lib -lRlapack -lRblas -lpgf90 -lpgf90_rpm1 -lpgf902 -lpgf90rtl -lpgftnrtl -lnspg...

...lled Matrix. The compilation failed with the following messages, ....... CC -library=stlport4 -G -L/opt/csw/lib/sparcv9 -L/opt/solstudio12.2/prod/lib/v9 -o Matrix.so CHMfactor.o Csparse.o TMatrix_as.o Tsparse.o init.o Mutils.o chm_common.o cs.o cs_utils.o dense.o dgCMatrix.o dgTMatrix.o dgeMatrix.o dpoMatrix.o dppMatrix.o dsCMatrix.o dsyMatrix.o dspMatrix.o dtCMatrix.o dtTMatrix.o dtrMatrix.o dtpMatrix.o factorizations.o ldense.o lgCMatrix.o sparseQR.o abIndex.o CHOLMOD.a COLAMD.a AMD.a -L/apps/sparcv9/R-2.12.0/lib/R/lib -lRlapack -L/apps/sparcv9/R-2.12.0/lib/R/lib -lRblas -lifai -lsunimath -lfai -lfai...

Dear all, I'm trying to store/extract the mean& standard error of the fixed effects parameter and the variance of the random effects parameter from "lmer" procedure from mlmre4 package developed by bates n pinheiro. while storing fixed effects parameter is straight forward, the same is not true for storing the variance parameter of the random effects. kindly help me ~prabhu

...rror occurs during building of the recommended packages and also if Matrix is compiled separately: [...] CC -G -lCstd -L/opt/sun/sunstudio12/lib/amd64 -o Matrix.so CHMfactor.o Csparse.o TMatrix_as.o Tsparse.o init.o Mutils.o chm_common.o cs.o cs_utils.o dense.o dgCMatrix.o dgTMatrix.o dgeMatrix.o dpoMatrix.o dppMatrix.o dsCMatrix.o dsyMatrix.o dspMatrix.o dtCMatrix.o dtTMatrix.o dtrMatrix.o dtpMatrix.o factorizations.o ldense.o lgCMatrix.o sparseQR.o CHOLMOD.a COLAMD.a AMD.a -L/home/user/R/R-devel/lib -lRlapack -L/home/user/R/R-devel/lib -lRblas -R/opt/sun/sunstudio12/lib/amd64:/opt/sun/sunstudio12/l...

...nclude -I/usr/include -I/usr/sfw/include -KPIC -O -I/opt/SUNWspro/prod/include -c dgeMatrix.c -o dgeMatrix.o cc -m64 -I/home/choid/bin/R281/lib/R/include -I./UFconfig -I/usr/local/include -I/opt/SUNWspro/prod/include -I/usr/include -I/usr/sfw/include -KPIC -O -I/opt/SUNWspro/prod/include -c dpoMatrix.c -o dpoMatrix.o cc -m64 -I/home/choid/bin/R281/lib/R/include -I./UFconfig -I/usr/local/include -I/opt/SUNWspro/prod/include -I/usr/include -I/usr/sfw/include -KPIC -O -I/opt/SUNWspro/prod/include -c dppMatrix.c -o dppMatrix.o cc -m64 -I/home/choid/bin/R281/lib/R/include -I./UFconfig -I/usr/loc...

...r [1:10] "(Intercept)" "FreqABCD.log.std" "LogitABCD.neg.log.std" "MIABCD.neg.log.std" ... ? .. ..$ : chr [1:5] "Estimate" "Std. Error" "df" "t value" ... ?$ sigma?????? : num 0.239 ?$ vcov??????? :Formal class 'dpoMatrix' [package "Matrix"] with 5 slots ? .. ..@ x?????? : num [1:100] 1.15e-03 3.40e-05 -5.12e-05 2.18e-05 3.65e-06 ... ? .. ..@ Dim???? : int [1:2] 10 10 ? .. ..@ Dimnames:List of 2 ? .. .. ..$ : chr [1:10] "(Intercept)" "FreqABCD.log.std" "LogitABCD.neg.log....

LS, I'm estimating multilevel regression models, using the lme-function from the nlme-package. Let's say that I estimated a model and stored it inside the object named 'model'. The summary of that model is shown below: Using summary(model)$tTable , I receive the following output: > summary(model)$tTable Value Std.Error DF t-value

I'm using FC4 and R 2.3.1 to fit a mixed effects logistic regression. The response is 0/1 and both the response and the age are the same for each pair of observations for each subject (some observations are not paired). For example: id response age 1 0 30 1 0 30 2 1 55 2 1 55 3 0 37 4 1 52 5 0 39 5 0 39 etc. I get the

r-devel list, I have recently written an R package that solves a linear least squares problem, and computes the multivariate normal density function. The bulk of the code is written in C, with interfacing code to the BLAS and Lapack libraries. The motivation here is speed. I ran into a problem computing the determinant of a symmetric matrix in packed storage. Apparently, there are no explicit

Hi helpeRs, I've been working through some examples in Pinhiero & Bates( 2000) trying to understand how to translate to the new Lme4 syntax but without much luck. Below is what I think I should do, but either the answers don't come out the same or I get errors. In the Oxide problems I'm particularly interested in obtaining the levels coeficients but this options no longer seems

The bulk is on my forums - the final post for today is: Results to date: A. It looks like I am going to need a newer compiler for C - xlc/xlC V11 apparently does not understand this code: "/data/prj/cran/R-3.2.3/src/main/memory.c", line 2149.31: 1506-046 (S) Syntax error. I will have to check if R-devel has different code before asking for assistence. +2139 #ifdef