search for: censoring

Dear R-experts: I am using survreg() to estimate the parameters of a Weibull density having right-censored observations. Some observations are weighted. To do that I regress the weighed observations against a column of ones. When I enter the data as 37 weighted observations, the parameter estimates are exactly the same as when I enter the data as the corresponding 70 unweighted observations.

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Is it possible to use the instant messaging feature of Polycom phones in Asterisk? At the moment I'm seeing this in the SIP messaging when I try to send one from a Polycom 450. <--- SIP read from UDP:<CENSORED POLYCOM IP>:5060 ---> INVITE sip:0100@<CENSORED>:5060;user=phone SIP/2.0 Via: SIP/2.0/UDP <CENSORED POLYCOM

Dear r-helpers, This is my first time to run survival analysis. Currently, I have a data set which contains two variables, the variable of time to event (or time to censoring) and the variable of censor indicator. For the indicator variable, it was coded as 0 and 1. 0 represents right censor, 1 means event of interest. Now I try to use "survfit" in the package of "survival". I wrote the following code: > rptsurv <- survfit(surv(time,censor)~1,d...

Dear R-devel, I'm working on Prof. Loader's new version of locfit to try to get it pass R CMD check. I'm almost there, but I have a problem with some Rd files that I hope some one can help me resolve. Here's an example: In the package there's a function called locfit.censor(). This function can be used in a few different ways: locfit.censor(x, y, cens, ...)

I have the complete data like id time censor 1 10 0 1 20 0 1 30 0 2 10 0 2 20 1 2 30 0 2 40 0 3 10 0 3 20 0 3 30 1 .... for id 1, i want to select the last row since all censor indicator is 0; for id 2, i want to select the row where censor ==1; for id 3, i also want to select the row where censor==1. So if there is a 1 for censor, then I want to select such a row, otherwise I want to select the

Dear all, Are there functions or packages in R that can handle interval-censored data? I have looked in various packages (such as survival5 or event), but it seems that only right-censored data can be analysed. More generally, are there methods to analyse both interval-censored observations and right-censored observations in the same data set? Thanks in advance. Emmanuel Paradis

Can survreg() handle interval-censored data like the documentation says? I ask because the command: survreg(Surv(start, stop, event) ~ 1, data = heart) fails with the error message Invalid survival type yet the documentation for Surv() states: "Presently, the only methods allowing interval censored data are the parametric models computed by 'survreg'"

I auto-upgraded puppet (as usual) using "yum upgrade -y", and now all agents are getting error 400. # puppet --version 3.0.0 On agent (any!) # puppet agent --test Warning: Unable to fetch my node definition, but the agent run will continue: Warning: Error 400 on SERVER: Could not retrieve facts for build6.[censored]: Failed to find facts from PuppetDB at abyss.[censored]:8081:

Hello, I have interval censored data, censored between (0, 100). I used the tobit function in the AER package which in turn backs on survreg. Actually I'm struggling with the distribution. Data is asymmetrically distributed, so first choice would be a Weibull distribution. Unfortunately the Weibull doesn't allow for zero values in time data, as it requires x > 0. So I tried the

Terry, I recently noticed the censor argument of survfit. For some analyses it greatly reduces the size of the resulting object, which is a nice feature. However, when combined with the id argument, only 1 prediction is made. Predictions can be made individually but I'd prefer to do them all at once if that change can be made. Chris ##################################### # CODE # create

I didn't know that instant messaging was a feature with the Polycom's. Do you have any documentation, how-to's, etc. that you can point me to? That would be just way too cool. Thanks; John -----Original Message----- From: asterisk-users-bounces at lists.digium.com [mailto:asterisk-users-bounces at lists.digium.com] On Behalf Of Michael Englehorn Sent: Monday, January 12, 2015 12:19

In your particular case I don't think that censoring is an issue, at least not for the reason that you discuss. The basic censoring assumption in the Cox model is that subjects who are censored have the same future risk as those who were a. not censored and b. have the same covariates. The real problem with informative censoring are the cova...

...'survreg'" > > Any pointers as to what I'm missing? > You are specifying left-truncated, right-censored data, not interval-censored. You need the type= argument, eg > Surv(c(1,1,NA,3),c(2,NA,2,3),type="interval2") [1] [1, 2] 1+ 2- 3 specifies interval censoring at [1,2], right-censoring at 1, left-censoring at 2, and an observed event at 3. -thomas

Dear R users, I tried to analysis the hazard function of some data by coxph function in survival package. The type of the data include "left-censored", "right-censored", "both right-censored and left-censored" (btw, does this has a technical term?), and "complete" ones. I noticed that event (one parameter in "Surv()") might be an indicator for the

One set of data has censored (less-than detection limits) water chemistry concentrations for 80-100% of all observations. My initial trial-and-error attempts to apply the cenboxplot() method suggests that it has an upper limit to the percentage of censored observations. I do not see this limit in Dennis Helsel's second edition. Has anyone experience plotting censored data and can provide

Hello, I would like to know if there is a function in R that will test for normality and handle censored data sets. Currently, I evaluate each censored data set by the extent to which a normal scores plot approximate a straight line. For complete data sets I use shapiro.test(). Below is an example of a censored data set. data1<-c(0.00, 0.00, 0.00, 5.86, 5.17, 8.17, 5.12, 4.92, 7.08,

Hi, I need to get mean and variance of right censored data. How can I do that? I have a vector of values (called a) and vector of booleans (whether value is censored) (called b). What to do with this? Sorry, I'm R beginner. Thank you! Pete

Hi, I have a right-censored (positive) random variable (e.g. failure times subject to right censoring) that is observed for N subjects: Y_i, I = 1, 2, ..., N. Note that Y_i = min(T_i, C_i), where T_i is the true failure time and C_i is the censored time. Let us assume that C_i is independent of T_i. Now, I would like to generate another random variable U_i, I = 1, 2, ..., N, which is correlated...

...on, I got problem with estimation, where it was too high" overestimation". Is there any one cen help me to know my problem when can be or recommend me to read some good references that can be useful? Best regards.  Below is the example,   cur=curd=cen=cens=array(1,100) Cur=RealCur=realcensoring=realcured=array(1,20) ExpCure=Bias=RealCure=array(1,21) ################################## Z1=c(rbinom(100,1,0.5)) Z2=c(rbinom(100,1,0.5)) dat1<-data.frame(time=rlnorm( 100,2,0.8),Censored=rbinom(100,1,0.9),Cured=rbinom(100,1,.3)) dat2<-dat1[order(dat1[,1]),] # order the dat...

...3 0 [6,] 3 B 3 1 [7,] 4 A 1 1 I need to create 2 variables. The first one is a time variable. Observe that for id=2, the variable mygroup=B was observed 3 times. In the solution we see in row 2 that id=2 has a mytime variable of 3. Next, I need to create a censoring variable. Notice id=2 goes through has values of B, C, D for mygroup. This means the change from B to C and C to D is observed. There is no change from D. I need to indicate this with a 'censoring' variable. So B and C would have values 0, and D would have a value of 1. As another example...