similar to: Random Effects for One-Way Anova

He, I used the "coxph" function, with four covariates. Let's say something like that > model.1 <- coxph(Surv(Time,Event)~X1+X2+X3+X4,data=DATA) So I obtain the 4 coefficients B1,B2,B3,B4 such that h(t) = h0(t) exp(B1*X1+ B2*X2 + B3*X3 + B4*X4). When I use the function on the same data > predict.coxph(model.1,type="lp") how it works in making the prediction?

Ciao a tutti ho appena iniziato ad utilizzare R per ora per attuare un'analisi geostatistica di dati. Volevo sapere come poter creare un oggetto gstat partendo da un file testo(che ho gia importato con read.table)e che contiene 3 colonne: x,y,value. Mi servirebbe far questo per costruire un variogramma. So che la domanda molto probabilmente per voi sara' banale....vi ringrazio

Salve a tutti; Sono un neofita di Ruby, sinceramente Rails ancora non lo ho utilizzato. Sono uno studente universitario, ora sto preparando una tesina su un framework per sviluppare applicazioni web non so se lo conoscete e lo Spring webflow. nella tesina ho cercato di convincere il prof sul fatto che Rails ? migliore. Solo che lui mi ha fatto una domanda che mi ha spiazzato alla quale al

Dear all, I am having quite a hard time in trying to figure out how to correctly spell out a model in R (a repeated-measures anova with a within-subject covariate, I guess). Even though I have read in the posting guide that statistical advice may or may not get an answer on this list, I decided to try it anyway, hoping not to incur in somebody's ire for misusing the tool. For the sake of

Hi all I am running a 2-way repeated measure anova with 1 between-subjects factor (Group=treatment, control), and 1 within-subject factor (Time of measurement: time1, time2). I extract the results of the anova with: summary(aov(effect ~ Group*Time + Error=Subj/Time, data=mydata)) Now, this must be clearly a dumb question, but how can I quickly extract in R all the post-hoc t-tests for the

Dear All, discussing with a statistician of a pharmaceutical company I received this answer about the statistical package that I have planned to use: As R is not a validated software package, we would like to ask if it would rather be possible for you to use SAS, SPSS or another approved statistical software system. Could someone suggest me a 'polite' answer? TIA Giovanni -- dr.

Dear All, I am looking for an implementation in R of the Appell Hypergeometric function. Any suggestions will be more than appreciated! GP -- dr. Giovanni Parrinello External Lecturer Medical Statistics Unit Department of Biomedical Sciences Viale Europa, 11 - 25123 Brescia Italy Tel: +390303717528 Fax: +390303717488 email: parrinel at med.unibs.it

Hello, I am planning a study with the main point to evaluate the interaction of two treatments, but for ethical reasons one cell is empty, that with patients receaving no treatment at all Treatment B

Hi all I am running a 2-way repeated measure anova with 1 between-subjects factor (Group=treatment, control), and 1 within-subject factor (Time of measurement: time1, time2). I extract the results of the anova with: summary(aov(effect ~ Group*Time + Error=Subj/Time, data=mydata)) Now, this must be clearly a dumb question, but how can I quickly extract in R all the post-hoc t-tests for the

Dear All, I got these results from the example in the function "dbeta": >x <- seq(0, 1, length=21) > dbeta(x, 1, 1) [1] 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 Any Idea? TIA Giovanni dr. Giovanni Parrinello Section of Medical Statistics Department of Biosciences University of Brescia 25127 Viale Europa, 11 Brescia Italy Tel: +390303717528 Fax: +390303701157

Last week Giovanni Parrinello posted a message asking why various packages were loaded when he loaded an .Rdata file. Brian Ripley replied saying he thought it was because the saved workspace contained a reference to the namespace of ipred. (Correspondence copied below). This begs the question: how did the reference to the namespace of ipred come to be in the .Rdata file? Brian did say it is

Good morning, i'm drawing a barplot (data are stored in a matrix,nrow=3,ncol=10). So i get 10 groups each containing 3 sub-bars. I need to change the range of the y-axis when i plot respectively the first and the other two sub-bars in each group. I can't manage. Could you help me? Thank you. Angela

Dear All, I am trying to apply the function 'cox.zph' of the library survival, but I receive this error message: not found the object 'residuals.coxph'. I have re-installed the library 'survival' without any change and also a search with RSiteSearch was unsuccessful.. Any suggestion? TIA Giovanni PS: R 2.2.0 Windows XP HE dr. Giovanni Parrinello Section of Medical

Dear R Fellow-Travellers: What is your recommended way of dealing with a left-censored response (non-detects) in (linear Gaussian) mixed effects models? Specifics: Response is a numeric positive measurement (of volume, actually); but when it falls below some unknown and slightly random value (depending on how the sample is prepared and measured), it cannot be measured and is recorded as 0.

Dear all, I am unsure about how to specify a model in R and I thought of asking some advice to the list. I have two groups ("Group"= A, B) of subjects, with each subject undertaking a test before and after a certain treatment ("Time"= pre, post). Additionally, I want to enter the age of the subject as a covariate (the performance on the test is affected by age),

Thanks Denis, I was looking for the option 'dns:x' in the wiki but I didn't find it. Now it works. I used    log level = 3 auth:3  dns:0 auth_audit:3 gives me unknown class message But where I can find a complete list of classes for log level? I'll also give a try on the last version of samba with json. Thanks again Giuseppe On 1/18/2018 4:52 PM, Denis Cardon wrote:

Ciao, the rot filesystem of my thumper is a ZFS with a single disk: bash-3.2# zpool status rpool pool: rpool state: ONLINE scrub: none requested config: NAME STATE READ WRITE CKSUM rpool ONLINE 0 0 0 c5t0d0s0 ONLINE 0 0 0 spares c0t7d0 AVAIL c1t6d0 AVAIL c1t7d0

dear all, This is a probably a silly question. If I type > grep("x",c("a.x" ,"b.x","a.xx"),value=TRUE) [1] "a.x" "b.x" "a.xx" Instead, I would like to obtain only "a.x" "b.x" How is it possible to get this result with grep()? many thanks for your attention, best, vito --

On 18/06/2014 04:17, Charles Lepple wrote: > On Jun 13, 2014, at 2:53 AM, Giuseppe Corbelli <giuseppe.corbelli at copanitalia.com> wrote: > >> As said in previous mail, I just finished a first working version of a driver for the UPS found on ASEM PB1300 device >> (http://www.asem.it/prodotti/industrial-automation/box-pcs/performance/pb1300/) >> Linux only, accessed

Dear all, It appears that MASS::polr() and Design::lrm() return the same point estimates but different st.errs when fitting proportional odds models, grade<-c(4,4,2,4,3,2,3,1,3,3,2,2,3,3,2,4,2,4,5,2,1,4,1,2,5,3,4,2,2,1) score<-c(525,533,545,582,581,576,572,609,559,543,576,525,574,582,574,471,595, 557,557,584,599,517,649,584,463,591,488,563,553,549) library(MASS) library(Design)