similar to: Function comparable to cutpt.coxph from "Survival Analysis using S"

The title says it all really; I am looking for a function along the lines of cutpt.coxph as described in "Survival Analysis Using S" (Tableman and Kim), Chapter 6. As may be guessed, the function optimises the cutpoint of a continuous variable for cox proportional hazard modelling. I can't find it, or any similar function, on CRAN. Alternatively, perhaps there is a way of extracting

Hi Everybody, please can you help me how to call BayesX in R in order to see the graph already exist in BayesX Thanks ---------- Forwarded message ---------- From: <r-help-request@r-project.org> Date: Sat, Jan 22, 2011 at 5:00 AM Subject: R-help Digest, Vol 95, Issue 22 To: r-help@r-project.org Send R-help mailing list submissions to r-help@r-project.org To subscribe or

Hi everyone, Does anyone know if there is a package to compute raking weights using R? What I need is to create a variable with weights base in some demographic variables (e.g. sexo, age group, area) using the raking procedure. Thank you in advance! -- Sebasti?n Daza

Hi everyone, Does anyone know how to get "black and white theme" (grey scale,, I would say) graphs using lattice or ggplot2, as it is shown in this webpage: http://lmdvr.r-forge.r-project.org/figures/figures.html? I am using Sweave, and I cannot get that color configuration. I have added the following option: trellis.device(color=FALSE) but I got a pdf file with color graphs. Thank

Hello, I have got a cluster sample using an election dataset where I already had the final results of a county-specific election. I am trying to figure out what would be the best sampling design for my data. The structure of the dataset is: 1) polling station (in general schools where people vote, for a county, for example, there are 15 polling stations) 2) inside each polling station, there

Hi everyone, I have just got different samples from a dataframe (independent and exclusive, there aren't common elements among them). I want to create a variable that indicate the sampling selection of the elements in the original dataframe (for example, 0 = no selected, 1= sample 1, 2=sample 2, etc.). I have tried to do it with ifelse command, but the problem is that the second line

Hello everyone, I am new in R and I am trying to create a id number (a correlative sequence of numbers) by group, and a correlative sequence of numbers inside each group (my idea is to get statistics by group without having to aggregate the database). Here an example: group id_group correlative_group A 1 1 A 1 2 A 1 3 A 1 4 B 2 1 B 2 2 B 2 3 C 3 1 C 3 2 C 3 3 C 3 4 C 3

Dear List, Does anybody no how to compare mean survival times for two (more) groups in R? What test statistics should I use? Thank you very much! Joe [[alternative HTML version deleted]]

Hi everyone, I have a dataset of friendship with this format: ego alter 4746 1 2 9742 1 3 14738 1 NA 4747 2 NA 9743 2 3 14739 2 1 4748 3 13 9744 3 5 14740 3 14 4749 4 NA 9745 4 NA 14741 4 NA 4750 5 NA 9746 5 13 14742 5 10 4751 6 12 9747 6 7 ... NA means that individuals don't select any

Dear R users, I still feel new to R so please apologize if I am doing something stupid here. My use of the polygon() function produces a result that I cannot comprehend: In a plot, I would like to shade the area below a normal distribution. However, I do not want the entire area to be shaded, but just the area on the right side of a vertical line that I draw through the distribution (in

Hi everyone, I am trying to build a table putting standard errors horizontally. I haven't been able to do it. library(memisc) berkeley <- aggregate(Table(Admit,Freq)~.,data=UCBAdmissions) berk0 <- glm(cbind(Admitted,Rejected)~1,data=berkeley,family="binomial") berk1 <- glm(cbind(Admitted,Rejected)~Gender,data=berkeley,family="binomial") berk2 <-

Skipped content of type multipart/mixed-------------- next part -------------- A non-text attachment was scrubbed... Name: not available Type: application/pgp-signature Size: 189 bytes Desc: Digital signature Url : http://westfish.xiph.org/pipermail/vorbis-dev/attachments/20040618/e7bee431/attachment.pgp

Hi everyone, I am just trying to figure out how to do a xyplot where in addition to dots and lines I can change dots' colors according to an individual variable (e.g., marital disruption across time, a dummy 0/1). When I use "groups" specification (see below), I get two different lines for each individual based on groups, and what I want is to get one line connecting dots, and

Hi everyone, I am trying to build a table putting standard errors horizontally. I haven't been able to do it. library(memisc) berkeley <- aggregate(Table(Admit,Freq)~.,data=UCBAdmissions) berk0 <- glm(cbind(Admitted,Rejected)~1,data=berkeley,family="binomial") berk1 <- glm(cbind(Admitted,Rejected)~Gender,data=berkeley,family="binomial") berk2 <-

Hi everyone, I need to get a between-component variance (e.g. random effects Anova), but using lmer I don't get the same results (variance component) than using random effects Anova. I am using a database of students, clustered on schools (there is not the same number of students by school). According to the ICC1 command, the interclass correlation is .44 > ICC1(anova1) [1] 0.4414491

How does one use the vcut from Vorbis-tools package? I tried to enter samples, miliseconds, seconds, h:mm:ss as cut point, but nothing worked (I even tried to enter 1 as cutpoint, but I always got "Cutpoint not within stream." message)... I'm doing this on Windows, if it matters... -- Jernej Simoncic, jernej.simoncic@guest.arnes.si http://www2.arnes.si/~sopjsimo/ ICQ: 26266467

Dear all, I'm trying to manage with user defined split function in rpart (file rpart\tests\usersplits.R in http://cran.r-project.org/src/contrib/rpart_3.1-34.tar.gz - see bottom of the email). Suppose to have the following data.frame (note that x's values are already sorted) > D y x 1 7 0.428 2 3 0.876 3 1 1.467 4 6 1.492 5 3 1.703 6 4 2.406 7 8 2.628 8 6 2.879 9 5 3.025 10 3 3.494

Dear List: I have searched the archives and my R books and cannot find a method to transform a continuous variable into a binary variable. For example, I have test score data along a continuous scale. I want to create a new variable in my dataset that is 1=above a cutpoint (or passed the test) and 0=otherwise. My instinct tells me that this will require a combination of the transform

<mailto:r-help@r-project.org>Hi every one, I want to calculate the combination function in R, the value not all the possible choices. I mean cmbn(5,2)=10. Is there any function unless using factorial? Regards, Amir [[alternative HTML version deleted]]

Dear all, I have two variables, y and x. It seems that the relationship between them is Piecewise Linear Functions. The cutpoint is 20. That is, when x<20, there is a linear relationship between y and x; while x>=20, there is another different linear relationship between them. How can i specify their relationships in R correctly? # glm(y~I(x<20)+I(x>=20),family = binomial, data =