search for: censor

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. This...

-----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 IP>;branch=z9hG4bK484dcd1fDD872ECE From: "Michael" <sip:3109@<CENSORED>>;tag=D2DAE96E-D8618427 To: <sip:0100@<CENSORED>;user=phone> CSeq:...

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)~...

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, ...) locfit.censor(formula, ...) locfit(formula, ..., lfproc=locfit.censor) What I did in locfit.censor.Rd is have something like: \synopsis{ locfit.censor(x, y, cens, ..., iter=3, km=FALSE) } \usage{ locfit(formul...

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 wa...

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 P...

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...

...t (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: Connection refused - connect(2) Info: Retrieving plugin Info: Loading facts in /var/lib/puppet/lib/facter/root_home.rb Info: Loading facts in /var/lib/puppet/lib/facter/puppet_vardir.rb Info: Loading facts in /var/lib/puppet/lib/facte...

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...

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...

...sterisk-users] Polycom instant messages -----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 IP>;branch=z9hG4bK484dcd1fDD872ECE From: "Michael" <sip:3109@<CENSORED>>;tag=D2DAE96E-D8618427 To: <sip:0100@<CENSORED>;user=phone> CSeq:...

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 c...

Can anybody give me a neat example of interval censored data analysis codes in R? Given that suvreg(Surv(c(1,1,NA,3),c(2,NA,2,3),type="interval2")~1) works why does survreg(Surv(data[,1],data[,2],type="interval2")~1) not work where data is : T.1 T.2 Status 1 0.0000000 0.62873036 1 2 0.000...

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 censored type. However, whenever I set event to 2 o...

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 e...

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, 5.73...

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...

Hi , I am trying to get some estimator based on lognormal distribution when we have left,interval, and right censored data. Since, there is now avalible pakage in R can help me in this, I had to write my own code using Newton Raphson method which requires first and second derivative of log likelihood but my problem after runing the code is the estimators were too high. with this email ,I provide simple example f...

...D 2 11 F D 2 12 F D 2 13 F D 2 14 M A 3 15 M A 3 16 M A 3 17 M B 3 18 M B 3 19 M B 3 20 M A 4 Here is the reshaping I am seeking (explanation below). id mygroup mytime censor [1,] 1 A 1 1 [2,] 2 B 3 0 [3,] 2 C 3 0 [4,] 2 D 6 1 [5,] 3 A 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, t...