Displaying 18 results from an estimated 18 matches for "tstop".
Did you mean:
stop
2009 Mar 10
2
simple question beginner
Hi there,
I am beginner in R and I have some basic question. Suppose I run a common procedure such as a t test or cox model like below:
out<-coxph( Surv(tstart,tstop, death1) ~ x1+x1:log(tstop+1) , test1,method=c("breslow"))
Which yields the following result:
Call:
coxph(formula = Surv(tstart, tstop, death1) ~ x1 + x1:log(tstop +
1), data = test1, method = c("breslow"))
coef exp(coef) se(coef) z p
x1 ...
2009 Sep 21
2
cox memory
Hi there,
I have a rather large data set and perform the following cox model:
test1 <- list(tstart,tstop,death1,chemo1,radio1,horm1)
out1<-coxph( Surv(tstart,tstop, death1) ~ chemo1+chemo1:log(tstop+1)+horm1+horm1:log(tstop+1)+age1+grade1+grade1:log(tstop+1)+positive1+positive1:log(tstop+1)+size1+size1:log(tstop+1), test1)
out1
Up to here everything works fine (with each covariate having a length...
2018 May 29
1
Difficulty in writing R code for one pool dynamic model
...nswer sheet provide the value of pool could be (23.8, 24.8, 25.0 or 24.7). Thevalue I got does not correspond to any of the choices given above. ?I don't know where I got it wrong.?
Below is the Ro code I tried to estimate thepool size of A at time 30.
#Step 1: Specify time?
Tstart <- 1
Tstop <- 30
Tstep <- 1
#Step 2: Set variables corresponding to pools
inA<-5
A <- inA
t <- Tstart
#Step 3: Set up rate parameters
A<-5
kinA <- 5
koutA <- 0.2
#Step 4: Write arrays that track pool size overtime
val <- array(1, dim=c(1, length(seq(Tstart,Tstop, by=Tst...
2010 Jun 21
1
Question about wine qcap v4l's MediaSampleTime
in wine/dlls/qcap/v4l.c:
ReadThread()
--> OutputPin_GetDeliveryBuffer((OutputPin *)capBox->pOut, &pSample, NULL, NULL, 0);
--> IMediaSample_SetTime(*ppSample, tStart, tStop);
in wine/dlls/quartz/memallocator.c:
StdMediaSample2_GetTime(IMediaSample2 * iface, REFERENCE_TIME * pStart, REFERENCE_TIME * pEnd)
--> if tstart/stop is NULL, leave pStart/pEnd not setted.
in wine/dlls/quartz/videorenderer.c:
VideoRenderer_Sample(LPVOID iface, IMediaSample * pSample)
-->...
2013 Nov 20
4
How to stop Kaplan-Meier curve at a time point
Hello R users
I have a question with Kaplan-Meier Curve with respect to my research. We
have done a retrospective study on fillings in the tooth and their survival
in relation to the many influencing factors. We had a long follow-up time
(upto 8yrs for some variables). However, we decided to stop the analysis at
the 6year follow up time, so that we can have uniform follow-up time for
all the
2004 Feb 18
0
[PATCH] fix tftp verbose mesage
...^^^^ ^^^^^^^^^
By the way, why don't you use timersub() or something like it?
--- tftp-hpa-0.36.org/tftp/tftp.c 2004-01-09 05:47:00.000000000 +0900
+++ tftp-hpa-0.36/tftp/tftp.c 2004-02-18 22:17:27.000000000 +0900
@@ -422,8 +422,8 @@
{
double delta;
- delta = (tstop.tv_sec+(tstop.tv_usec/100000.0)) -
- (tstart.tv_sec+(tstart.tv_usec/100000.0));
+ delta = (tstop.tv_sec+(tstop.tv_usec/1000000.0)) -
+ (tstart.tv_sec+(tstart.tv_usec/1000000.0));
if (verbose) {
printf("%s %lu bytes in %.1f seconds", direction, amount, delta);
printf(" [%....
2010 Nov 13
0
using if statment and loops to create data layout of recurrent events
...s:
ID, patient?s identification (which is the sequence number of the subject)
Visit, visit number (with value k for the kth potential MI recurrence)
TStart, time of the (k?1)th recurrence for Visit=k, or the entry time 0 if
VISIT=1, or the follow-up time if the (k?1)th recurrence does not occur
TStop, time of the kth recurrence if Visit=k or follow-up time if the kth
recurrence does not occur
Status, event status of TStop (1=recurrence and 0=censored)
For instance, patient # 3 with only one recurrence time at month 3 who was
followed until month 14 will have values for Visit, TStart, TStop,...
2010 Aug 09
2
recurrent events
...The exposure variable ("ESHKOL") is the socio-economic of the neighborhood
while adjusting for additional variables such as age, sex, etc. In addition
each patient has geo- statics -code for his neighborhood ("stat_unit").
I ran an Andersen Gill model:
Ag<-coxph(Surv(TStart,TStop,Status)~ ESHKOL+VAR1+
VAR2+cluster(id),data=recmi)
My questions are:
· I would like to check the unique effect of the "stat_unit"
variable. How can I write it in code?
· In your opinion does the AG model is the best model for this type
of data? Or should I use models th...
2010 Aug 19
1
memory problem
...memory at
0x#######, the memory cannot be "read".
and then i have to close R.
what is the problem and how can i solve it?
thanks in advance
Avi
my code
# frailtypack
library(frailtypack)
cgd.ag <- read.csv("C:/rfiles/RE/cgd.csv")
cgd.nfm <-frailtyPenal(Surv(TStart, TStop,
Status)~cluster(Center)+subcluster(ID)
Treatment,data=cgd.ag,Frailty=TRUE,n.knots=8,kappa1=50000,
cross.validation=TRUE,recurrentAG=TRUE)
--
View this message in context: http://r.789695.n4.nabble.com/memory-problem-tp2330510p2330510.html
Sent fro...
2018 Feb 16
0
Competing risks - calibration curve
...df$cause),0,df$cause)
table(df$cause)
df.w<-crprep("time","cause", data=df, trans=c(1,2), cens=0, id="nid", keep=c("x", "x.1"))
with(df.w,table(failcode,status))
ddist<-datadist(df.w)
options(datadist='ddist')
mod<-cph(Surv(Tstart,Tstop,status==1)~rcs(x,3)+x.1,data=df.w, weight=weight.cens, subset=failcode==1, x=T, y=T, surv=T, time.inc = 2.5) mod2<-cph(Surv(Tstart,Tstop,status==1)~(x.1+rcs(x,3))^2,data=df.w, weight=weight.cens, subset=failcode==1, x=T, y=T, surv=T, time.inc = 2.5)
mod2
########################## To develop n...
2013 Oct 09
1
frailtypack
I can't comment on frailtypack issues, but would like to mention that coxme will handle
nested models, contrary to the statement below that "frailtypack is perhaps the only ....
for nested survival data".
To reprise the original post's model
cgd.nfm <- coxme(Surv(Tstart, Tstop, Status) ~ Treatment + (1 | Center/ID), data=cgd.ag)
And a note to the poster-- you should reprise the original message to which you are
responding.
Terry Therneau
On 10/09/2013 05:00 AM, r-help-request at r-project.org wrote:
> Hello,
>
> I am encountering very similar problems with...
2024 Feb 07
2
Difficult debug
...dy allowed me to find a few other (mostly dumb)
> errors that led
> to an out of bounds access, e.g., the wrong endpoint variable in a for( )
> loop. What
> would others advise as a next step?
>
> Here is the last part of the screen
> > fit2 <- coxph(list(Surv(tstart, tstop, bstat) ~ 1,
> + c(1:4):5 ~ age / common + shared), id= id,
> istate=bili4,
> + data=pbc2, ties='breslow', x=TRUE)
> > surv2 <- survfit(fit2, newdata=list(age=50), p0=c(.4, .3, .2, .1, 0))
> > test2 <- mysurv(fit2, pbc2$bili4, p0...
2018 Apr 24
0
Compile issues when building survival package from source
...gw_32/bin/gcc -I"D:/R/R-35~1.0/include" -DNDEBUG -O3
-Wall -std=gnu99 -mtune=generic -c agfit4.c -o agfit4.o
agfit4.c: In function 'agfit4':
agfit4.c:334:20: warning: 'dtime' may be used uninitialized in this
function [-Wmaybe-uninitialized]
if (tstop[p] < dtime) break; /* no more to add */
When I build using R CMD build however, the file agfit4.c is never even
built. That explains the undefined references, but I have no clue what the
difference is between both methods.
I suspected a newer version so I downloaded the source from CRAN manual...
2006 Sep 07
0
counting process form of a cox model (cluster(id))?
Hi,
I am currently analysising a counting process form of a cox model allowing for the inclusion of time dependent covariates. An example model I have fitted is
modlqol<-coxph(Surv(Tstart,Tstop,cens.time)~tmt.first+risk +lqol+cluster(id),data=cat)
summary(modlqol)
My question is quick. I am looking at 1 event (death), and repeated measurements (the time dependent covariate 'lqol') are frequently taken on a subject, so I assume that measurements on the same subject will be correl...
2024 Feb 07
2
Difficult debug
...The
valgrind approach had already allowed me to find a few other (mostly dumb) errors that led
to an out of bounds access, e.g., the wrong endpoint variable in a for( ) loop.??? What
would others advise as a next step?
Here is the last part of the screen
> fit2 <- coxph(list(Surv(tstart, tstop, bstat) ~ 1,
+??????????????????? c(1:4):5 ~ age / common + shared), id= id, istate=bili4,
+?????????????? data=pbc2, ties='breslow', x=TRUE)
> surv2 <- survfit(fit2, newdata=list(age=50), p0=c(.4, .3, .2, .1, 0))
> test2 <- mysurv(fit2, pbc2$bili4, p0= 4:0/10, fit2, x0 =50)
=...
2011 Sep 17
0
Warning in 'probtrans'-function ('mstate'-package)
...),
status = c(NA, "s0", "s1", "s2"),
data= S1, trans=tmat, keep=covs)
# add 'transition-specific' covariates
SLV <- expand.covs(SL, covs, append=TRUE, longnames=TRUE)
head(SLV)
names(SLV)
# fit stratified Markov model
cox1 <- coxph(Surv(Tstart, Tstop, status) ~
(v1.1 + v2.1) * v3.1 +
(v1.2 + v2.2) * v3.2 +
(v1.3 + v2.3) * v3.3 +
strata(trans), SLV,
method = "breslow")
summary(cox1)
# 'cloned' data for making predictions: specific covariate 'v3' fixed
at 0 (& 1):
SLV0 <- SLV[1:3,]
SLV0$v3.1 <- 0
SLV0$v3.2...
2005 Feb 15
3
Autostart Asterisk on Slackware?
Maybe trivial question, but I cannot find an answer:
How to autostart Asterisk (daemon) on Slackware 10? I know that I should
put something in /etc/rc.d, but what?
2012 Feb 03
1
coxme with frailty--variance of random effect?
Dear all,
This probably stems from my lack of understanding of the model, but I
do not understand the variance of the random effect reported in coxme.
Consider the following toy example:
#------------------------------- BEGINNING OF CODE
------------------------------------------------
library(survival)
library(coxme)
#--- Generate toy data:
d <- data.frame(id = c(1:100), #