similar to: problem with strata in boot

Hello, I am new to R, and am having trouble getting the output I want from a stratified bootstrap. I didn't recieve a reply the first time I posted this question so I have tried to make it more clear. My data frame (denboot) is set up as follows: SITE cswa parea treat 1 BeanA 3 1.20 m 2 BeanBC 3 1.05 m 3 BeanD 1 0.93 m 4 BlackB 1 1.01

Hello all, Is AIC calculated incorrectly in lmer? It appears as though it uses AIC = -2*logLik - 2*#parms, instead of -2*LogLik + 2*#parms? Below is output from one of many models I have tried: Generalized linear mixed model fit using PQL Formula: cswa ~ pcov.ess1k + (1 | year) Data: ptct50.5 Family: poisson(log link) AIC BIC logLik deviance 224.8466 219.19 -114.4233 228.8466

Hello, I'm trying to plot a fitted lm() line on a plot when the one explanatory variable is log transformed and log="x". I get different lines using abline and predict.lm(). #Example x <- 1:100 y <- rnorm(100) plot(y ~ x, log="x") abline(lm(y ~ log(x))) lines(x, predict(lm(y ~ log(x))), lwd=2) I'm sure I'm missing something but could someone tell me which

Hello, I would be very grateful if someone could give me a hand with my split plot design problems. So here is my design : I am studying the crossed-effects of water (wet/dry) and mowing (mowed/not-mowed = nm) on plant height (PH) within 2 types of plant communities (Xerobromion and Mesobromion) : - Within each type of communities, I have localised 4 blocks - In each block, I have defined

Hello All. I must be missing the really obvious here: mm <- function(d, i) median(d[i]) b1 <- boot(gravity$g, mm, R = 1000) b1 b2 <- boot(gravity$g, mm, R = 1000, strata = gravity$series) b2 Both b1 and b2 seem to have done (almost) the same thing, but it looks like the strata argument in b2 has been ignored. However, str(b1) vs str(b2) does show that the strata have been noted

Just browsing the documentation, and searching the list came up short... I have some unbalance data and was wondering if, in a "0" v "1" classification forest, if these options might yield better predictions when the proportion of one class is low (less than 10% in a sample of 2,000 observations). Not sure how to specify these terms... from the docs, we have: classwt: Priors

Hello, I was hoping that someone could answer a few questions for me (the background is given below): 1) Can the coxph accept an interaction between a covariate and a frailty term 2) If so, is it possible to a) test the model in which the covariate and the frailty appear as main terms using the penalized likelihood (for gaussian/t frailties) b)augment model 1) by stratifying on the variable that

Hi all, I am aware that crr lacks the "friendly" command structure of functions such as cph. All is clear to me about including covariates until I want to include a stratification term in the competing risk framework (no nice strat command). I am still a bit of a novice in R - I am looking for an example to help me with this, but can't seem to find one. Any advice appreciated (no

[code]>require("survival") > coxph(Surv(futime,fustat)~age + strata(rx),ovarian) Call: coxph(formula = Surv(futime, fustat) ~ age + strata(rx), data = ovarian) coef exp(coef) se(coef) z p age 0.137 1.15 0.0474 2.9 0.0038 Likelihood ratio test=12.7 on 1 df, p=0.000368 n= 26, number of events= 12 > coxph(Surv(futime,fustat)~age, ovarian, subset=rx==1)

Dear All, I got a warning message "X matrix deemed to be singular" in Cox model with a time dependent coefficient. In my analysis, the variable "SEX" is a categorical variable which violate the PH assumption in Cox. I first used the survSplit() function to break the data set into different time intervals, and then fit the model. The procedures can be described as follows:

> -----Original Message----- > From: Brian Roy [mailto:mister.roy@gmail.com] > I do this with a Toshiba Strata DK280. The model before the CTX. You > shouldn't have any problem doing this with the CTX either. Assuming > you don't have SMDI integration with your voicemail (you didn't state > that you did), you just need the phones to send trailing digits, and > you

Hi everyone, Been reading up on Asterisk, and very interested in learning more. I've googled and read the archives and haven't found anything definitive on support for this phone system. We have a fairly large investment in the system itself and the phones, but would love to get away from the voicemail system it forces on us. Can anyone provide any feedback on using this system with

I've been using strata(sampling) and found that if the dataframe to be sampled ("data") consists of repeated measures that are not sorted in order the run will fail on a given stratum. Ordering the dataframe prior to using strata() eliminates this problem. Thanks, Eric -- Eric Vander Wal Ph.D. Candidate University of Saskatchewan, Department of Biology, 112 Science Place,

How can I get the log odds associated with the levels in strata() within a clogit() model? I'm running R-1.9.0 on a Linux platform. I am using clogit() to run a Rasch model in Item Response Theory in psychometrics. Symbolically, the model is: logit(p_{j,k}) = \log({\Pr(p_{j,k}) \over \Pr(1-p_{j,k})}) = \theta_j - \alpha_k, That is, the log odds of answering an test item correctly is

Dear list members, I am working with a Cox ph model for the duration of unemployment. The event of interest in my analysis is getting employed. I have various background variables explaining this event: age, sex, education etc. I have multiple unemployment spells per person. I use a model with person-specific frailty terms in order to take into account the correlation of spells by the same

Sorry for the repost, but I've really been looking, and can't find any syntax direction on this issue... Just browsing the documentation, and searching the list came up short... I have some unbalanced data and was wondering if, in a "0" v "1" classification forest, some combo of these options might yield better predictions when the proportion of one class is low (less

Dear experts, when I call the step function for a coxph model with n covariates and a dicotomous variable included as strata, the first term removed by step is always the strata variable. This is not what I want and then I do a manual step updating the model minus the least significant covariate and testing with anova, until I have minimized the model. Is there a package were this can be done? or

hie l would like to create about ten strata l have tried using sample(c(1,2),1,replace=TRUE) but all l get is a vector of 1 and 2 if i change it to sample(c(1,10),1,replace=TRUE) it gives me a vector of 1 and 10 how do l go about it --------------------------------- [[alternative HTML version deleted]]

Hello, I'm trying to sample in Random Forests by a factor, but it is a regression problem and I can't figure out how to do this (I can only see how to sample by strata in classification). Thanks Jesse Lasky

My data.frame table consist of 3 variables (x,y and z) where each variable has 1000 units. I need to create 5 equal size strata according to one of the variable (let's say x) whereas units of x variable with a higher value have higher probability to be selected in a strata with a higher number (max strata number is 5). I've been trying different things so far and since I am fairly new to