similar to: Questions concerning survival analysis

Hello everybody, does anybody know how the function plot.survfit exactly works? I'd like to plot the log of the cummulative hazard against the log time by using plot.survfit(...fun="cloglog") which does not work correctly. The scales are wrong and there is an error message about infinit numbers. It must have something to do with the censored data, doesn't it? #Example:

... and the mailing list is picky about attachments... whatever you attached did not conform to the stringent requirements mentioned in the Posting Guide. Pasting the code right into the email is usually safest, though you DO have to post using plain text (as the Posting Guide indicates) or your code may get mangled by the automatic html format removal. On May 15, 2018 7:04:31 AM PDT, Bert Gunter

Unless there is good reason not to, always cc the list -- there are lots of smarter folks than I on it who can help. I may or may not have time to look at this. Hopefully someone else will. -- Bert Bert Gunter "The trouble with having an open mind is that people keep coming along and sticking things into it." -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip

OK, Let's try this again! Here is the reproducible script; it is long because I had to copy the panel dataset here. My question is related to systemfit; I don't know how to get the result for the entire panel. #Reproducible script Empdata<- read.csv("/Users/ngwinuiazenui/Documents/UPLOADemp.csv") View(Empdata) install.packages("systemfit")

Hi r-users,   I would like to extract the data that match.  Attached is my data: I'm interested in matchind the value in column 'intg' with value in column 'rand_no' > cbind(z=z,intg=dd,rand_no = rr)             z  intg rand_no    [1,]  0.00 0.000   0.001    [2,]  0.01 0.000   0.002    [3,]  0.02 0.000   0.002    [4,]  0.03 0.000   0.003    [5,]  0.04 0.000   0.003    [6,] 

Dear R-People I did some two way variance analysis and have now the following problem: I''d like to do some multiple Comparisons between different factor levels (two factors, each with some 5 or 6 levels), because there are some significant interactions between these two factors. In S-plus 5, there is the function "multicomp", but in R, there ist no matchable funktion, as far

Hello I'd like to try out some functions in the package randomForest. Therefore, I did install this package. However, it is not possible to load the library, although I have R-Version 2.1.1 (i.e. later than 2.0.0). The commands I used and the Answers/Error from R is as follows: > install.packages("C://Programme//R//rw2011//library//randomForest_4.5-16.zip",

Sadly you failed to set your email program to send plain text and the data is corrupted at my end. I also think you need to reduce the size of the data set... the intent here is to increase your understanding, not debug your particular analysis. I will say that I am having a very challenging time understanding what you are trying to accomplish though. What are the equations that you think need

Dear all,I am a new user of WinBUGS and need your help. After running the following code, I got parameters of beta0 through beta4 (stats, density), but I don't know how to get the prediction of the last value of h, the variable I set to NA and want to model it using the following code.Does anyone can given me a hint? Any advice would be greatly appreciated.Best

Hello all, I am working with a 1000x1000 matrix, and I would like to return a 1000x1000 matrix that tells me which value in the matrix is greater than a theshold value (1 or 0 indicator). i have tried mat2<-as.matrix(as.numeric(mat1>0.25)) but that returns a 1:100000 matrix. I have also tried for loops, but they are grossly inefficient. THanks for all your help in advance. Lanre

First I would really like to thank the mailing list for help I got in the past, as a new to R I am really needing some support on hoe to code the following problem. I am trying to sort some data I have in a big file. The file has 4 columns and 19000 rows. An example of it looks like this:- G 0.892 A 0.108 G 0.883 T 0.117 T 0.5 C

I would like to use a tclArray: mytkarray <- tclArray() as the variable for a table: table1 <- tkwidget(f1, 'table', variable= mytkarray) but if I include character strings with embedded spaces, I get braces appearing in the table. I can remove them using a non-R tclarray, (the difference between the first example of http://www.sciviews.org/_rgui/tcltk/Tktable.html and

I'm getting an error message when using the new version of sqldf, > library(sqldf) > str(kdv) 'data.frame': 71 obs. of 3 variables: $ dpss: num 0.117 0.144 0.164 0.166 0.165 ... $ npdp: num 0.1264 0.0325 0.0109 0.0033 0.0055 ... $ logk: num 1.12 1.29 1.41 1.41 1.42 ... > test=sqldf("select * from kdv") Error in get("fun", env = this, inherits =

Dear all, I have been running some tests of my package RSiena on different platforms and trying to reconcile the results. With Mac, the commands options(digits=4) round(1.81652, digits=4) print 1.817 With Windows, the same commands print 1.816 I am not bothered which answer I get, but it would be nice if they were the same. A linux box agreed with the Mac. Mac sessionInfo(): R version

Hi. I'm writing an article on network backups, and instead of using my old ssh1 software, I decided to go with openssh all the way. I got the hang of the openssh way of doing protocol 2 public key authentication, but ssh is failing to terminate when a pipe is broken. I am ssh-ing to a remote host and doing a cat or zcat of a dump file, then on the localhost, I'm using restore to extract

Is there a way to have printing data.frames with POSIXct to display milliseconds if digits.secs is set as a default? You can use the digits argument in print, such as print(df, digits = 3) to get the intended output, but I assumed it was done with the option digits.secs set. Tibbles by default do this printing, which is shown below, but I was unsure if digits.secs should affect printing

hello: my data looks like: time1  time2   event  catagoria 2004    2006        1            C 2004    2005        0            C 2005    2010        1            E 2007    2009        1            C 2006    2007        0            E 2008    2010        0            C 2008    2010        1            E ... and the census interval is 1 year I have tried  this

After upgrading our CentOS 5 box to the latest revisions last week (including Dovecot 1.1), we're seeing the following error message in the log files. Sieve was working fine with Dovecot 1.0. I have yet to turn up anything via Google for this particular error. SELinux is not logging any error messages at the moment, so I'm pretty sure that we've properly allowed all the SELinux

No reason to tell me about security, I just need chmod for my tests. Where to change it ? On Thu, Jan 25, 2018 at 12:00 PM, <dovecot-request at dovecot.org> wrote: > Send dovecot mailing list submissions to > dovecot at dovecot.org > > To subscribe or unsubscribe via the World Wide Web, visit > https://dovecot.org/mailman/listinfo/dovecot > or, via email,

Hola, Estoy de vacaciones hasta el 1 de Agosto. Para dar soporte sobre la centralita de telefonia: voip at mildmac.es Perdonen las molestias. Ruth Llaneza Lapausa - Tecnico de VoIP. ruth.llaneza at mildmac.es Tlf: 902 199 384 Mildmac SA ? www.mildmac.es ? mildmac at mildmac.es C/ Hnos. Garc?a Noblejas 41, 6? planta. 28037 - Madrid Tlf: +34 91 501 33 02 Fax: +34 91 501 57 45