search for: bialowieza

Full_Name: Kamil Barto? Version: 2.8.0 OS: windows xp Submission from: (NULL) (212.33.92.187) According to the "deparse" function documentation "nlines" is the *maximum* number of lines to produce. But, when "nlines" argument is supplied, it produces exactly nlines of result, and the result contains empty elements at the end. Example: >

Recently, I have been working with some data that look like two overlapping gaussian distributions. I would like to either 1) determine the mean and SD for each of the two distributions OR 2) get some (bayesian ?) statistic that estimates how likely an observation is to belong to the left-hand or right-hand distribution In case I'm using the wrong language, my data looks something like

Hi, I have data for one dependent variable and multiple independent variables y = b0 + b1*x1 + b2*x2 + ... I want to a list of all models that have some subset of the independents (just x1 x2, and not x3, etc.) and their corresponding BIC values. Is there a pre-existing function that does this? I saw that you can calculate individual BIC values using 'lm' and something like AIC(lm1, k

...s that you're not using the whole vector. From your example: Delta <- c(m1 = 0, m2 = 1.8, m3 = 4.2, m4 = 6.2) exp(-0.5*Delta)/sum(exp(-0.5*Delta)) m1 m2 m3 m4 0.63529363 0.25829111 0.07779579 0.02861947 In general the dRedging package at http://www.zbs.bialowieza.pl/users/kamil/r/ can do these problems (I hate to recommend this package because it offers the danger of thoughtless convenience, but if you really know that you want to enumerate models and do IC-based model averaging it can save a lot of time). At the moment, though, it doesn't work with gl...

I have spent a fair amount of time looking for a package that is automated to run glm (binomial) regression models with all possible subsets of my independent variables. Something akin to Lumley's "leaps" package, but can be applied to glms, not just lms; or something similar to Stata's brute force "tryem" function? If anyone can point me in the right direction I

...s that you're not using the whole vector. From your example: Delta <- c(m1 = 0, m2 = 1.8, m3 = 4.2, m4 = 6.2) exp(-0.5*Delta)/sum(exp(-0.5*Delta)) m1 m2 m3 m4 0.63529363 0.25829111 0.07779579 0.02861947 In general the dRedging package at http://www.zbs.bialowieza.pl/users/kamil/r/ can do these problems (I hate to recommend this package because it offers the danger of thoughtless convenience, but if you really know that you want to enumerate models and do IC-based model averaging it can save a lot of time). At the moment, though, it doesn't work with gl...