[R] bootstrap

2007-02-28 Thread Indermaur Lukas
Hi, I would like to evaluate the frequency of the variables within the best selected model by AIC among a set of 12 competing models (I fit them with GLM) with a bootstrap procedure to get unbiased results. So I would ike to do the ranking of the 12-model-set 10'000 times separately and

Re: [R] fitting of all possible models

2007-02-27 Thread Indermaur Lukas
coefficients (order of CV=SE/coefficient). Any hint for me? Cheers Lukas Indermaur Lukas wrote: Hi, Fitting all possible models (GLM) with 10 predictors will result in loads of (2^10 - 1) models. I want to do that in order to get the importance of variables (having an unbalanced variable

[R] fitting of all possible models

2007-02-26 Thread Indermaur Lukas
Hi, Fitting all possible models (GLM) with 10 predictors will result in loads of (2^10 - 1) models. I want to do that in order to get the importance of variables (having an unbalanced variable design) by summing the up the AIC-weights of models including the same variable, for every variable

[R] testing slopes

2007-02-20 Thread Indermaur Lukas
Hello Instead of testing against 0 i would like to test regression slopes against -1. Any idea if there's an R script (package?) available. Thanks for any hint. Cheers Lukas °°° Lukas Indermaur, PhD student eawag / Swiss Federal Institute of Aquatic Science and Technology ECO -

[R] dataframe operation

2007-01-24 Thread Indermaur Lukas
hi i have a dataframe a which looks like: column1, column2, column3 10,12, 0 NA, 0,1 12,NA,50 i want to replace all values in column1 to column3 which do not contain NA with values of vector b (100,200,300). any idea i can do it? i appreciate any hint regards lukas °°° Lukas

[R] extract standard errors, write them to file

2007-01-12 Thread Indermaur Lukas
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

[R] batch job GLM calculations

2007-01-11 Thread Indermaur Lukas
Hello I want to batch job the calculation of many GLM-models, extract some values and store them in a file. Almost everything in the script below works (read file, extract values and write them to file) except I fail in indexing the GLM with the modelstructure it should run. Running GLM's