similar to: how to obtain "par(ask=TRUE)" with trellis-plots

Dear list members, the following hlm was constructed: hlm <- groupedData(laut ~ design | grpzugeh, data = imp.not.I) the grouped data object is located at and can be downloaded: www.anicca-vijja.de/lg/hlm_example.Rdata The following works: library(nlme) summary( fitlme <- lme(hlm) ) with output: ... AIC BIC logLik 425.3768 465.6087 -197.6884 Random effects:

Dear altogether, is it possible to integrate "weights" arguments within lmer to incorporate statements to handle heteroscedasticity as it is possible with lme? I searched the R-archive but found nothing, insofer I assume it is not possible, but as lmer is under heavy develpoment, maybe something changed or is solved differently. Thus my question: While encountering heavy

Dear alltogether, I am searching for a way to determine "effect size" in multi-level models by using lme(). Coming from Psychology, for ordinary OLS there are measures (for meta-analysis, etc.) like CohensD <- (mean_EG - mean_CG) / SD_pooled or (p)eta^2 <- SS_effect / (SS_effect + SS_error) I do not intend to lead a discussion of the usefulness of such measures as long as

Dear alltogether, I experience very strange behavior of imputation of NA's with the NORM library. I use R 2.2.0, win32. The code is below and the same dataset was also tried with MICE and aregImpute() from HMISC _without_ any problem. The problem is as follows: (1) using the whole dataset results in very strange imputations - values far beyond the maximum of the respective column, >

Hi all, does anyone know how to automatically update starting values in R? I' m fitting multiple nonlinear models and would like to know how I can update starting values without having to type them in. thank all --- On Fri, 12/26/08, r-help-request@r-project.org <r-help-request@r-project.org> wrote: From: r-help-request@r-project.org <r-help-request@r-project.org> Subject:

Dear alltogether, two lme's, the data are available at: http://www.anicca-vijja.de/lg/hlm3_nachw.Rdata explanations of the data: nachw = post hox knowledge tests over 6 measure time points (= equally spaced) zeitn = time points (n = 6) subgr = small learning groups (n = 28) gru = 4 different groups = treatment factor levels: time (=zeitn) (n=6) within subject (n=4) within smallgroups

Hi all: I am a user of "JM" package. Here's the problem of "sample size". The warning is: Error in jointModel(fitLME, fitSURV_death, timeVar = "time", method = "piecewise-PH-GH") : sample sizes in the longitudinal and event processes differ. According to the suggestion of "missing data",I use the same data set(data_JM) without any

Dear altogether, I tried local regression with the following data. These data are a part of a bigger dataset for which loess is no problem. However, the plot shows extreme values and by looking into the fits, it reveals very extreme values (up to 20000 !) although the original data are > summary(cbind(x,y)) x y Min. :1.800 Min. :2.000 1st Qu.:2.550

I contacted the package developer and that lead to me removing events at time 0 (or subjects with only 1 longitudinal measurement). I then still had the error message "Can't fit a Cox model with 0 failures" which I have managed to avoid by adding 1.8*10^(-15) to all my survival times, any number greater than this also works but nothing smaller! Any explanation of this would help!

Hello All, While trying to install Matrix package following error message came out: /usr/bin/ld: cannot find -lblas-3 collect2: ld returned 1 exit status make: *** [Matrix.so] Error 1 ERROR: compilation failed for package 'Matrix' ** Removing '/usr/lib/R/library/Matrix' ** Restoring previous '/usr/lib/R/library/Matrix' It seems that someting is missing. Have anybody

I have used LME to fit a mixed effects model on my data. The data has 274 subjects with 1 to 6 observations per subject. Time is not linearly associated with the outcome, so I used ns to fit a natural cubic spline with 3 auto knots. Subject and the natural cubic time of spline are both treated as random effects. This model has run without any problem, but now I would like to plot trajectories for

Hello alltogether, how is it possible to assign the colnames of a data.frame to a function called by apply, e.g. for labeling a plot? Example: I want to plot several qqnorm-plots side by side and there should be a maintitle for each qqnorm-plot which is identical to the respective colname. I checked, but the column which is processed by the function called by apply does not contain a colname

Yes, keep the $(lsb-release -cs) there. This allows you to change the version of the OS without having to change anything else (for instance we are migrating all of our docker to 24.04 LTS right now, so we only have to modify the FROM entry and all of our Dockerfile will still work). From: Gerber, Lauren J <lauren.gerber at helsinki.fi> Date: Thursday, July 11, 2024 at 2:21?AM To:

Dear R-list, I try to transform a mathematica script to R. #######relevant part of the Mathematica script (* p_sv *) dd = NN (DsD - DD^2); lownum = NN (L-DD)^2; upnum = NN (H-DD)^2; low = lownum/(2s^2); up = upnum/(2s^2); psv = NIntegrate[1/(s^NN) Exp[-dd/(2s^2)] (Gamma[1/2,0,up] + Gamma[1/2,0,low]),{s,sL,sH}, MinRecursion->3]; PSV = psv/Sqrt[2NN]; Print["------------- Results

Hi, I am trying to extract the coefficients matrix from a gls summary. Contrary to the lm function, the command fit$coefficients returns only the estimates of the model, not the whole matrix including the std errors, the t and the p values. example: ctl <- c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14) trt <- c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69) group <-

Hi Marco, Thank you for the follow up email. In the line below, do I replace $(lsb-release -cs) with jammy (Ubuntu 22.04) or leave it alone: Add-apt-repository ?deb https://cloud.r-project.org/bin/linux/ubuntu jammy-cran40/? It?s good to know about installations for Bioconductor. I am not currently using it, but I may need it in the future for other projects. Br, Lauren -- Lauren J. Gerber

Lauren, On 15 October 2024 at 13:51, Gerber, Lauren J wrote: | I just tried fixing my script and added library(BiocManager) before library | (ggbreak), then rebuilt the image, and am still getting the same error. You need _two_ distinct commands. First install.packages("BiocManager") and second once that is done use it to install 'ggbreak' from BioConductor

Lauren, little oversight from my Dockerfile, if you want to follow the CRAN steps and get the R from the cloud.r-project.org server, you want to replace && add-apt-repository -y "ppa:marutter/rrutter4.0" With && add-apt-repository "deb https://cloud.r-project.org/bin/linux/ubuntu $(lsb_release -cs)-cran40/" The ppa repo is Michael Rutter personal repo

Don't you have to load BiocManager before using it? -----Original Message----- From: R-SIG-Debian <r-sig-debian-bounces at r-project.org> On Behalf Of Gerber, Lauren J Sent: Tuesday, October 15, 2024 9:15 AM To: r-sig-debian at r-project.org Cc: dirk at eddelbuettel.com Subject: Re: [R-sig-Debian] R Package: 'ggbreak' in Dockerfile Hi All, Does anyone have experience

Hi All, Does anyone have experience successfully installing and loading the R package ?ggbreak? in a Docker container? I ask because I have successfully used the ?ggbreak? package in a script outside the Docker container (on my local computer), but for some reason, I get the following error message when trying to load it using library(ggbreak) or library(?ggbreak?) inside the same script that