similar to: BMA and Poisson regression

Hi I'm trying to read SAS-data on VAX/VMS to Windows R using foreign pakage read.xport and experience some problems. Following lines are used in SAS to create XPORT file LIBNAME a ''xxx; LIBNAME b XPORT ''; PROC COPY IN=a OUT=b; RUN; and I succeed in getting file that looks like correct xport file.However when typing in R following

Hello everybody, In order to predict income for different time points, I fitted a linear model with polynomial effects using BMA (bicreg(...)). It works fine, the results are consistent with what we are looking for. Now, we would like to predict income for a future time point t_next and of course draw the prediction interval around the estimated value for this point t_next. I've found the

Hi, I'm trying to using pspline in bic.surv{BMA}. ############################# library(BMA) library(survival) data(veteran) test.bic.surv<- bic.surv(Surv(time,status) ~ karno+pspline(age,df=3)+diagtime+prior, data = veteran, factor.type = TRUE) summary(test.bic.surv, conditional=FALSE, digits=2) ############################# The results are:

Dear friends, In R, the help of "bic.glm" tells the difference between postmean(the posterior mean of each coefficient from model averaging) and condpostmean(the posterior mean of each coefficient conditional on the variable being included in the model), But it's still unclear about the results explanations, and the artile of Rnews in 2005 on BMA still don't give more detail on

R lets one lock an environment with both an R function, base::lockEnvironment, and a C function, R_LockEnvironment, but, as far as I can tell, no corresponding function to unlock an environment. Is this omission on principle or just something that has not been done yet? I ask because several packages, including the well-used R6 and rlang packages, fiddle with some bits in with SET_ENVFLAGS and

Hi Toni, This error is the result of a test that samba does in the module "server_validate", (in password.c) to check the password server for a bug where NT 4 (some versions) would not correctly set the guest bit. There is currently no smb.conf parameter to allow you to control this behavior, so to change it , you would need to actually hack the password.c module to disable it. Not

To whom it may be concernedDear all: i am a statistics of Humboldt university in berlin. I need R-2.4 for my data analysis, but unfortunately, it does not work. The installation works smoothly, but when I double click on the shortcut or from the start menu, the error warning jumps out. My computer is in german system and so do the error information. The following is the error warning and its

Hello! I have data containing a large number of probabilities (about 60) of nonzero coefficients to predict 10 different independent variables (in 10 different BMA models). i've arranged these probabilities in a matrix like so: (IV1) (IV2) (IV3) ... p(b0) p(b0) p(b0) p(b1) p(b1) p(b1) p(b2) p(b2) p(b2) ... where p(b1) for independent variable 1 is p(b1 !=

Hello Chris,  I 've found two other issues  with MAE and CRPS, giving warning  when running examples. I've the same issue on my data.  Hope that you could find some time to take a look here. Thank you Anna > library(ensembleBMA) Loading required package: chron > example(MAE) MAE>   data(ensBMAtest) MAE>   ensMemNames <-

HI, List I want to know weather any body has used BMA Package for model averaging for prediction of future values.....What are the paprmeters we have to suplly to the model. Any script for the same. thanks in advance ANIL KUMAR( METEOROLOGIST) LRF SECTION NATIONAL CLIMATE CENTER ADGM(RESEARCH) INDIA METEOROLOGICAL DEPARTMENT SHIVIJI NAGAR PUNE-411005 INDIA MOBILE +919422023277

I am getting an odd situation when backing up a number of ext3 filesystems and was wondering if it could be caused by journalling. Over the space of a minute the load average will jump from 2 to over 40 and the system will be unresponsive for anywhere from 8 to 25 minutes. I am going to be trying a number of things, but was wondering if anyone could see the reason for the high load given the

I think we can mark this bug closed and check in the patch. Anyone have any concerns/questions before I do? Roy -------- Original Message -------- Subject: [ wxruby-Bugs-6144 ] error: cast from ?void*? to ?int? loses precision Date: Sun, 15 Oct 2006 08:37:31 -0400 (EDT) From: <noreply at rubyforge.org> To: noreply at rubyforge.org Bugs item #6144, was opened at 2006-10-15 00:04

I am trying to fit the 4-parameter logistic model with SSfpl. The estimation is terminated with a message "step factor 0.000488281 reduced below 'minFactor' of 0.000976562". I tried to change the minFactor with nls.control(minFactor=yy) with no apparent success. A call of nls.control() without arguments always shows the value of 0.000976562. Also the estimation terminates as

Hello, I have a little problem which appeared when I installed first 2.2 generation sambaserver on linux and same effect when updated 2.0.3 to 2.2.2 on hp-ux. The actual problem is that my security servers log (NT4 server) get's failure message every time someone connects to samba share which is quite log and it get's frustrating to search real information from logs when this failure

hi, you cannot use the root-principal unless it has administrative privileges, this works for me: kdestroy kinit administrator@YOURREALM net ads join bye, guenther On Thu, Nov 07, 2002 at 08:48:44PM +0200, Jaakko Niemi wrote: > > Starting program: /usr/bin/net ads join > root password: > net: /home/liiwi/cvs-foo/openldap-2.1.8/libraries/libldap/getvalues.c:97:

On 2019-12-27 04:34, Duncan Murdoch wrote: > On 26/12/2019 11:14 p.m., Spencer Graves wrote: >> Hello, All: >> >> >> ? ????? The default "simulate" method for lm and glm seems to ignore the >> sampling variance of the parameter estimates;? see the trivial lm and >> glm examples below.? Both these examples estimate a mean with formula = >>

The article on BMA (Bayesian model averaging) presents most valuable tools for model selection, but I find one detail confusing in Example 1. In page 4 of RNews 5/2, second paragraph says that the probability of Time variable not being in the model is 0.445. It seems to me that the figure should be 1 - 0.445 = 0.555, because p!=0.445 is the prob. of Time variable being in the model. The plot in

Hello R Group, Below is the code I submit: library("BMA") X <- read.table("C:/Documents and Settings/Administrator/Desktop/coding.txt",header=TRUE) Y <- read.table("C:/Documents and Settings/Administrator/Desktop/1DCS.txt",header=TRUE) IGout<- iBMA.glm(X, Y, glm.family= gaussian(), verbose = TRUE, thresProbne0 = 5 ) summary(IGout) IGout I

This is maybe not directly an R problem but I have used R to try to solve it so I think somebody may be able to help. I have a mixture model with three components and a quadratic Scheffe polynomial p1x1+p2x2+p3x3+p12x1x2+p13x1x3+p23x2x3 fitted to the response. Now I'd like to compute the mixture corresponding the maximum response. Model for Y1 has the parameters p1=124.02 p2=60.973 p3=41.479

Hi, I used the bic.surv function, S-PLUS functions developed by Chris Volinsky http://www.research.att.com/~volinsky/bma.html, without problems with S-PLUS. I have to use it with R but I am face with a problem: this function call a fortran routine named "leaps" (answer <- .Fortran("leaps", arguments)). I loaded the leaps library, and the leaps function work well with my R,