search for: parrinello

...f 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 Statistics Department of Biotecnologies Viale Europa, 11 25100 Brescia Italy Tel: +3930/3717528 Fax: +3930/3717488 [[alternative HTML version deleted]]

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

...tistical 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. Giovanni Parrinello External Lecturer Medical Statistics Unit Department of Biomedical Sciences Viale Europa, 11 - 25123 Brescia Italy Tel: +390303717528 Fax: +390303717488 email: parrinel@med.unibs.it [[alternative HTML version deleted]]

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 nam...

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 nam...

...- c ------- I am looking for functions in R to estimate the sample size and/or to conduct the analysis. I have just found an article from Byar in Statistics in Medicine for a 2^3 incomplete factorial design, but I would like not to discover again the wheel.. TIA dr. Giovanni Parrinello Section of Medical Statistics Department of Biosciences University of Brescia 25127 Viale Europa, 11 Brescia Italy Tel: +390303717528 Fax: +390303701157 [[alternative HTML version deleted]]

Dear All, I am looking for the library QVF [ NONLINEAR MEASUREMENT ERROR PROGRAMS FOR SAS AND SPLUS ]. TIA Giovanni Parrinello dr. Giovanni Parrinello Medical Statistics Unit Department of Medical Sciences University of Brescia Via Valsabbina, 19 25127 Brescia Italy tel.: +39303717528 fax:+39303701157 E-mail: parrinel at med.unibs.it url: http://www.med.unibs.it/dip/stat/index.html -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-....

Dear All, I would like to find something ( references, code,..) to implement a comparison of three treatments in an observational study using the 'Propensity Score'. Any help is much appreciated. Thanks! Giovanni -- dr. Giovanni Parrinello Department of Biotecnologies Medical Statistics Unit University of Brescia Viale Europa, 11 25123 Brescia email: parrinel at med.unibs.it Phone: +390303717528 Fax: +390303717488

...#39;) --- Please select a CRAN mirror for use in this session --- Error in as.list(read.dcf(file = file)[1, ]) : indice fuori limite # index out of bound TIA Giovanni PS: I have already tried to solve the problem reinstalling R or changing the options of update.packages.. -- dr. Giovanni Parrinello Department of Biotecnologies Medical Statistics Unit University of Brescia Viale Europa, 11 25123 Brescia email: parrinel at med.unibs.it Phone: +390303717528 Fax: +390303717488

Dear All, can anyone explain this message Error in load("matched2.RData") : Value of SET_STRING_ELT() must be a 'CHARSXP' not a 'builtin'? Have I the possibility to retrieve the data? TIA Giovanni -- dr. Giovanni Parrinello Department of Biotecnologies Medical Statistics Unit University of Brescia Viale Europa, 11 25123 Brescia email: parrinel at med.unibs.it Phone: +390303717528 Fax: +390303717488

...acchetto richiesto: 'survival' The following object(s) are masked from package:Hmisc : untangle.specials Carico il pacchetto richiesto: class Carico il pacchetto richiesto: nnet ######### So I have many unrequired packages loaded. Any idea? TIA Giovanni dr. Giovanni Parrinello Section of Medical Statistics Department of Biotecnologies Viale Europa, 11 25100 Brescia Italy Tel: +3930/3717528 Fax: +3930/3717488 [[alternative HTML version deleted]]

...:tools" [5] "package:stats" "package:graphics" "package:grDevices" "package:utils" [9] "package:datasets" "package:methods" "Autoloads" "package:base" R version 2.7.2 -- dr. Giovanni Parrinello Department of Biotecnologies Medical Statistics Unit University of Brescia Viale Europa, 11 25123 Brescia email: parrinel at med.unibs.it Phone: +390303717528 Fax: +390303717488

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

Hello to All. I'd want to use a one-way ANOVA. This means that I have only one factor, with, lets say, 5 levels. I made a dataframe, called "DATA", with two Columns: A, that is my response, and it is "class numerical". B, that defines the different levels of my factor, and it is "class factor". If I want to use a fixed effect model, I know that the formula I have

Dear All, I am looking an R version of the "Computer program for sample size and power calculations in the design of multi-arm and factorial clinical trials with survival endpoints". Best regards, Giovanni Parrinello P.S.: in the meantime I am preparing a summary for my preceeding question about time-varying covariates in the Cox model and I thank Frank Harrell, Chuck Cleland, Patrick Connolly and Perperoglou A for their interesting answers. -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-...

Dear All, I have found in the poly help this sentence: The orthogonal polynomial is summarized by the coefficients, which can be used to evaluate it via the three-term recursion given in Kennedy & Gentle (1980, pp. 343–4), and used in the predict part of the code. My question: which type of orthogonal polynomials are used by this function? Hrmite, legendre.. TIA Giovanni [[alternative HTML

channel <- odbcConnectExcel("nuova tabella terapia occupazionale mod.xls") > ## list the spreadsheets > sqlTables(channel) TABLE_CAT TABLE_SCHEM TABLE_NAME 1 c:\\TABELLE DEFINITIVE\\nuova tabella terapia occupazionale mod <NA> 'emi tot 2006 OAI_60g TO FU1$' 2 c:\\TABELLE DEFINITIVE\\nuova

...e issues raised are still of interest... Heather >>> Prof Brian Ripley <ripley at stats.ox.ac.uk> 01/18/06 01:31pm >>> On Wed, 18 Jan 2006, Heather Turner wrote: [Lines wrapped for legibility and R-help removed as it is not an appropriate list.] > 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 d...

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, reading an article by Rodolphe Thiebaut and Helene Jacqmin-Gadda ("Mixed models for longitudinal left-censored repeated measures") I have found this program in SAS proc nlmixed data=TEST QTOL=1E-6; parms sigsq1=0.44 ro=0.09 sigsq2=0.07 sigsqe=0.18 alpha=3.08 beta=0.43; bounds $B!](B1< ro < 1, sigsq1 sigsq2 sigsqe >= 0; pi=2*arsin(1); mu=alpha+beta*TIME+a i+b i*TIME;