similar to: R-help post Bayesian CART

Hi All Please see the Rprofile file which i have modified as follows and after that when I start R then I see that R says to me "TRUE" for all the packages implying that all loaded at once. But when i try to use commands as simple as help("lm"), it doesnt work nor any of the menu "Packages" is not working. Although the regression using lm ( Y ~ X ) is working

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

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

Hello everyone, at the moment I'm using the tgp package for modelling a nonstationary data set on a two dimensional area D and I'm interested in the prediction and the estimated covariance matrix. For this purpose I'm using the function btgp. As far as I understand, btgp uses a MCMC algorithm to split up D along lines parallel to the coordinate axes and estimates independent

Hello, I was using the functino kgs of maptree and I've realized maptree has been removed from CRAN repository. Anybody knows why has it been removed? or where can I find the function kgs that computes the Kelley-Gardner-Sutcliffe penalty function for a hierarchical cluster tree Thank you. [[alternative HTML version deleted]]

I am using the latest version of R 2.13.1 and need to load the following packages; library(maps) library(mapdata) library(mapproj) library(lattice) library(tgp) library(spatstat) library(akima) I get the follwing error when loading the packages, "maps", "mapdata" and "mapproj" Error: package 'maps' is not installed for 'arch=i386' Error: package

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:

Hi ! I have been using BMA (bayesian model Averaing) package for modeling purposes, but was faced with a problem of incorporating offset term in the Poisson regression of disease rates. It looks like bic.glm does not accept offset keyword like glm ? Any ways to solve the problem using wt option in BMA? Janne Pitk?niemi -- Department of Public Health P.O.Box 41 (Mannerheimintie 172) 00014

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

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

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 all. Are there any R functions around that do quick logistic regression with a Gaussian prior distribution on the coefficients? I just want posterior mode, not MCMC. (I'm using it as a step within an iterative imputation algorithm.) This isn't hard to do: each step of a glm iteration simply linearizes the derivative of the log-likelihood, and, at this point, essentially no

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 have created a tree and want to save some of the data so that I can create a html table from it. I would like to save the output from data.ltr (see example below) to a file, but haven't found a way to do that, keeping the nice format that typing data.ltr gives me (see output below). Is there a way to do this? Example: library (maptree) library (tree)

Hi All: Recently I am researching my dissertation about the quantile model selection by bayesian approach. I have the dependent variable(return) and 16 independent variables and I need to select the best variable for each quantile of return. And the method I used is the bayesian approach, which is based on calculating the posterior distibution of model identifier. In other words, I need to obtain