similar to: Cox model with interaction

Dear all, I'm working with R 1.2.2 on win98 (256Mb RAM) with ESS 5.1.18. It happens that the start-up is very slow. I don't know exactly how to measure, but with a simple watch I recorded 15 seconds, while the startup of S+ 2000 (Sqpe+4 under ESS) takes around 3 seconds. Tha fact is that sometimes it takes almost a minute (58.30 sec recorded) for the prompt to appear! And ESS gives me an

Hi all. I'm "trying" to install R on a linux box (mandrake 7.1) but when I do make install it stops with the following errors: /usr/bin/install: zero.so: File doesn't exist (translated fron italian!!) make[2]: *** [install-objects] Error 1 make[2]: *** [install] Error 2 make: *** [install] Error 1 Where I'm wrong? The only change I've done is ./configure

Hi all. A simple question. Is there a function to compute the Odds Ratio and its confidence intervall, from a logistic model (glm(.......,family=binomial....). I've written my own, but certainly someone did a better job. Thank you in advance, Stefano *********************************************** Stefano Calza Istituto di Statistica Medica e Biometria Universit? degli Studi di Milano Via

Hi all. I'm wondering why there isn't yet the debian version of R 1.2.3. Am I missing something? TIA, Stefano ******************************* Stefano Calza Istituto di Statistica Medica e Biometria Universit? degli Studi di Milano Via G. Venezian, 1 20133 Milano tel. +39 02/2361302, 02/70600908 fax 02/2362930

Dear all, how can I get the R-squared of the sequence of models fitted by step(....,direction="forward")? TIA, Stefano ******************************* Stefano Calza Istituto di Statistica Medica e Biometria Universit? degli Studi di Milano Via G. Venezian, 1 20133 Milano tel. +39 02/2361302, 02/70600908 fax 02/2362930

Hi, is there a function/library like print.display by F.Harrel for R? BTW Harrel clearly states in the reame that "you may not port the library to R" ... Thanks, Stefano. -- *********************************************** Stefano Calza Istituto di Statistica Medica e Biometria Universit? degli Studi di Milano Via G. Venezian, 1 20133 Milano tel. +39 02/2361302, 02/70600908 fax

I am interested Andrea is whether you ever established why your R2 was 1. I have had a similar situation previously. My main issue though, which I'd be v grateful for advice on, is why I am obtaining such negative values -0.3 for Somers Dxy using validate.cph from the Design package given my value of Nagelkerke R2 is not so low 13.2%. I have this output when fitting 6 variables all with

Hi, I have tried to install the R package Rcplex on windows xp without success. I have only cplex_studio124.win-x86-32 version. I have modified the makevars.win file as indicated in the installation guide ( http://cran.r-project.org/web/packages/Rcplex/INSTALL ), then I zipped the whole folder to order an installation from a zipped file (utils :menuInstallLocal()) The message I had is file

Dear all, a question about the package boot: is it possible to bootstrap a function of arguments of different length? is there an alternative package? below an example of a simple try, many thanks paolo > library(boot) > blue<-c(4, 3, 4, 1, 2, 3, 1) > red<-c(4, 3, 4, 1, 2, 3, 1, 2, 3, 1, 3, 2) > > test1<-function(x,y,d){ + return(sum(x) - sum(y)) + } > b=boot(red

Hi Cristian, I don't know about a R routine for exact CIs, but I found a function called "ageadjust" using a gamma distribution approximation some time ago -take a look at http://medepi.org/epitools/rfunctions/index.html. That Webpage seems a bit outdated now, but the calculations given in the ageadjust functions should still work, I think. Maybe this helps; regards Heinrich.

Hi! How can I fill with colors a portion of a graph (e.g.: I want fill in red the area within two confidence intervals)? Thank you very much Cristian ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Cristian Pattaro ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ Unit of Epidemiology & Medical Statistics University of Verona Tel +39 045 8027668 fax +39 045 505357 cristian at biometria.univr.it

A collegue of mine, in a different office branch of my organization, suffer continuous troubles on roaming profiles, troubles thta in my site, with rather similar setup, does not happens. Trouble last some month, but 'till now we have thinked was the network, and we had worked checking cables, switches, configurations. But trouble continues... and seems frequent but not with a specific

Dear all, I'm using the package "Genetics", and I'm interested in the computation of D' statistics for Linkage Disequilibrium, for which the LD() command has been realised. Unfortunately I don't find any reference on "how" the D' is computed by the LD() function. In the package documentation it is generally referred as "MLE" estimation, but

Hi! After installing R 2.0.0, I have problem with the "help.search" as below: > help.search("table") Error in help.search("table") : could not find package 'file23525' > help.search("mean") Error in help.search("mean") : could not find package 'file23525' I have had the same problem with other request... What is

Dear all, I need to read an ASCII file with diffent length lines. This is what is contained in the file gene.txt: 1st line ID description snp_id genotype 2nd line 10003 Low rs152240 3rd line 10003 Moderate rs189011 TC 4th line 10004 Conservative rs152240 GC 5th line 10004 Bad rs154354 6th line 10013 Bad rs152240 7th line 10019 Conservative rs152240 AC etc... This is what I would like to obtain

Dear all, I have a surprising problem with the representation of frequencies in a histogram. Consider, for example, the R code: b<-rnorm(2000,3.5,0.3) hist(b,freq=F) When I plotted the histogram, I expected that values in the y-axis (the probability) varied between 0 and 1. Instead, they varied within the range 0-1.3. Have you got any suggestion for obtaining a correct graph with

Dear Thierry, thank you for the quick reply. I have only one question about the approach you proposed. As you suggested, imagine that the model we end up after the model selection procedure is: mod2.1 <- lmer(dT_purs ~ T + Z + (1 +T+Z| subject), data =x, REML= FALSE) According to the common procedures specified in many manuals and recent papers, if I want to compute the p_values relative to

Dear all, is there a automatic method to obtain Odds Ratio estimates and their confidence intervals from a GLM Logistic model? Thanks Cristian ============================================= Cristian Pattaro ============================================= Unit of Epidemiology & Medical Statistics Department of Medicines and Public Health University of Verona cristian@biometria.univr.it

Hi All, I was wondering if someone could help me to solve this issue with lmer. In order to understand the best mixed effects model to fit my data, I compared the following options according to the procedures specified in many papers (i.e. Baayen <http://www.google.it/url?sa=t&rct=j&q=&esrc=s&source=web&cd=1&ved=0CDsQFjAA

Hi all! Do you know if there is any R function/package that can be used to estimate "logit" models with panel data and forecasting? Thanks, Moysés. -- Moysés Nascimento Bacharel em Estatística/UFES Mestre em Estatística Aplicada e Biometria/UFV Doutorando em Estatística e Experimentação Agropecuária/UFLA moysesnascim@gmail.com [[alternative HTML version deleted]]