search for: geraci

...tes that Rgui.exe is using 25% of the CPU speed. There is no other process that is slowing down R. Does anyone know if it's a limitation of R? If so, can I set up the program for the CPU usage like I do for the memory usage? Is there a version of R for servers? Thank you in advance, Marco Geraci --------------------------------- [[alternative HTML version deleted]]

...1:M){ lme(..) } and for each "i" i get the warning message >Fewer observations than random effects in all level 1 groups >in: ... I know I'm using "fewer observations...", I just don't want to see the message printed on the screen during the loop. Thanks, Marco Geraci > sessionInfo() R version 2.2.1, 2005-12-20, i386-pc-mingw32 attached base packages: [1] "methods" "stats" "graphics" "grDevices" "utils" "datasets" [7] "base" other attached packages: nlme "3.1...

...likelihood Sum_i {w_i * logLik_i}, since 'weight' describes the with-in heteroscedasticity structure. I think I need something like the option 'iweight' (importance weight) for the command 'xtreg' of Stata. Any suggestion for R? Thanks in advance, Marco Geraci > sessionInfo() R version 2.2.1, 2005-12-20, i386-pc-mingw32 attached base packages: [1] "methods" "stats" "graphics" "grDevices" "utils" [6] "datasets" "base" other attached packages: nlme &q...

Good evening, I am Marta Colombo, student at Milan's Politecnico. Thank you very much for your kindness, this mailing list is really useful. I am using the function kde2d for two-dimensional kernel density estimation and I'd like to know something more about this kind of density estimator. In particular I'd like to know: what bandwidth is used ? Thank you in advance for your attention

..., but I guess that this problem is related to different evaluation frames, the one invoked by 'myfunc' and the one invoked by 'update' (both for 'lme' and 'groupedData' objects). Any suggestion to overcome this problem would be greatly appreciated Marco Geraci R version: 2.2.1 nlme version: 3.1-68.1 OS: Window XP --------------------------------- Photo Books. You design it and we’ll bind it! [[alternative HTML version deleted]]

...near quantile regression models with random effects. Commands for independent and count data are also available. A few bugs have been fixed since lqmm 1.0. Updates of this package are (usually) made available on https://r-forge.r-project.org/R/?group_id=1396 prior to submission to CRAN. -- Marco Geraci PhD CStat Lecturer in Biostatistics Institute of Child Health, University College London London WC1N 1EH, UK Web: http://marcogeraci.wordpress.com [[alternative HTML version deleted]] _______________________________________________ R-packages mailing list R-packages at r-project.org https://...

...near quantile regression models with random effects. Commands for independent and count data are also available. A few bugs have been fixed since lqmm 1.0. Updates of this package are (usually) made available on https://r-forge.r-project.org/R/?group_id=1396 prior to submission to CRAN. -- Marco Geraci PhD CStat Lecturer in Biostatistics Institute of Child Health, University College London London WC1N 1EH, UK Web: http://marcogeraci.wordpress.com [[alternative HTML version deleted]] _______________________________________________ R-packages mailing list R-packages at r-project.org https://...

Hi, Dear R users I have problem with the following code. The matrix result must be a matrix (3x3). But I have obtained a matrix(3x1) and I don't know why. So, I need your help Best regards ##################################################################### taille <- function (delta, t, prob = 0.2) { niv.conf <- c(0.90, 0.95, 0.99) if(niv.conf <- 0.90) { t <- 1.645 }

Dear R Users I would like to use optim function to optimize a function. I read help but I couldn't find what I need: is it possible to get information after each iteration, for example as there is in MATLAB: Gradient's Iteration Func-count f(x) Step-size infinity-norm 0 24 388.976

Hi All, I have a list of matrices: > x [,1] [,2] [1,] 1 4 [2,] 2 5 [3,] 3 6 > y [,1] [,2] [,3] [,4] [,5] [,6] [1,] 18 21 24 27 30 33 [2,] 19 22 25 28 31 34 [3,] 20 23 26 29 32 35 > z =list(x,y) I want to create a second list that is has a subset each matrix in the list subsetting so I get the 2nd and 3rd row of each (and

Dear R users, I am estimating Poisson mixed models using glmmPQL (MASS) and lmer (lme4). We know that glmmPQL do not provide the correct loglikelihood for such models (it gives the loglike of a 'pseudo' or working linear mixed model). I would like to know how the loglike is calculated by lmer. A minor question is: why do glmmPQL and lmer give different degrees-of-freedom for the same