similar to: SAS code in R

Dear All, I have analysed time to event data for continuous variables by considering the multivariable fractional polynomial (MFP) model and comparing this to the untransformed and log transformed model to determine which transformation, if any, is best. This was possible as the Cox model was the underlying model. However, I am now at the situation where the assumption that the competing risks

I am having trouble formatting some survival data to use in a time dependent cox model. My time dep. variable is habitat and I have it recorded for every day (with some NAs). I think it is working properly except for calculating the death.time. This column should be 1s or 0s and as I have it only produces 0s. Any help will be greatly appreciated.

Hola a todos, -Estoy estudiando el efecto de dos genotipos (~tratamientos) en la aparición de síndrome metabólico (MetS) con datos longitudinales recogidos a tiempo 0,7,10,15,20 y 25 años. -He hecho un dataframe con las siguientes variables MetS: Síndrome Metabólico (Si=1,No=0) bmi: Indice de masa corporal (IMC) cuando se produce la conversión a MetS+ . Para los que permancen MetS-, esta variable

Dear R-help, I'm very sorry to ask 2 questions in a week. I am using the package 'crr' and it does exactly what I need it to when I use the dataset a. However, when I use dataset b I get the following error message: Error in drop(.Call("La_dgesv", a, as.matrix(b), tol, PACKAGE = "base")) : system is computationally singular: reciprocal condition number =

Dear R users, My question is more methodology related rather than specific to R usage. Using time on study as time in a cox model, eg: library(Design) stanf.cph1=cph(Surv(time, status) ~ t5+id+age, data=stanford2, surv=T) #In this case the 1000-day survival probability would be: stanf.surv1=survest(stanf.cph1, times=1000) #Age in this case is a covariate. #I now want to compare the above

Hi, I have the binned data (observed and generated from model)  that I would like to test using the anderson-darling goodness of fit test.  But I'm not sure which package in R to use. I tried ad.test(...) but it does not recognise the test by Vito Ricci in FITTING DISTRIBUTIONS WITH R   > ad.test(hist_hume_beec[,1],hist_hume_beec[,2]) Error: could not find function "ad.test"

Dear useRs, There is a new version 0.1-2 of the package surv2sample available on CRAN. Users of the previous versions should update because a bug in the function cif2.ks has been fixed. General information about the package: surv2sample provides various two-sample tests for right-censored survival data. Three main areas and corresponding methods are: * comparison of two survival

Dear useRs, There is a new version 0.1-2 of the package surv2sample available on CRAN. Users of the previous versions should update because a bug in the function cif2.ks has been fixed. General information about the package: surv2sample provides various two-sample tests for right-censored survival data. Three main areas and corresponding methods are: * comparison of two survival

You asked about survival curves with age scale versus follow-up scale. > fit1 <- coxph(Surv(time/365.25, status) ~ t5 + id + age, data=stanford2) > surv1<- survfit(fit1) > surv1 n events median 0.95LCL 0.95UCL 157.000 102.000 1.999 0.898 3.608 > summary(surv1, times=3) time n.risk n.event survival std.err lower 95% CI upper 95% CI 3 46 85 0.451

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.

?I've hit a roadblock debugging a new update to the survival package.?? I do debugging in a developement envinment, i.e. I don't create and load a package but rather? source all the .R files and dyn.load an .so file, which makes things a bit easier. ? Running with R -d "valgrind --tool=memcheck --leak-check=full" one of my test files crashes in simple R code a dozen lines

Hi all, I'm planning to work through the book "Analyzing Linguistic Data" by R.H. Baayen, which is an introduction to R used for, well, what the title says. ;-) On the first page of the book, Baayen says that in order to work with the book, R needs to download and install a number of packages from CRAN. The problem is that one of these packages, "Design", has apparently

Hi, I am trying to run a conditional logistic model on a nested case-control study using cph() and then estimate survival based on the model. The data came from Prof Bryan Langholz website where he also has the SAS code to this, so I am trying to replicate the SAS results. The data attached. Basically, the variables are: rstime: risk set age rsentry: fake entry time, just before rstime setno:

The combination of survfit, coxph, and factors is getting confused. It is not smart enough to match a new data frame that contains a numeric for sitenew to a fit that contained that variable as a factor. (Perhaps it should be smart enough to at least die gracefully -- but it's not). The simple solution is to not use factors. site1 <- 1*(coxsnps$sitenew==1) site2 <-

Hi all, is it possible to configure dovecot in that way that special folders do not create automatically by some imap clients? i.e: outlook or thunderbird creates a folder junk-mail i would like that such folders are not created automatically or by users Best regards Stefan

Can anyone offer suggestion about getting readline to work on Solaris? (v1.5.1) When I use the up arrow or ctrl-p to try to repeat the last command my session just beeps. It appears that the .Rhistory file that should appear in my startup directory is not being created. Thanks! The screen shot below shows cled=TRUE and no .Rhistory file created: $ /apps/R-1.5.1/bin/R R : Copyright 2002, The

Hi @ll, i use an intel S5000VSA server board with two xeon cpus. I installed a debian squeeze amd64 system with debian xen and qemu packages. Before i installed xen i check out successfully the virtualization flag of cpus by grep "vmx" /proc/cpuinfo After installation of xen (http://wiki.debian.org/Xen#Installation_on_squeeze) i can not run a hvm domu. xm dmesg | grep -i vmx shows

# doveconf -a # 2.2.13: /etc/dovecot/dovecot.conf doveconf: Fatal: Error in configuration file /etc/dovecot/conf.d/10-ssl.conf line 6: Unknown setting: ssl doveconf: Error: managesieve-login: dump-capability process returned 89 doveconf: Fatal: Error in configuration file /etc/dovecot/conf.d/10-master.conf line 35: Invalid size: $default_vsz_limit # grep default_vsz_limit

Hi Stefan, My apologies for the long delay. Been ill. I tried this and had an error, # mv /etc/dovecot/conf.d/10-ssl.conf /etc/dovecot/conf.d/10-ssl.conf.bak # cp -p /usr/share/dovecot/conf.d/10-ssl.conf /etc/dovecot/conf.d/ # /etc/init.d/dovecot restart /var/log/daemon.log Jun 4 11:05:55 mail systemd[1]: Starting Dovecot IMAP/POP3 email server... Jun 4 11:05:55 mail systemd[1]: Started

Hi Aki, # dovecot --build-options Build options: ioloop=epoll notify=inotify ipv6 openssl io_block_size=8192 Mail storages: shared mdbox sdbox maildir mbox cydir imapc pop3c raw fail SQL driver plugins: mysql postgresql sqlite Passdb: checkpassword ldap pam passwd passwd-file shadow sql Userdb: checkpassword ldap(plugin) nss passwd prefetch passwd-file sql -- On June 5, 2017 8:59:08 AM