similar to: nested anova diagnostics

Hello all. This may be a trivially simple question to answer, but I'm a little bit stumped with respect to the calculation of the F statistics in nested anovas in R. If I understand correctly, the F statistic for the among-subgroups but within groups hypothesis is calculated as MS_subgroups/MS_error, while the F statistic for the factor is calculated as MS_factor/MS_subgroups (I'm

Hi all, I have a questions regarding mixed effects models: I'm trying to plot e.g. fitted vs residuals for each level of the random effects, and i'm getting the same error. I guess this might be a problem of the graphic capabilities of R. Is there any way to obtain those plots? library(nlme) attach(ergoStool) names(ergoStool) [1] "effort" "Type"

Dear R users, I have more of a statistical/econometric question than straight R one. I have a data set with the discrete hazard rate of small firms survival on 400 counties over a period of 9 years. This data was generated using census information from the VAT registration number of each one of these business. I would like to analyze the effect of regional factors (deprivation index, real

Hi, Does anybody have a good tip of how to treat bimodal data to perform statistical analyses? My data set ranges from -1 to 1 (any values are posssible in between) and most data are either close to -1 or close to 1. They are the results of a two choice experiment where individuals could choose more than once in either direction and scores were calculated. Simone Simone Immler

I have been trying to make use of the odesolve library on my university's Linux grid - currently R version 2.0.1 is installed and the system runs 64-bit Scientific Linux based on Redhat. I cannot seem to get lsoda working when I define the model as a shared C library. For example, the following snippet uses the mymod.c example bundled with the package: ### START rm(list=ls())

Dear all Somebody may have asked this before but I could not find any answers in the web so let me ask a question on lme. When I have a fixed factor of, say, three levels (A, B, C), in which each level has different size (i.e. no. of observations; e.g. A>B>C). When I run an lme model, I get the same degree of freedom for all the contrast t-tests (e.g. AvsB or BvsC). I have tried this to

Dear Dr Skaug and R users, just discovered glmm.admb in R, and it seems a very useful tool. However, I ran into a problem when I compare two models: m1<-glmm.admb(survival~light*species*damage, random=~1, group="table", data=bm, family="binomial", link="logit") m1.1<-glmm.admb(survival~(light+species+damage)^2, random=~1, group="table", data=bm,

Dear all - I have followed the thread the reply to which was lead by Thomas Lumley about using pweibull to generate fitted survival curves for survreg models. http://tolstoy.newcastle.edu.au/R/help/04/11/7766.html Using the lung data set, data(lung) lung.wbs <- survreg( Surv(time, status)~ 1, data=lung, dist='weibull') curve(pweibull(x, scale=exp(coef(lung.wbs)),

Hi, I want to remove the entire row from the dataset if any of the row contains blank or NA my dataset look like State Year Y X2 X3 X4 X5 X6 State1 1960 27.8 397.5 42.2 50.7 78.3 65.8 State2 1960 29.9 413.3 38.1 52 79.2 66.9 State1 1961 29.8 439.2 40.3 54 79.2 67.8 State2 1961 30.8 459.7 39.5 79.2 69.6 State1 1962 31.2 492.9 37.3 54.7 77.4 68.7 State2 1962 528.6 38.1 80.2 73.6 State1 1963

Hi all - R2.0.1 on OSX;MASS library;nlme library I am trying to emulate the solution to a problem set that has normally been run in Genstat, using R. The problem that I am having at the moment is with the following glmm question (using glmmPQL from the MASS library): "We have two different forest habitats (first rotation thicket, and high forest) which we want to survey for the

Hi I'm getting the intercepts of the Random effects as 0. Please help me to understand why this is coming Zero This is my R code Data<- read.csv("C:/FE and RE.csv") Formula="Y~X2+X3+X4 + (1|State) + (0+X5|State)" fit=lmer(formula=Formula,data=Data) ranef(fit). My sample Data State Year Y X2 X3 X4 X5 X6 S2 1960 27.8 397.5 42.2 50.7 78.3 65.8 S1 1960 29.9 413.3 38.1

Hi I'm getting the intercepts of the Random effects as 0. Please help me to understand why this is coming Zero This is my R code Data<- read.csv("C:/FE and RE.csv") Formula="Y~X2+X3+X4 + (1|State) + (0+X5|State)" fit=lmer(formula=Formula,data=Data) ranef(fit). My sample Data State Year Y X2 X3 X4 X5 X6 S2 1960 27.8 397.5 42.2 50.7 78.3 65.8 S1 1960 29.9 413.3 38.1

Dear R users, I have been having problems getting believable estimates from anova on a model fit from lmer. I get the impression that F is being greatly underestimated, as can be seen by running the example I have given below. First an explanation of what I'm trying to do. I am trying to fit a glmm with binomial errors to some data. The experiment involves 10 shadehouses, divided between

Hi all - R2.0.1, OS X Perhaps while there is some discussion of lme going on..... I am trying to execute a glmm using glmmPQL from the MASS libray, using the example data set from McCullagh and Nelder's (1989, p442) table 14.4 (it happens to be the glmm example for GENSTAT as well). The data are binary, representing mating success (1,0) for crosses between males and females from two

Let's try this again. I have a Quintum AX Tenor gateway sending calls to Asterisk from BT analogue lines connected to FXO. The agents hear an echo on their side but incoming callers hear the conversation fine. I can't seem to find the problem. Anyone seen this issue before? <p style="margin: 0; padding: 0; border-collapse: collapse; font-family: Tahoma, Arial, Sans-Serif;

Dear all, I am trying to repeat an example from Sokal and Rohlfs "Biometry" -- Box 11.4, example of a randomized-complete-blocks experiment. The data is fairly simple: series genotype weight 1 pp 0.958 1 pb 0.985 1 bb 0.925 2 pp 0.971 2 pb 1.051 2 bb 0.952 3 pp 0.927 3 pb 0.891 3 bb 0.892 4

Hi i am having a problem with the 'integrate' function the function i want to integrate has the form sum(vector^x) i have defined the function with a for loop first - integrandtotest <- function(x) {a<-rep(0,len=2) for (i in 1:2) {a[i]<-t[i]^x} sum(a) } the results gives errors ########### Error in integrate(integrandtotest, lower = 0.1, upper =

Dear all, I have a large data set and would like to fit a logistic regression model using the bigglm function. I need to include an offset in the model but when I do this the bigglm function seems to ignore it. For example, running the two models below produces the same model and the offset is ignored bigglm(y~x,offset=z,data=Test,family=binomial(link = "logit"))

  Hi r-users,   I would like to use gev and my data (annual rainfall ) is as follows:   > head(dat,20) A B C D E F G H I J 1 45.1 41.5 58.5 50.1 46.0 49.1 37.7 49.1 59.8 54.0 2 50.3 39.8 49.4 56.4 49.4 48.8 42.1 49.8 49.4 58.3 3 41.7 39.3 44.6 39.1 35.7 41.5 40.8 40.8 38.5 45.6 4 50.7 33.9 48.4 28.2 35.5 39.1 61.4 17.0 30.7 38.3 5 39.3 30.6 46.9 23.8 25.8

Dear all, I am trying to calculate the proximity matrix for a data set with 16 variables and 6804 observations using random forests. I have a Pentium 4, 3.00GHz processor with 1 GB of RAM. When I use the command randomForest(data.scale,proximity=T) I get the warning message Error: cannot allocate vector of size 361675 kb Is this because I have reached the limit of what my computer is