similar to: Case-Cohort Analysis

Dear all, I have a couple of questions regarding the survival:::cch function. 1) I notice that Prentice and Self-Prentice functions are giving identical standard errors (not by chance but by programming design) while their estimates are different. My guess is they are both using the standard error form from Self and Prentice (1986). I see that standard errors for both methods are

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 different. Can we

Dear all, I am reading the cch function source code. But I can not understand the following codes. Please help me. What's the fitter here? fitter <- get(method) out <- fitter(tenter = tenter, texit = texit, cc = cc, id = id, X = X, ntot = nn, robust = robust) [[alternative HTML version deleted]]

I am interested in developing a backend for the PIC microcontrollers. Specifically, I plan to write for the 8-bit MCUs, but would hopefully also support the 16 and 32 bit models (eventually). It is my understanding that a backend for PIC16 was previously a part of the project but was dropped in version 2.9. Is there any development on this front currently? Is there any interest in developing this

I am attempting to fit a full zero-inflated Poisson model then use backward elimination to arrive at the best-fitting model. When I try to fit the model with zeroinfl() I get this error: Error in while (abs((ll_old - ll_new)/ll_old) > control$reltol) { : missing value where TRUE/FALSE needed Any suggestions? Thanks for your help! Paige Barlow MS Student Virginia Tech Dept Fish

Hi, I am reading a "|" delimited text file into R using read.table(). I am using colClasses= to specify some variables as factors. Some of these variables include missing values coded as "NA". Unfortunately the R code I am using (pasted bellow) includes "NA" as one of the factor levels. Is it possible to remove the "NA" level from a factor with in

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

? stato filtrato un testo allegato il cui set di caratteri non era indicato... Nome: non disponibile Url: https://stat.ethz.ch/pipermail/r-help/attachments/20061124/6ac3d11f/attachment.pl

Hello Gurus: When I plot a survival fit using rpart for the classification tree, for each node, there is a decimal based number above the event/total. I tried to see if it's the exponential ratio or logrithmics, neither seem to be the case. I'm wondering if anyone knows what they are. Thanks, Karen _________________________________________________________________ Find hidden words,

Case cohort function cch() is in survival package. In cch(), the prentice method is implemented like this: Prentice <- function(tenter, texit, cc, id, X, ntot,robust){ eps <- 0.00000001 cens <- as.numeric(cc>0) # Censorship indicators subcoh <- as.numeric(cc<2) # Subcohort indicators ## Calculate Prentice estimate ent2 <- tenter ent2[cc==2] <-

Dear R users I would like to calculate the Cumulative incidence for an event adjusting for competing risks and adjusting for covariates. One way to do this in R is to use the cmprsk package, function crr. This uses the Fine & Gray regression model. However, a simpler and more classical approach would be to implement the Kalbfleisch & Prentice method (1980, p 169), where one fits cause

survexp in survival5 doesn't seem to work for me. see below: > library(survival5) Attaching Package "package:survival5": The following object(s) are masked from package:base : sort.list > library(chron) > data(ratetables) > survexp(~ratetable(year=julian(6,1,1991), + sex=1,age=35*365.24),times=(0:30)/6*365.24) Error in as.character(as.date(c(min(R[, 3]),

I am just beginning to learn R. Please help. I have two matrices. One 20 columns by 16 rows represent student responses on a 20 question test. The second is 20 by 16 representing the correct responses 16 times. Subtracting one from the other should yield 0's if a response is correct. Yet Response minus Bigkey doesn't work. Documents indicate vectors should add and subtract, and

Dear All, A method for the estimation is univariate unimodal densities (with unknown mode) is described in "Statistical Inference under Order Restrictions" by Barlow et al.. Would anyone know whether there is an R-implementation (preferably with reference) for the estimation of univariate discrete unimodal densities (with unknown mode)? Thanks in advance for your help. Kind

Hi, everyone. I'm new to the list. This question may have come up before, but I couldn't find an answer to it by searching the centos website. Could be I was using bad keywords for my search. I have 3 HP xw9300 workstations (dual AMD Opteron 64-bit processors) that I'm upgrading from RHEL3WS to CentOS 4.4. I've upgraded two already, with one left to go. On the system with

I have a study best described as a retrospective case-cohort design: the cases were all the events in a given time span surveyed, and the controls (event-free during the follow-up period) were selected in 2:1 ratio (2 controls per case). The sampling frequency for the controls was about 0.27, so I used a weight vector consisting of 1 for cases and 1/0.27 for controls for coxph to adjust

Is anyone using CentOS 4.4 with an EMC Clariion SAN? In particular, I want to know if PowerPath and Naviagent will work on a CentOS 4 system connected to a san using Qlogic or Emulex cards. I tried a couple years ago with Fedora Core 1 or 2, and couldn't get the kernel modules from PowerPath to work since the kernel wasn't the official RHEL one. I'm most interested in

Hi, I am writing a basic function to extract the z scores for some linear regression coefficients: zscore<-function( y, x) { lm<-lm( y ~ x ) z <- coef(lm)/sqrt(diag(vcov(lm))) return(z) } I would like to pass a dataframe to the function as a argument so the function call changes from zscore(df$y1,df$x1) to zscore(y1,x1,data=df) but I am not sure how to reference the data

-----Original Message----- From: centos-bounces at centos.org [mailto:centos-bounces at centos.org] On Behalf Of Yiorgos Stamoulis Sent: Thursday, October 19, 2006 3:33 PM To: CentOS mailing list Subject: [CentOS] is everyone aware of this? BRUCE STANLEY wrote: > I am having some issues (e.g. System lock ups) with my ATI 9250 Video > card. > > I am thinging about replacing it with

I''m a new puppet user. I can''t puppetmasterd to start. I searched this group for the error message I get, but found nothing. Any ideas? [root@puppet ~]# puppetmasterd undefined method `name'' for nil:NilClass Using the --debug switch and ruby''s --debug option shows the following: [root@puppet ~]# puppetmasterd --debug debug: Parsing /etc/puppet/puppet.conf