similar to: censor=FALSE and id options in survfit.coxph

Hi, You may try: dat1 <- structure(list(SubID = 1:8, CSE1 = c(6L, 6L, 5L, 5L, 5L, 5L, 3L, 3L), CSE2 = c(5L, 4L, 5L, 4L, 6L, 4L, 6L, 6L), CSE3 = c(6L, 7L, 5L, 3L, 7L, 3L, 6L, 6L), CSE4 = c(2L, 2L, 5L, 4L, 5L, 6L, 3L, 3L), WSE1 = c(6L, 6L, 5L, 4L, 6L, 4L, 6L, 6L), WSE2 = c(2L, 6L, 5L, 4L, 4L, 3L, 5L, 5L), WSE3 = c(2L, 2L, 4L, 5L, 4L, 7L, 2L, 4L), WSE4 = c(4L, 3L, 5L, 2L, 1L, 3L, 1L, 7L)),

Hello R Help! I would like to make a legible boxplot of tree growth rates for each of seven tree species at each of seven different sites. It's a lot of data to put on one figure, I know. I made a beautiful, interpretable figure using color, but my target journal can't deal with color figures. I can use seven shades of grey to fill the boxes, but the figure then becomes uninterpretable -

Hi, I am not sure I understand your question correctly. dat1<- read.table(text=" id??????????? responsed_at???????????????? number_of_connection????????????????? scores 1????????????????? 12-01-2010?????????????????????????????????? 1????????????????????????????????????????????? 2 1????????????????? 15-02-2010??????????????????????????????????

HI, Try this: dat1<- read.table(text=" ?V1,V2,V3,V4,V5,V6,V7 ?chr1,564563,564598,564588 564589,1336,+,134 ?chr1,564620,564649,564644 564645,94,+,10 ?chr1,565369,565404,565371 565372,217,+,8 ?chr1,565463,565541,565480 565481,1214,+,15 ?chr1,565653,565697,565662 565663,1031,+,28 ?chr1,565861,565922,565883 565884,316,+,12 ",sep=",",header=TRUE,stringsAsFactors=FALSE)

Hello, I am trying to filter a data set like below so that the peaks in the Phase value are more obvious and can be identified by a peak finding function following the useful advise of Carl Witthoft. I have written the following for(i in length(data$Phase)){ newphase=if(abs(data$Phase[i+1]-data$Phase[i])>6){ data$Phase[i+1] }else{data$Phase[i] } } I get the following error which I have not

Hello, Thank you to all those great folks that have helped me in the past (especially Dennis Murphy). I have a new challenge. I often generate time-series data sets that look like the one below, with a variable ("Phase") which has a series of flat-topped peaks (sample data below with 5 "peaks"). I would like to calculate the phase value for each peak. It would be great to

Hello, I have a column called max_date in my data frame and I only want to keep the bigger values for the same activity. How can I do that? data frame: activity max_dt A 2013-03-05 B 2013-03-28 A 2013-03-28 C 2013-03-28 B 2013-03-01 Thank you for your help -- View this message in context:

Hi everyone: This may be trivial but I just have not been able to figure it out. Imagine the following dataframe: a b c d TRUE TRUE TRUE TRUE FALSE FALSE FALSE TRUE FALSE TRUE FALSE FALSE I would like to create a new dataframe, in which TRUE gets 0 but if false then add 1 to the cell to the left. So the results for the example above should be something like: a b c

Hi, I'm trying to copy the first row of one data frame to another. This is the statement I am using : df2[1,]<-df1[1,]; I have printed them out separately: df1[1,] = A C D E F But after copying: df2[1,] = 96 29 88 122 68 Why isn't it copying? They are both data frames, and "as.character" isn't working either. Thanks for your input :) [[alternative HTML version

Hi,. We have 2 values (first formant F1, second formant F2) for a given phoneme for six languages. We want to see whether the languages are significantly different one from another for this given phoneme. We have done a manova on our data and it works well, but we doesn't allow us to see which pair of languages are different. If we have only one formant for the phoneme, we would use

Thank you very much for your reply. Your code work well with this example. I modified a little to fit my real data, I got an error massage. Error in split.default(x = seq_len(nrow(x)), f = f, drop = drop, ...) : Group length is 0 but data length > 0 On Thu, Jan 31, 2013 at 12:21 PM, arun kirshna [via R] < ml-node+s789695n4657196h87@n4.nabble.com> wrote: > Hi, > Try this: >

The following fit does not make sense to me, please, correct me if I have a logical error. > moddowsn Call: coxph(formula = Surv(start, stop, resp) ~ sn + matfac2, data = coxsn1, method = "efron") coef exp(coef) se(coef) z p sn2 0.0497 1.051 0.02030 2.450 1.4e-02 sn3 -0.0532 0.948 0.02038 -2.610 9.0e-03 sn4 -0.0410 0.960

Dear Dr Skaug and R users, just discovered glmm.admb in R, and it seems a very useful tool. However, I ran into a problem when I compare two models: m1<-glmm.admb(survival~light*species*damage, random=~1, group="table", data=bm, family="binomial", link="logit") m1.1<-glmm.admb(survival~(light+species+damage)^2, random=~1, group="table", data=bm,

Dear Mao Jianfeng, "r-help-owner" is not the place for help, but: r-help at r-project.org (CC-ed here) In any case, strsplit() does the job, i.e.: > unlist(strsplit("BCPy01-01", "-")) [1] "BCPy01" "01" You can work with the whole variable, like: splitpop <- strsplit(df1$popcode, "-") then access the first part with >

I have a model with quant vars only and the error message does not make sense: (mod1 <- coxph(Surv(time=strt,time2=stp,event=(resp==1))~ +incpost+I(amt/1e5)+rate+strata(termfac), subset=dt<"2010-08-30", data=inc,method="efron")) Call: coxph(formula = Surv(time = strt, time2 = stp, event = (resp == 1)) ~ +incpost + I(amt/1e+05) + rate + strata(termfac),

Hi r-helpers. Please forgive my ignorance, but I would like to plot a smoothing spline (smooth.spline) from package "stats", and show the knots in the plot, and I can't seem to figure out where smooth.spline has located the knots (when I use nknots). Unfortunately, I don't know a lot about splines, but I know that they provide me an easy way to estimate the location of local

Could you please help me on the following ANCOVA issue? This is a part of my dataset: sampling dist h 1 wi 200 0.8687212 2 wi 200 0.8812909 3 wi 200 0.8267464 4 wi 0 0.8554508 5 wi 0 0.9506721 6 wi 0 0.8112781 7 wi 400 0.8687212 8 wi 400 0.8414646 9 wi 400 0.7601675 10 wi 900 0.6577048 11 wi 900

Hi, I've been trying to work with the gnls() function in the "nlme" package. My decision to use gnls() was so that I could fit varPower and such to some of the data. However, in working with a small dataset, I've found that the results given by gnls() don't seem to make any sense and they differ substantially from those produced by nls(). I suspect that I am just

Hello R helpers, I would like to automate this code for many files of the same type. But I don´t know how to make it. In particular, i don´t know how to read many files each one as an r object with the name of the file. Then a for loop would be sufficient, right? Many thanks and a happy new year. Dominic datos <- read.table('global2001.csv',head=T,sep=';',stringsAsFactors=F)

Hi all, I'm seeing this failure on my PPC G4 box running TOT with llvm-gcc 4.2. Is anyone else seeing this? I'm sure it's related to the byval stuff that's recently gone into LLVM. I'm attaching the output of this command: $ llvm-gcc -emit-llvm -O3 -S -o - -emit-llvm /Users/wendling/llvm/ llvm.src/test/CFrontend/2008-01-25-ByValReadNone.c As you can see in it, there