similar to: conflict of package "Matrix" and summary of lme object

Using ranef() (package nlme, version 3.1-75) with an 'lme' object I can obtain random effects for intercept and slope of a certain level (say: 1) - this corresponds to (say level 1) "residuals" in MLWin. Maybe I'm mistaken here, but the results are identical. However, if I try to get the standardized random effects adding the paramter "standard=T" to the

I wonder whether the functions for skewness and kurtosis in the e1071 package are based on correct formulas. The functions in the package e1071 are: # -------------------------------------------- skewness <- function (x, na.rm = FALSE) { if (na.rm) x <- x[!is.na(x)] sum((x - mean(x))^3)/(length(x) * sd(x)^3) } # -------------------------------------------- and #

Trying to select a subset of cases (rows of data) I encountered several problems: Firstly, because I did not read the help to read.spss() thoroughly enough, I treated the data read as a data frame. For example, dr2000 <- read.spss('myfile.sav') d <- subset(dr2000,RBINZ99 > 0) and thus received an error message (Object "RBINZ99" not found), because dr2000 is not a

I receive R-help messages in digest form that makes it difficult to answer a post. I noticed that my answer is not added to the thread (instead, a new thread is started) although I use the same subject line (starting with "Re: ") as the original post. Is there a solution (I prefer the digest to separate messages for several reasons and don't want to change my email reader)? The

Although I tried to find an answer in the manuals and archives, I cannot solve this (please excuse that my English and/or R programming skills are not good enough to state my problem more clearly): I want to write a function with an indeterminate (not pre-defined) number of arguments and think that I should use the "..." construct and the match.call() function. The goal is to write

Just out of curiosity (and I hope that my question does not confuse those who intuitively understood the par(new=TRUE) statement right in the beginning): I would like to know whether beginners (learning R) had the same difficulty as I had: Intuitively I thought that par(new=TRUE) would draw a new plot (and not into an already existing plot) and that par(new=FALSE) would not draw a new plot

The update of VR by downloading VR_7.1-6.zip and using install.packages (from local zip files) fails with the following error message: Error in file(file, "r") : unable to open connection In addition: Warning message: cannot open file `VR/DESCRIPTION' Other packages can be installed without problems, except of dse_2003.4-1 with a similar error message. Why? Operating System:

Is it possible to avoid the loop in the following function (or make the function otherwise more efficient) and can someone point me to a possible solution? (It would be great if hours could be reduced to seconds :-). # --------------------------------------------- RanEigen=function(items=x,cases=y,sample=z) { X=matrix(rnorm(cases*items),nrow=cases,byrow=F) S=crossprod(X-rep(1,cases) %*%

I'm sorry, but my previous message, as often happens, some brackets were wrong: Here are the correct formulas: sd = 1/n * sum((x-mean(x))^2) # (1) sd = 1/(n-1) * sum((x-mean(x))^2) # (2) This also occured in the last paragraph. Dirk ************************************************* Dr. Dirk Enzmann Institute of Criminal Sciences Dept. of Criminology Edmund-Siemers-Allee 1 D-20146

Can someone show me how to specify zero inflated poisson and zero inflated negative poisson models in R? I would like to replicate the example given in Long (1997: Regression Models for Categorical and Limited Dependent Variables) in Chapter 8.5 (pp. 242-247). TIA Dirk ************************************************* Dr. Dirk Enzmann Criminological Research Institute of Lower Saxony

Begin forwarded message: > Subject: R: Studdy Missing Data, differentiate between a percent > with in the valid answers and with in the different missing answers > take a look at the Website from Dirk Enzmann: > http://www2.jura.uni-hamburg.de/instkrim/kriminologie/Mitarbeiter/ > Enzmann/Software/Enzmann_Software.html > He provides a function called "freq.r" - maybe

Fine with me as long as eliminating the inconveniences associated with it can be put on the roadmap. The alias instability and the fact that the user has no way to know if s/he should do ?`foo,numeric-method` or ?`foo,numeric,ANY-method` to find the method has been a long-standing problem. H. On 3/21/19 21:29, Michael Lawrence wrote: If we started over, I'd try to avoid this sort of

Hi All. I need 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

Hi, I have sparse (tridiagonal) matrices, and I use the Matrix package for handling them. For certain computations, I need to either set the first row to zero, or double the last row. I find that neither operation preserves sparsity, eg > m <- Diagonal(5) > m 5 x 5 diagonal matrix of class "ddiMatrix" [,1] [,2] [,3] [,4] [,5] [1,] 1 . . . . [2,] . 1

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

Hi, I got the following error. Sorry but this time I couldn't reproduce it with a simple chunk of code: .TM.repl.i.2col(): drop 'matrix' case ... Error in .nextMethod(x = x, i = i, j = j) : 'i' has no integer column number should never happen; please report In addition: Warning messages: 1: Ambiguous method selection for "%*%", target

I hesitate to post this question in the light of recent threads, indeed I have hesitated for several weeks, however I have come to a full stop and really need some help if I am going to progress. I am a new user of R for medical statistics. I have attempted to read all the relevant documents, but would welcome any suggestions as to what I have missed. I am trying to contruct "table 1"

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

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%