similar to: Error using newdata argument in survfit

Dear List, Is the test for proportional hazards valid when the model contains a cluster variable? The output looks strange with the cluster variable. My intervals are based on calendar time and the clustering variable is related to the season the event occurs in. model1<-coxph(Surv(Start,Stop,Event)~LagAOO+I(LagAOO^2)+ FirstSeen+TSLE+strata(CumPOO.rc)+cluster(quarter),data=data8, x=T)

Full_Name: Hendrik Weisser Version: 2.6.1 OS: Linux Submission from: (NULL) (139.19.102.218) The following computation never finishes and locks R up: > values <- list(x=10:30, y=c(23.85, 28.805, 28.195, 26.23, 25.005, 20.475, 17.33, 14.97, 11.765, 8.857, 5.3725, 5.16, 4.2105, 2.929, 2.174, 1.25, 1.0255, 0.612, 0.556, 0.4025, 0.173)) > y.max <- max(values$y) > model <- nls(y ~

day<-structure(c(7.7, 7.7, 7.7, 7.7, 7.7, 7.71, 7.7, 7.71, 7.71, 7.7, 7.7, 7.7, 7.7, 7.69, 7.68, 7.68, 7.67, 7.67, 7.67, 7.66, 7.65, 7.65, 7.65, 7.64, 7.64, 7.63, 7.63, 7.63, 7.62, 7.62, 7.62, 7.62, 7.63, 7.63, 7.63, 7.63, 7.63, 7.64, 7.64, 7.65, 7.65, 7.65, 7.66, 7.66, 7.67, 7.67, 7.67, 7.68, 7.68, 7.69, 7.69, 7.69, 7.69, 7.7, 7.7, 7.7, 7.7, 7.7, 7.71, 7.7, 7.7, 7.71, 7.71, 7.7, 7.7, 7.7, 7.7,

#I would like to interpolate a straight line between 06/08/06 04:16:00 - 06/08/06 20:31:00 with values and then sum them. This is an estimate of ecosystem #respiration and I will be using this in a larger context(48 days of these diurnal curves), but for right now I am just trying to figure out how to do it for this one #day example. I have some other code for Ecosystem (stream) Metabolism that

dear professor: I am a doctor of urinary,and I am developing a nomogram of bladder tumor.Now I have a problem about this. I have got the result like this through analysing the dataset "exp11.sav" through multinominal logistic regression by SPSS 17.0.(the Sig. is high,that is good ,it is just aexperimental data ) Parameter Estimates Ya B Std. Error Wald df Sig. Exp(B) 95%

R Help - I'm trying to run a correlation matrix with a covariate of "age" and will at some point will also want to covary other variables concurrently. I'm using the "psych" package and have tried other methods such as writing a loop to extract semi-partial correlations, but it does not seem to be working. How can I accomplish this? library(psych) > set.cor(y =

Hello all, Here goes one of my first functions. I want to make a nonparametric multiple sample comparison with unequal sample sizes (see Zar?s Biostatistical Analysis, 3rd. Ed., pg. 201 Example 10.11, pg. 288 Example 11.10). In the real world, I want to compare samples of fish length captured with different fishing gears. After using the Kruskal-Wallis test I want to check the differences

Hi all This is probably a simple problem but somehow I am having much trouble with finding a solution, so I seek your help! I have a data-set with continuous response variables. The explanatory variably is 4xpH treatments (so 8.08, 7.94, 7.81 and 7.71) so also continuous and not technically factorial. However I have decided to do Anova's (as well as regression) to explore the effect of

Hi all, I have this problem with lattice that xyplot() won't draw some of my axis labels if the type (i.e. the relation argument) of scales is set as free. For example, in the plot below, I would want it to also show: 1. the labels E1,...E6 below the 10th panel (i.e. 3rd row, 2 col)....just as it is now done below the 12th panel.... 2. as well as the labels (2,4,6,8) on the top of panels 1

I am struggling to generate p values for comparisons of levels (post-hoc tests) in a glm with a negative binomial distribution I am trying to compare cell counts on different days as grown on different media (e.g. types of cryogel) so I have 2 explanatory variables (Day and Cryogel), which are both factors, and an over-dispersed count variable (number of cells) as the response. I know that both

Hello, I am doing simulations, and I generate a list at each iteration (with three component matrices in the example below), saving the results in a list. For example, after two iterations, I have something like > str(sim.theta) List of 2 $ :List of 3 ..$ : num [1:6, 1:4] -3.67 -1.07 -2.99 -18.38 -3.26 ... ..$ : num [1:6, 1:6] -7.56 -3.14 -4.99 1.03 2.79 ... ..$ : num [1:6, 1:4]

I tried toy problems and there doesn't seem to be a basic problem between lm() and ts objects: X = data.frame(x=c(1,2,7,9), y=c(7,2,3,1)) lm(y ~ x, X) X <- lapply(X, function(x) ts(x, frequency=12, start=c(1994,7))) lm(y ~ x, X) and this works fine - whether you do an lm() before or after making ts objects, it's okay. But I have a situation where things aren't okay.

Hola eRReros. Os lo explico con un ejemplo: # Cargamos los paquetes y el ejemplo install.packages("SPEI") library("SPEI") data(wichita) # los primeros 12 datos head(wichita,12) # mi subset de los primeros 12 datos meu<-wichita[c(1:12),] meu # como veis los valores de TMED son iguales en ambos dataframes. # ahora viene el problema # calculamos

Thanks Don and Jim, "Get" did the trick! That command is new to me, and is exactly what this problem needed. The syntax that worked is actually slightly simpler than what was suggested. This is what worked: all.comps <- ls(pattern="^res") for(i in all.comps){ obj = get(i); ... Thanks again. --Kelly V. -----Original Message----- From: MacQueen, Don [mailto:macqueen1 at

Hi folks, On Fri, Apr 24, 2015 at 08:01:54AM +0800, Jones Syue wrote: > Hello Ervin, > > > I'll notify to you about the result. (First, I need to upgrade > > that server.) > > Sure! That will be great help me to understand if this patch > works on thast compiling case too, > thank you. I"ve downloaded the Ubuntu's Samba 4.1.6 deb source, and patch it

X2R is a bundle of three software libraries for passing complicated data structures from Fortran, C/C++, or AD Model Builder to R. An update has been sent to CRAN and should be available from mirrors shortly. From the menu at the left of the CRAN home page, look under Software / Other. We also have updated FishGraph, a compatible set of R functions for examining output from fish population

X2R is a bundle of three software libraries for passing complicated data structures from Fortran, C/C++, or AD Model Builder to R. An update has been sent to CRAN and should be available from mirrors shortly. From the menu at the left of the CRAN home page, look under Software / Other. We also have updated FishGraph, a compatible set of R functions for examining output from fish population

X2R is a bundle of three software libraries allowing the user to pass structured data easily from Fortran, C/C++, or AD Model Builder to R. An update to X2R has been sent to CRAN and should be available at mirrors shortly. We also have uploaded FishGraph, a set of R functions to generate graphics from fish population models. From the menu at the left of the CRAN home page, look under Software /

X2R is a bundle of three software libraries allowing the user to pass structured data easily from Fortran, C/C++, or AD Model Builder to R. An update to X2R has been sent to CRAN and should be available at mirrors shortly. We also have uploaded FishGraph, a set of R functions to generate graphics from fish population models. From the menu at the left of the CRAN home page, look under Software /

Hello all, I would like to connect two IPv4-only networks through global IPv6-only network with tinc. Could anyone tell me how to configure tinc.conf and the other config files? Thanks in advance, Ichiro Tinc: Discussion list about the tinc VPN daemon Archive: http://mail.nl.linux.org/lists/ Tinc site: http://tinc.nl.linux.org/