similar to: cox memory

  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 =

Hi everyone, I was trying to mode the following exercise using R. The question: Set up a one pool model using numericintegration. The model will run from time 1 to time 30 using a time step of 1.The pool (A) will be fed by flux "inA" at a rate of 5 units per hour anddrained by flux "outA" at a rate of 20% per hour. At time 0, A has 5units. At time 30, what is the pool size of

Hello, I have a cohort with approx 1,200 patients at the ages of 30-65 that had their first myocardial infarction during 1992: · They were in a follow up until 2005. · About 400 of them died during this period of time (right censored) · Each one of them had up to 4 mi recurrent events. I am using the semi-parametric model in order to assess the relationship of

Hi, when i run the following code i get this massege: "The instruction at 0x######## reference 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,

--- include/share/alloc.h | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/include/share/alloc.h b/include/share/alloc.h index 914de9ba..63878db0 100644 --- a/include/share/alloc.h +++ b/include/share/alloc.h @@ -168,7 +168,7 @@ static inline void *safe_realloc_add_2op_(void *ptr, size_t size1, size_t size2) free(ptr); return 0; } - return realloc(ptr, size2); + return

Dear R experts, I am trying to arrange multiple plots, creating one graph for each size1 factor variable in my data frame, and each plot has the median price on the y-axis and the size2 on the x-axis grouped by clarity: library(ggplot2) df <- data.frame(price=matrix(sample(1:1000, 100, replace = TRUE), ncol = 1)) df$size1 = 1:nrow(df) df$size1 = cut(df$size1, breaks=11)

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

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

Hi, I would like to get from a plot the size of the symbols plotted. Imagine I have the following plot function : plot(1:2,1:2, pch=15, cex=4) I would like the get the values SIZE1 and SIZE2 so that if I plot the following rectangle : rect(1.5,1.5, 1.5+SIZE1, 1.5+SIZE2) then the size of this square is exactely the same as the one of the symbols that have been plotted. Thanks for any idea. --

tftp-hpa-0.35 and tftp-hpa-0.36 show wrong transmission time. >tftp -V tftp-hpa 0.36, with readline >tftp -m octet -v localhost -c get aaa /tmp/aaa Connected to localhost.localdomain (127.0.0.1), port 69 getting from localhost.localdomain:aaa to /tmp/aaa [octet] Received 33554432 bytes in -3.3 seconds [-81407477 bit/s] ^^^^ ^^^^^^^^^ By the way, why

Hi, Does anyone know of a function in R that can handle a fractional variable as the dependent variable? The catch is that the function has to be inclusive of 0 and 1, which betareg() does not. It seems like GLM might be able to handle the fractional logit model, but I can't figure it out. How do you format GLM to do so? Best, Rachael [[alternative HTML version deleted]]

Hi , I have a data set with recurrence time (up to four) of myocardial infarction (MI). Part of the file is showing below: Num1 Trt Sex Time T1 T2 T3 T4 1011 1 1 9 1211 0 1 59 3020 1 2 14 3 1245 0 1 18 12 16 3069 1 2 26 6 12 13 2051 0 1 53 3 15 46 51 The data consist of the following eight variables: Num1 , patient number Trt, treatment group (1=placebo and 2=drug) Sex,

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

Here is our proposal to extend/enhance the x86-64 compact unwind descriptors to fully describe the prologue/epilogue for asynchronous unwinding.  I believe there are missing/lacking CFI directives as well, but I'll save that for another thread. Asynchronous Compact Unwind Descriptors Ron Brender, VMS Software, Inc. Revised January 25, 2018 1  Introduction This document proposes means to

Hi Nick, It is a pleasure to be in contact with the creator of the compact unwind approach! I can see how an array of 32-bit unwind blocks could be used to describe each distinct point within a function (within a prolog in particular). But then you end up with six or seven or more such blocks for a large percentage of functions, don't you? Seems like a lot of additional space for something

I have a matrix with 3 years of data (2006, 2009, 2010). I am trying to split this matrix by year so that I have 3 separate matrices. My matrix looks like this:     Q16.1 Year Gender Grade1       3 2006      1     52       2 2006      0     53       3 2006      0     54       3 2006      0     55       3 2006      1     56       2 2006      0     57       3 2006      0     58       3 2006      1  

John and Ron, I developed the original compact unwind implementation for macOS 10.6 back in 2009. I tried to leave space in the design to support finer grain exception handling such as for asynchronous or for the shrink wrap optimization. The idea I had at the time was instead of having just one 32-bit compact unwind info per function, there could be an array of them each covering a different

Hi, Sorry not to provide R-code in my previous mail. R code is below #install.packages("rms") require(rms) #install.packages("mstate") library(mstate) require(splines) library(ggplot2) library(survival) library(splines) #install.packages("survsim") require(survsim) set.seed(10) df<-crisk.sim(n=500, foltime=10, dist.ev=rep("lnorm",2),

I'm looking at an issue reported by the Coverity static analyzer. In iconvert() in src/share/utf8/iconvert.c on line 152 there is newbuf = safe_realloc_add_2op_(utfbuf, ...); If the request size is not valid, the function will free utfbuf and return 0. This is followed by goto fail and utfbuf is freed for the second time. A simply fix would be to set utfbuf to 0 if newbuf is 0. However, this

I haven't done any R memory debugging lately, but https://www.mail-archive.com/rcpp-devel at lists.r-forge.r-project.org/msg10289.html shows how I used to have gdb break where valgrind finds a problem so you could examine the details. Also, running your code after running gctorture(TRUE) can help track down memory problems. -Bill On Wed, Feb 7, 2024 at 12:03?PM Therneau, Terry M., Ph.D.