similar to: likelihood ratio tests using glmmPQL

I am analysing binomial data using a generalised mixed effects model. I understand that if I use glmmPQL it is not appropriate to compare AIC values to obtain a minimum adequate model. I am assuming that this means it is also inappropriate to use AIC values from lmer since, when analysing binomial data, lmer also uses PQL methods. However, I wasn't sure so please could somebody clarify

Can you please tell me the minimum system spec for R on a windows machine? Thanks _______________________________________ Matthew Longley IT Administrator Zoological Society of London Regent's Park London NW1 4RY Charity No: 208728 Tel: 020 74496412 Fax: 020 74496294 email: Matthew.Longley at zsl.org WWW: www.zsl.org _______________________________________

Hello, I was wondering if anyone knew of an R function/R code to plot bivariate (3 dimensional) partial dependence plots in random forests (randomForest package). It is apparently possible using the rgl package (http://esapubs.org/archive/ecol/E088/173/appendix-C.htm) or there may be a more direct function such as the pairplot() in MART (multiple additive regression trees)? Many

Hi, I have been trying to plot density plots using the example on: http://addictedtor.free.fr/graphiques/graphcode.php?graph=139 I used to use this function, but I cannot get any old code or even the example to work. library("geneplotter") require("RColorBrewer") x1 <- matrix(rnorm(1e4), ncol=2) x2 <- matrix(rnorm(1e4, mean=3, sd=1.5), ncol=2) x <-

Douglas Bates wrote: >From 0.62 onward you should not have to create a symbolic link in >/usr/local/bin. It should be that you can run > cd $RSOURCE > ./configure --prefix=/usr/local > make install >and you will end up with the R script installed in /usr/local/bin and >all the files needed to run R in /usr/local/lib/R. > >Can you tell what the prefix is set to after

I was wondering if someone could help me with a problem getting locfit running. I am running Digital Unix 3.2D on an alphastation, with R 0.61.1 I downloaded the CRAN distribution of locfit. I had no problems installing it with R INSTALL (I installed it into my library rather than the default base library; the output is attached at the end of this email message). I then tried to test it with the

Hi there I'm trying to construct a model of mortality risk in 2D space that requires numerical integration of a hazard function, for which I'm using the integrate function. I'm occasionally encountering parameter combinations that cause integrate to terminate with error "Error in integrate... the integral is probably divergent", which I'm not sure how to interpret. The

The problem is the following # make docs Make: Don't know how to make ../src/library/*/man/*.Rd. Stop. *** Exit 1 (ignored) # I am using a DEC alpha 200/4/233 with Digital Unix 3.2D My previous version of R was R-0.50-a1. I have installed around 3 or 4 previous versions on the same system and never had a problem (but, then the Makefile was very different). The first make to compile

The problem is the following # make docs Make: Don't know how to make ../src/library/*/man/*.Rd. Stop. *** Exit 1 (ignored) # I am using a DEC alpha 200/4/233 with Digital Unix 3.2D My previous version of R was R-0.50-a1. I have installed around 3 or 4 previous versions on the same system and never had a problem (but, then the Makefile was very different). The first make to compile

I am trying to set up R-0.62.1 on an alphastation running digital unix 3.2D. I had R-0.61.1 working fine, but have had a number of problems getting 0.62.1 going.... The next problem is that plotting to any device causes R to hang. Looking at processor usage it would appear that it gets stuck in an infinite loop rather than pausing for input. I have managed to track the problem down to plot.new().

Hi All I have a column/variable called time difference. It has a whole list of numbers from 0 through to the hundreds eg 236. I want to assign a corresponding "name" to each variable from a predefined list: Month or less, 1 -2 months, 2-3 months etc So the result would look something like: Time Difference Month 1 Month or

Dear R Community, I'm trying to exploit the elegance of R by doing the following pseudocode routine without a WHILE or FOR loop in R: for i = 1 to length-1 for j = (i+1) to length print a[i], a[j] That is, I want i and j to be the indices of a half-matrix 1 2, 1 3, 1 4, ..., 1 length, 2 3, 2 4, ..., 2 length, 3 4, ..., 3 length 1. Can this be done with the 'whole

Hi all, I am trying to predict a target variable that takes values 0 or 1 using the rpart command. In my initial dataset I have few positive observations of the target variable; therefore I have oversampled the rare event by a multiple of 6 (i.e. from 762 to 4572). However, in my results, I end up with a number of positives in one of the terminal nodes that is not divisible by 6. As I have the

Hi- I am trying to fit a log function to my data, with the ultimate goal of finding the second derivative of the function. However, I am stalled on the first step of fitting a curve. When I use the following code: FG2.model<-(nls((CO2~log(a*Time)+b), start=setNames(coef(lm(CO2 ~ log(Time), data=FG2)), c("a", "b")),data=FG2)) I get the following error: Error in

I am trying to make the R-0.62.1 distribution on a DEC alpha running Digital Unix 3.2d. I have run into a few problems..... (1) the Makefile in src/graphics need to have INCLUDES= -I../include not -I ../include (2) I have the standard DEC fortran 90 compiler which needs to link R.binary with a -nofor_main flag. The appropriate slot in config.site appears to be ignored, which means Makeconf

Hey, Can anyone answer this question. I am working with really large datasets and most of the programs I have been running take quite some time. I heard that R may be faster in Unix. I sthis true and if so can anyone reccomend which system and requirements may allow things to go faster for? Thanks!! Elizabeth Lawson --------------------------------- [[alternative

Hi All, could anybody help me to understand what is this error means ? mydata=read.table("C:/Documents and Settings/angieb/Desktop/CommercialGL/cl_ilf_claimdata.csv",header=TRUE,sep=",") > names(mydata) [1] "ILFTable" "liabLimit" "AnnAggLimit" "DedAmt" "Loss" "TIL" >

Hey, I am estiamting a partially linear model y=X\beta+f(\theta) where the f(\theta) is estiamted using wavelets. Has anyone heard of methods to test if the betas are significant or to address model fit? Thanks for any thoughts or comments. Elizabeth Lawson __________________________________________________ [[alternative HTML version deleted]]

Hello all, I would like to perform SVM-RFE (Guyon et al. 2002) in R and have only found one implementation of this algorithm. The function belongs to the MCRestimate package but when I try to use it I encounter a problem - the function appears to be missing a required package or other function that I simply cannot find available anywhere. Here is my session info followed by a simple example

Ms.Karunambigai M PhD Scholar Dept. of Biostatistics NIMHANS Bangalore India ----- Forwarded Message ---- From: karuna m <m_karuna2002@yahoo.com> To: R-help@r-project.org Sent: Thu, 26 May, 2011 10:55:02 AM Subject: likelihood ratio test Dear R-help, Can anybody tell me which R package has Lo-Mendell Rubin LR test and Bootstrap LR test to compare the model fit between k class and k+1