similar to: Installing CAIC

Dear all, For a project I need to calculate the conditional AIC of a mixed effects model. Luckily, I found a reference in the R help forum for a function to be used: CAIC <- function(model) { sigma <- attr(VarCorr(model), 'sc') observed <- attr(model, 'y') predicted <- fitted(model) cond.loglik <- sum(dnorm(observed,

Hello, I estimated a model using SEM package in R, which was fit to a raw moment matrix, and includes an intercept term. The only goodness of fit statistics that are output are Model Chisquare, AIC, AICc, BIC, CAIC, and normalized residuals. How can I get the other goodness of fit statistics, like adjusted goodness of fit, RMSEA, and R-squared? And how can I get the final value of the

Hi all, This was posted originally on r-sig-mixed-models, but I thought I would post here as well as it might be of more general interest. With a colleague, I have been trying to implement the Conditional AIC described by Vaida and Blanchard 2005 Biometrika, "Conditional Akaike information for mixed-effects models". This quantity is derived in a way analogous to the AIC, but is

Dear R users, Would it be difficult to change the code of stepAIC (from the MASS library) to use AICc instead of AIC? It would be great to know of someone has tried this already. Best wishes Christoph.

Please help me in writing the R code for this problem. I've been solving this for 4 days. It was hard for me to solve it. It's a simulation problem in R. The problem is My true model is a normal mixture which is given as 0.5 N(-0.8,1) + 0.5 N(0.8,1). This model has two components. I will get a random sample of size 100 from this model. I will do this 300 times. That means, I will have

Dear R users: I'm generating the following survival data: set.seed(123) n=200 #sample size x=rbinom(n,size=1,prob=.5) #binomial treatment v=rgamma(n,shape=1,scale=1) #gamma frailty w=rweibull(n,shape=1,scale=1) #Weibull deviates b=-log(2) #treatment's slope t=exp( -x*b -log(v) + log(w) ) #failure times c=rep(1,n) #uncensored indicator id=seq(1:n) #individual frailty indicator

Hello, I have conducted an SEM in which the resultant standardized path coefficients are much higher than would be expected from the raw correlation matrix. To explore further, I stripped the model down to a simple bivariate relationship between two variables (NDVI, and species richness), where it's my understanding that the SEM's standardized path coefficient should equal the correlation

Dear Gustaf, I wish you a happy new year! I am a PhD Staff at LSHTM and I have been trying to use WHO anthro macros in R with the file you posted on internet but when I double check the WHZ, HAZ, WAZ output with output from ENA or STATA, it differs slightly. Do you have any idea why that might be the case? Many thanks Severine

Hello. I am trying to write an R function to plot the survival function (and associated hazard and density) for a Siler competing hazards model. This model is similar to the Gompertz-Makeham, with the addition of a juvenile component that includes two parameters---one that describes the initial infant mortality rate, and a negative exponential that describes typical mortality decline over the

Hi all, some years ago, I sent a question to the mailing list regarding the WHO anthro macros. Since I've now received three mails asking how I solved it, I thought I'd cc R-help in for future reference. Attaching a zip file with the relevant code parts that I used that I'm not sure gets through (if anyone has recommendations on how to manage such files for the list, I'd be

Hi all, I've got a project where I have to calculate weight-for-age Z-scores, preferably using the WHO standards. WHO have been very nice to publish macros for doing this in both STATA,SPSS, SAS and Splus formats (see http://www.who.int/childgrowth/software/en/), but for some reason have chosen not to use the free R alternative to Splus. In the Splus zipfile there are nine datafiles with a

I am using two mixed effect models. Firstly, what I am trying to do is to compare green roofs abundance with brownfield, green roof with green space abundance, and finally green space with brownfield abundance. I am unsure if I have done the correct model. I have to use a mixed effect model because my data is nested. This is the code and output >

Hi, i'm looking to test if two distance matrices are statistically different from each others. These two matrices have been computed on the same set of data (several population samples) 1. using a particular genetic distance 2. weighting that genetic distance with an extra factor (we can look at this as one set is computed before applying a treatment and the second one after applying

Hi, I'm new to R. Apologies if this is a simple query, I've searched the mailing lists and docs but can't find a solution to my problem. I'm trying to make some packages available on our intranet. During development the 'intranet' is a webserver running on localhost. * When I call "install.packages" I get a mesage about not being able to access 'index

Greetings. There's a coalescing group which is working to imitate Dirk's fine translation of CRAN to APT, in the Fedora/Red-Hat flavored land. One of the things we're trying to do is implement as much as is possible in R directly, and use as much of the existing package management infrastructure as we can. To this end, I humbly submit a few small patches to that infrastructure.

Greetings, and apologies for the braino leading to my first posting to the general list. I'm reposting here, and changing how I attached the files; they were hard to work with in the other post. There's a coalescing group which is working to imitate Dirk's fine translation of CRAN to APT, in the Fedora/Red-Hat flavored land. One of the things we're trying to do is implement as

Hi R help, What is the best way to plots surfaces in 3 dimensions? I also have the following availability problem with "plot3d" and "scatterplot3d", and "wireframe": install.packages("scatterplot3d") Warning: unable to access index for repository http://cran.ch.r-project.org/bin/macosx/leopard/contrib/2.10 Warning message: In getDependencies(pkgs,

Hello List, Could someone tell why I can't install the Himsc and survey packages for R version 2.13.0 (2011-04-13)? What am I doing wrong here? Thanks, Pradip > install.packages ("Hmisc", dependencies=TRUE) --- Please select a CRAN mirror for use in this session --- Warning: unable to access index for repository http://watson.nci.nih.gov/cran_mirror/bin/windows/contrib/2.13

? Hi: I get?a couple of warnings? when trying to download gridExtra: install.packages("gridExtra",repos=http://R-Forge.R-project.org) ?Warning: unable to access index for repository http://R-Forge.R-project.org/bin/windows/contrib/2.10 Warning message: In getDependencies(pkgs, dependencies, available, lib) : ? package ?gridextra? is not available I would like to download the binary

Dear Rxperts.. I am trying to install a few packages in R (version 2.14.2) on a *64-bit Window 7* workstation. Some of the errors are posted below.. Both 32-bit and 64-bit versions of R were installed. I did not get errors while installing packages through *32-bit R*, except that it could not update some of the packages (already installed with 64-bit R) *A*) When I start *64-bit R (2.14.2)*