similar to: Adding rows to a table with a loop

Dear List, I'm using glmulti for model averaging in R. There are ~10 variables in my model, making exhaustive screening impractical - I therefore need to use the genetic algorithm (GA) (call: method = "g"). I need to include random effects so I'm using glmulti as a wrapper for lme4. Methods for doing this are available here

Dear all, I had a R segmentation fault, and then invoked debug mode and ran step by step. When I reached "terms(Y~X1*X2*...*X16)", I would then have "segmentation" fault. However, if I just ran this under regular "R interactive" mode, it would be fine though taking long time. My questions are: 1. Is there a known limit of terms for a formula? 2. Why does the

Dear all, I had a R segmentation fault, and then invoked debug mode and ran step by step. When I reached "terms(Y~X1*X2*...*X16)", I would then have "segmentation" fault. However, if I just ran this under regular "R interactive" mode, it would be fine though taking long time. My questions are: 1. Is there a known limit of terms for a formula? 2. Why does the

Dear R, The glmulti package no longer loads through the library() command, apparently because of a problem with rJava. I have today reinstalled R from scratch (updated to v2.14.0) and reinstalled all packages from scratch and updated them all too. The problem is the same as I found on v2.13.2. See session below for the error. I tried install.packages(rJava) as advised by the error report but it

Dear List, I am using both the clm() and clmm() functions from the R package 'ordinal'. I am fitting an ordinal dependent variable with 5 categories to 9 continuous predictors, all of which have been normalised (mean subtracted then divided by standard deviation), using a probit link function. From this global model I am generating a confidence set of 200 models using clm() and the

Hello, I'm using the "glmulti" package to run models of all the possible combinations of my variables. However, I am only interested in a few interactions between my variables. I have tried the equivalent of: mod1<-lm(y~a+b+c+a:b) glmulti(mod1, level=1) mod2<-lm(y~a+b+c+a:b) glmulti(mod2, level=2) and glmulti("y", c("a", "b", "c"),

Dear all, I have a question about the glmulti package. I want to include some variables in all models. To that end I applied the wrapper function as shown in the examples (http://www.inside-r.org/packages/cran/glmulti/docs/glmulti). To include the variable "Geslacht" in all models: > glm.redefined = function(formula, data, always="", ...)

I'm running r 2. on a mac running 10.6.4 and a dual-core macbook pro. I'm having a funny time with multicore. When I run it with 2 cores, mclapply, R borks with the following error. The process has forked and you cannot use this CoreFoundation functionality safely. You MUST exec(). Break on __THE_PROCESS_HAS_FORKED_AND_YOU_CANNOT_USE_THIS_COREFOUNDATION_FUNCTIONALITY___YOU_MUST_EXEC__()

All, I am working on a multiple-regression meta-analysis and have too many alternative models to fit by hand. I am using the "metafor" package in R, which generates AIC scores among other metrics. I'm using a simple formula to define these models. For example, rma(Effect_size,variance, mods=~Myco_type + N.type +total, method="ML")->mod where Effect_size is the

hii guys: i get the message from the asterisk: Started music on hold, class 'default', on Local/s at skype-web-callback-dial-263to263-1775,1 [2008-11-11 14:32:41] WARNING[1781]: format_wav.c:156 check_header: Unexpected freqency 11025 [2008-11-11 14:32:41] WARNING[1781]: file.c:322 fn_wrapper: Unable to open format wav [2008-11-11 14:32:41] WARNING[1781]:

i call to my tollfree number buy my CLI send the next error: Aug 14 08:15:22 WARNING[25931]: format_wav.c:169 check_header: Unexpected freqency 22050 Aug 14 08:15:22 WARNING[25931]: file.c:441 ast_filehelper: Unable to open file on /var/lib/asterisk/sounds/procall3.wav Aug 14 08:15:22 WARNING[25931]: file.c:828 ast_streamfile: Unable to open procall3 (format ulaw): No such file or directory Aug

Hi, all I have a medium sized data, 6 years. Each observation is a case with a date variable, such as '2004-08-02'. Some of the months didn't occur a case. I want to plot the 6 years data by month, and the Y_axis is the freqency of cases for each month, meaning 12*6=72 bars or points in the figure. I though of a method, 1st, using the months function, then ploting. But I need to

Hi All, Many Thanks to Underwood for her excellent review of our big trouble which prevent LMS-based AEC algorithms to be used in most computer. Maybe it can be summaried as follows: 1. Different sample rate of sampling and rendering does exists in most low-cost soundcards (In my experiments over more than 20 soundcards, the differences range from 0.5Hz to more than 50Hz when sample rate is set

Dear all R users : I am a IT student with few statistical background and new R user for only have two month exprience. I have a data named medcost, import by read.table() as follow for example (real dataset has 500 cases), the heander id means case id, member means members in a family and cost is the family pay for medical cost every 6 months. id member cost 1 4

Hi, I am trying to fi a glm-poisson model to 400.000 records. I have tried biglm and glmulti but i have problems... can it really be the case that 400.000 are too many records??? I am thinking of using random samples of my dataset..... Many thanks, -- View this message in context: http://r.789695.n4.nabble.com/glm-poisson-fitting-400-000-records-tp3925100p3925100.html Sent from the R help

Hi, I'd like to extend rsync's --log-format to be more detailed and ultimately work in --dry-run mode. While i need both of these enhancements, i would rather not maintain such mods over time for our local use, so before i work on this, i'd like to hear back that this is indeed desirable/desired, as well as any advice, suggestion or feedback you may have in the hope to see this

A few general comments about stepwiseAIC and a suggestion of how to select models a) Apart form the problem, that stepwise selection is not a garanty to get the best model, you need to have a lot of data to avoid overfitting if your model includes 7 parameter plus interactions (> 10 observations per parameter is what you are ideally looking for). b) Have a look at Anderson and Burnham's

Dear All, wonder if you have some thoughts on running the with() function (and perhaps including the pool() function to get the results?) in glmulti? In other words, how to run glmulti with a data set that is produced by mice()? publicly available code: data <- airquality data[4:10,3] <- rep(NA,7) data[1:5,4] <- NA data <- data[-c(5,6)] library(mice) library(glmulti) the following

hello all, I have a simple linear model with 4/5 variables that I am trying to fit. I would like to find the lowest AIC value with any combination of all the variables. I would like to implement this with a while/for loop. Possibly I would like to generalize this so then I can use it when I have many more variables. I do not want to use step AIC. At the moment I am doing it manually but I

Hey folks, I stumbled upon this while looking for something else http://www.linux.com/archive/feature/151340 And it is something I could actually really make use of. But right on that site they list 3 different ones, and so I'm wondering what all is out there and what I should use. Is there one that is part of the standard CentOS? thanks, -Alan -- ?Don't eat anything you've