search for: cohort

Consider this data structure (df1) ... Group Year PctProf FullYr 1 Never RF 2004 87 88 2 Cohort 1 2004 83 84 3 Cohort 2 2004 84 86 4 Cohort 3 2004 87 87 5 Cohort 4 2004 73 74 6 Never RF 2005 85 86 7 Cohort 1 2005 81 82 8 Cohort 2 2005 81 81 9 Cohort 3 2005 78 79 10 Cohort 4 2005 72 74 11 Never RF 2006...

...ne > To: ggplot2 at googlegroups.com > Date: Friday, April 17, 2009, 3:11 PM > I have a data frame consisting of a group ID, year, and > percent of > students proficient on a state test. > > df1 <- > > Group Year PctProf > 1 Never RF 2004 87 > 2 Cohort 1 2004 83 > 3 Cohort 2 2004 84 > 4 Cohort 3 2004 87 > 5 Cohort 4 2004 73 > 6 Never RF 2005 85 > 7 Cohort 1 2005 81 > 8 Cohort 2 2005 81 > 9 Cohort 3 2005 78 > 10 Cohort 4 2005 72 > 11 Never RF 2006 83 > 12 Cohort...

Consider an analysis of covariance involving age and cohort. The goal is to assess whether the influence of cohort depends upon the age. The simplest case involves data as follows value Age Cohort x1 ????? 1?????? 3 x2?????? 1?????? 4 x3?????? 1?????? 5 x4 ????? 2 ????? 3 x5 ????? 2 ????? 4 x6 ????? 2 ????? 5 etc. Age is a factor. The numeric res...

...a set looks like: > str(data) 'data.frame': 1004 obs. of 5 variables: $ site : int 0 0 0 0 0 0 0 0 0 0 ... $ sex : Factor w/ 2 levels "0","1": NA NA NA NA 1 NA NA NA NA NA ... $ age : Factor w/ 2 levels "0","1": 1 1 1 1 1 1 1 1 1 1 ... $ cohort: Factor w/ 11 levels "1996","2000",..: 11 11 11 11 11 11 11 11 11 11 ... $ birth : Factor w/ 3 levels "5","6","7": 3 3 2 2 2 2 2 2 2 2 ... I know that, particularly for one level of variable "cohort" (2004 value), it should be a strong e...

Hi, I have a data set with 999 observations, for each of them I have data on four variables: site, colony, gender (quite a few NA values), and cohort. This is how the data set looks like: > str(dispersal) 'data.frame': 999 obs. of 4 variables: $ site : Factor w/ 2 levels "1","2": 1 1 1 1 1 1 1 1 2 2 ... $ gender: Factor w/ 2 levels "0","1": NA NA 2 1 2 NA 1 2 2 NA ... $ colony: Factor w/ 2...

I tried to compare if cch() and coxph() can generate same result for same case cohort data Use the standard data in cch(): nwtco Since in cch contains the cohort size=4028, while ccoh.data size =1154 after selection, but coxph does not contain info of cohort size=4028. The rough estimate between coxph() and cch() is same, but the lower and upper CI and P-value are a little differ...

Just trying to understand how geom_abline works with facets in ggplot. By way of example, I have a dataset of student test scores. These are in a data table dt with 4 columns: student: unique student ID cohort: grouping factor for students (A, B, . H) subject: subject of the test (English, Math, Science) score: the test score for that student in that subject The goal is to compare cohorts. ## Code to generate dt library(data.table) ## cohorts: list of cohorts with number of students in each...

...(datadist='dd') bp.cat2 <- all.data2.stand$bp.cat2 u <- cr.setup(bp.cat2) u b.mean <-rep(all.data2.stand$b.mean, u$reps) r.mean <-rep(all.data2.stand$r.mean, u$reps) mean.ova.energy <- rep(all.data2.stand$mean.ova.energy, u$reps) y <- (u$y) # constructed binary response cohort <- u$cohort attach(all.data2.stand[u$subs,]) dd <- datadist(dd, cohort) ord.cr <- lrm(y ~ cohort + mean.ova.energy + b.mean + r.mean, x=TRUE, y=TRUE, na.action=na.delete) summary(ord.cr) p.cr <- predict(ord.cr, all.data2.stand, type='mean', codes=TRUE) pred.mean2 <- data....

Now that we have case cohort model , we have 1000 people and 50 cases Let the first 10 cases occur at the same time second 10 " third 10 " fourth 10 " fifth 10 " How easy is it to randomly sample 50 different c...

Hi all, I have done a backward stepwise selection on a full binomial GLM where the response variable is gender. At the end of the selection I have found one model with only one explanatory variable (cohort, factor variable with 10 levels). I want to test the significance of the variable "cohort" that, I believe, is the same as the significance of this selected model: > anova(mod4,update(mod4,~.-cohort),test="Chisq") Analysis of Deviance Table Model 1: site ~ cohort Model 2:...

Hello, I have written code to plot an xyplot as follows: library(lattice) xyplot(len~ageJan1|as.factor(cohort),groups=sex,as.table=T,strip=strip.c ustom(bg='white',fg='white'),data=dat, xlab="Age (January 1st)",ylab="Length (cm)",main="Linear models for male and female cod, by cohort",type='p', lwd=1.5,auto.key=list(text=c("Male","Female&...

...-------- Original from the R Book: http://books.google.com/books?id=8D4HVx0apZQC&lpg=PA799&ots=rQgd_8ofeS&dq=r%20coxph%20crawley&pg=PA799#v=onepage&q&f=false --------------------------------- My result: > summary(model1) Call: coxph(formula = Surv(death, status) ~ strata(cohort) * gapsize, data = seedlings) n= 60, number of events= 60 coef exp(coef) se(coef) z Pr(>|z|) gapsize -0.001893 0.998109 0.593372 -0.003 0.997 gapsize:strata(cohort)cohort=September 0.717407 2.049112...

--------- begin included message --------- I tried to compare if cch() and coxph() can generate same result for same case cohort data Use the standard data in cch(): nwtco Since in cch contains the cohort size=4028, while ccoh.data size =1154 after selection, but coxph does not contain info of cohort size=4028. The rough estimate between coxph() and cch() is same, but the lower and upper CI and P-value are a little differ...

> Now that we have case cohort model , we have 1000 people and 50 cases > Let the first 10 cases occur at the same time > second 10 " > third 10 " > fourth 10 " > fifth 10 " > How easy is it to randomly...

I have a question about stepAIC and extractAIC and why they can produce different answers. Here's a stepAIC result (slightly edited - I removed the warning about noninteger #successes): stepAIC(glm(formula = (Morbid_70_79/Present_70_79) ~ 1 + Cohort + Cohort2, family = binomial, data = ghs_70_79, subset = ghs_70_full),direction = c("backward")) Start: AIC=3151.41 (Morbid_70_79/Present_70_79) ~ 1 + Cohort + Cohort2 Df Deviance AIC <none> 1797.6 3151.4 - Cohort 1 1826.2 3178.0 - Cohort2 1 1826.3...

----- begin included message In case cohort study, we can fit proportional hazard regression model to case-cohort data. In R, the function is cch() in Survival package Now I am working on case cohort analysis with time dependent covariates using cch() of "Survival" R package. I wonder if cch() provide this utility or not? The cch()...

Hi All, I am in the middle of doing an analysis of a Case-Cohort design. I had three questions about the analysis: a) Does any one know of some public code for developing the patient risk sets (indexed by failure time) or is there a better way to organize the data? b) I was planning to use the Barlow weighting method. Has this or any other weighting method (...

...as 4 levels. I see r-help inquiries regarding odds ratios for what looked like a continuous predictor variable. I was wondering how to get the pairwise odds ratios for comparisons of levels of a categorical predictor variable. I can't seem to get the correct output using: > sp.glm=glm(cohort$logreg~cohort$hab, family=binomial) > summary(sp.glm) This gives me the coefficient for the effect of hab overall, but not the coefficients for the comparisons of each of the levels. I suspect this has something to do with the contrasts statement, but I can't figure it out, and would be...

Buenos días, Necesito calcular la vida media de cada uno de nuestros modelos de producto pero tengo un histórico incompleto (para los productos más nuevos tengo la vida completa de las cohortes, pero para los más antiguos tengo sólo la parte final de la vida de las cohortes). He decidido calcular la esperanza de vida dinámica con un modelo actuarial, concretamente la librería 'demography' y tengo varias dudas sobre su uso: 1- Todas las funciones permiten el uso por 'c...

Hi. I've a dataframe with about 46000 women with about 500 cases (cancers). I want to define a case-cohort sampling scheme, matching by age and hospital centre. Is there anyone who has already written a code for that? It should be something similar to the stcacoh macro in Stata. TIA, Stefano -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read ht...