[R] extract standard errors, write them to file

From: Indermaur Lukas <Lukas.Indermaur_at_eawag.ch>
Date: Fri 12 Jan 2007 - 13:55:26 GMT


Hello

I want to repeatedly extract coefficients and standard errors from a GLM and write them into a file (1row=all coefficients of model A, 2 row=all coefficients of model B, etc.). I can extract coefficients but not standard errors. furthermore I fail to write extracted values line by line into the predifined matrix G.

I appreciate any idea to solve the problem

best regards

Lukas  

#---------start code------------------------------------------------------------------------------------------------------------------------------------

global                           <- formula(logHRS~Ri + E  + Co + LWD +Alwd +W  + T2 + A  + N + Sex + y)      #1
richness_evenness        <- formula(logHRS~Ri + E  + D1 + D2  +D3   +D4 + D5 + D6 + N + Sex + y)        #2

all_models <- c(global, richness_evenness) for (i in 1:length(all_models))
{
ts.model <- glm(all_models[[i]],family=gaussian,data=t.data)

G            <- matrix(NA,length(all_models),length(all_models))
G            <- coefficients(ts.model)             #regression coefficents (betas)
}
write.table(G, paste(t.url, file="Coefficients.txt"), sep="\t", quote=F)

#---------end code-------------------------------------------------------------------------------------------------------------------------------------
     

°°°
Lukas Indermaur, PhD student
eawag / Swiss Federal Institute of Aquatic Science and Technology ECO - Department of Aquatic Ecology
Überlandstrasse 133
CH-8600 Dübendorf
Switzerland  

Phone: +41 (0) 71 220 38 25
Fax : +41 (0) 44 823 53 15
Email: lukas.indermaur@eawag.ch
www.lukasindermaur.ch



R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. Received on Sat Jan 13 00:59:45 2007

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.1.8, at Fri 12 Jan 2007 - 14:30:26 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-help. Please read the posting guide before posting to the list.