similar to: lm: mark sample used in estimation

Hi, perhaps this question was answered previously however I could not find them. My problem is how how to extract a particular statistic from the result given by lm(). For e.g. ctl <- c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14) trt <- c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69) group <- gl(2,10,20, labels=c("Ctl","Trt")) weight <- c(ctl, trt) >

Dear R-helpers, I was curious why the output of summary (and many other functions in R) does not add a separator between the name of a factor and the label of a level (e.g., in the example below, 'group Trt'). If the user had lower case labels (e.g., 'trt'), the output would be particularly hard to read. ctl <- c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14) trt

Hi! In R 3.0.0 from current SVN, ?summary.lm says: > Value [...] > df degrees of freedom, a 3-vector (p, n-p, p*), the last > being the number of non-aliased coefficients. ?summary.glm says: > df a 3-vector of the rank of the model and the number of residual > degrees of freedom, plus number of non-aliased coefficients. It seems to me that the description is reversed: p is

Buenos días, Gracias Carlos, siguiendo el ejemplo que comentas, esto es lo que he introducido en el Scrip de RStudio: *library(xlsx)* *library(xlsxjars)* *library(rJava)* *library(aod)* *R<-read.csv("2002.CSV", sep=";", dec=",", header=T)* *attach(R)* *group<-gl(2,670,1340,labels= c("AVE", "Log.Imports.Value.in.1000.USD"))*

I'm running an ANOVA on some data for respiration in a forest. I am having a problem with my degrees of freedom. For one of my variables I get one fewer degrees of freedom than I should. I have 12 plots and I therefore expected 11 degrees of freedom, but instead I got 10. Any ideas? I have some code and output below: > class(Combined.Plot) [1] "character" >

Dear Everyone in list: I am writing some codes to automate the process of fitting linear models where the names of variables of models are produced and stored in character vectors. But I have problems to pass the vectors to the lm( ) because I don't know how to strip the quotation marks automatically. Here are the codes of the example of lm( ): ## Annette Dobson (1990) "An Introduction

<r-help@r-project.org> Hi, R users, How can I export an R object as a .txt file? As an example I have the result from a regression and I need to save this object in a .txt file ctl <- c(4.17,5.58,5.18,6.11,4.50,4.61,5.17,4.53,5.33,5.14) trt <- c(4.81,4.17,4.41,3.59,5.87,3.83,6.03,4.89,4.32,4.69) group <- gl(2,10,20, labels=c("Ctl","Trt")) weight <- c(ctl,

Hello, I probably have a syntax error in trying to generate an expected survival curve from a weighted cox model, but I can't see it. I used the help sample code to generate a weighted model, with the addition of a "weights=albumin" argument (I only chose albumin because it had no missing values, not because of any real relevance). Below are my code with the resulting error

Hi list, is there any direct way to obtain confidence intervals for the regression slope from lm, predict.lm or the like? (If not, is there any reason? This is also missing in some other statistics softwares, and I thought this would be quite a standard application.) I know that it's easy to implement but it's for explanation to people who faint if they have to do their own programming...

Hola a todos y todas, Gracias por vuestro apoyo en cantidad de preguntas anteriores, de nuevo os escribo para compartir una duda: Estoy trabajando con un modelo bien sencillo, es una regresión simple, pero me gustaría comprobar la significación estadística de cada uno de los coeficientes de regresión en el modelo. La idea es hacer un contraste de hipótesis. Me he descargado el paquete

Hi Everyone, I just don't know how to extract the information I want from the summary of a linear regression model fitting. For example, I fit the following simple linear regression model: results = lm(y_var ~ x_var) summary(results) gives me: Call: lm(formula = y_var ~ x_var) Residuals: Min 1Q Median 3Q Max -5.9859 -1.5849 0.4574 2.0163 4.6015 Coefficients:

Hola: Al tratar de utilizar la función lm R explota, no se si se trata de un problema de mi máquina o me he dado de bruces con un bicho. Después de intentar depurar mi código he probado lo más básico, abrir R (probado en consola y Rkward) y lanzar los ejemplos de lm, esto es lo que me encuentro: ---------- Corte desde la consola ------------------- R version 2.14.0 (2011-10-31) Copyright (C)

Dear all, Hope this question is not too trivial... In order to correct standard errors from an estimation of a fixed effects regression model y need to extract the vector of standard errors of the coefficients of a simple linear model estimated using lm(var1~var2+var3+var4) Unfortunately I can't find out the corresponding function. Thank you very much for your help. Uli -- Uli

Hi, Is there a way to get/extract a matrix of regression variable name, coefficient, and p values? (for lm and glm; which can be sort by p value?) thanks Dhruv [[alternative HTML version deleted]]

Dear all, Is there anybody who can help me to avoid scientific number in the summary of an lm model? Here there is an example of a usual output of the lm model. Thank you! Guillermo Example summary(lm(promiscuity.Index~allK)) Call: lm(formula = promiscuity.Index ~ allK) Residuals: Min 1Q Median 3Q Max -1.67094 -0.13126 0.06703 0.19913 0.40673 Coefficients:

Hello, I am interested in plotting my regression analysis results(regression coefficients and standard errors obtained through OLS and Tobit models) in the form of graphs.What is the best way to accomplish this? Thanks. Murali Kuchibhotla [[alternative HTML version deleted]]

I was hoping that there might be some way to attach significance code like the ones from summary.lm to a dataframe. Anyone know how to do something like that. Here is the function i'd like to add that functionality to: add1.coef <- function(model,scope,test="F",p.value=1,order.by.p=FALSE) { num <- length(model$coefficients) add <- add1(model,scope,test=test) sub <-

Dear all, I would like to ask how to extract the p-value for the whole model from summary(lm). This didn't help a lot summary.lm summary(lm(speed~dist, cars)) Thanks a lot! [[alternative HTML version deleted]]

hi, is there an option to calculate the 'within' & 'between' group variances for a simple ANOVA (aov) model (2 groups, 1 trait, normally distr.) ? or do I have to calculate them from the Sum Sq ? thanks for your time and greetings, gregor -- Gregor Rolshausen PhD Student; University of Freiburg, Germany e-mail: gregor.rolshausen at biologie.uni-freiburg.de tel. :

I think I have isolated a problem with integration between Sweave and beamer. Could you please see the file: http://www.mayin.org/ajayshah/tmp/bugdemo.Rnw Unfortunately, it uses some of my internal libraries, so you can't run it. When I put it through Sweave, I get: http://www.mayin.org/ajayshah/tmp/bugdemo.tex which is, of course, a generic latex file which you can read and