similar to: Predict in the package R2BayesX

Hi all, I have a problem with the library R2BayesX, when i try to use the command bayesx i get this error: dyld: Library not loaded: /usr/local/lib/libreadline.5.2.dylib Referenced from: /Library/Frameworks/R.framework/Versions/2.15/Resources/library/BayesXsrc/libs/i386/BayesX Reason: image not found I obtain this message also with the example in the bayesx help: ## generate some data

Hi, I wonder if anyone can help me with this issue. I am using R2BayesX. It seems that the model can maximally contain 20 interactions. When the number of interaction terms exceed 20, the code stops working. Here is a piece of toy code. rm(list=ls()) library(BayesX) library(R2BayesX) #data generating model f2<-function(x1,x2,x3,x4) { y<-2*sin(pi*x1)*1.5+exp(2*x2)/3+2 * sin(4 * pi * (x3

Hi all, I have a set of 54 files that I need to convert from ASCII grid format to .shp files to .bnd files for BayesX. I have the following R code to operate on those files: library(maptools) library(Grid2Polygons) library(BayesX) library(BayesXsrc) library(R2BayesX) readfunct <- function(x) { u <- readAsciiGrid(x) } modfilesmore <- paste0("MaxFloodDepth_", 1:54,

Dear All, I have two questions regarding the use of the R2BayesX package for Bayesian analysis. First, is it possible to generate predictions based on the fitted model? According to Gelman and Hill (2007, pp. 361-363), there are at least two ways to do this in BUGS: (1) generate additional data points with the dependent variable coded as missing (and all the independent variables fixed at

Hi, I am trying to find out the scale efficiency and optimal scale of banks by stochastic frontier analysis given the panel data of bank. I am free to choose any model of stochastic frontier analysis. The only approach I know to work with R is to estimate a translog production function by sfa or other related function in frontier package, and then use the Ray 1998 formula to find the scale

I am trying to use the CART in a survival analysis. I have three variables of interest (all 3 ordinal - x, y and z, each of them with 5 categories) from which I want to make smaller groups (just an example 1st category from X variable with the 2nd and 3rd categories from the Y category and 2, 3 and 4 categories from the Z category etc) based on their, let's say, association with mortality. Now

Dear R fans ( and wiki fans), I am just writing a draft to introduce confidence intervals of various "effect sizes" to my students. Surely, I'll recommend the package MBESS in R. Currently, it means I have to recommend R's interface at first. As a statistics teacher in a dept of psychology, I often have to reply why not to teach SPSS. Psychologists and their students hate to

If you mean this paper by Fahrmeir: http://biomet.oxfordjournals.org/cgi/content/abstract/81/2/317 I would recommend BayesX: http://www.stat.uni-muenchen.de/~bayesx/. BayesX interfaces with R and estimates discrete (and continuous) time survival data with penalized regression methods. If you are looking for a bona fide Bayesian survival analysis method and do not wish to spend a lot of time

Dear Mr. Simon Wood, dear list members, I am trying to fit a similar model with gam from mgcv compared to what I did with BayesX, and have discovered the relatively new possibility of incorporating user-defined matrices for quadratic penalties on parametric terms using the "paraPen" argument. This was really a very good idea! However, I would like to constraint the coefficients

Hi Everyone In the Contributed Documentation part of the R Project website there are dozens of various documents explaining this and that on R. Furthermore there is also the document "Introduction to R". In my thesis I have been using R here and there, so I would classify myself as an intermediate user after about 3 years of using it, but I am in no sense a professional. I am now on a

Hi, I am using GAMs (package mgcv) to smooth event rates in a penalized regression setting and I was wondering if/how one can select the order of the derivative penalty. For my particular problem the order of the penalty (parameter "m" inside the "s" terms of the formula argument) appears to have a larger effect on the AIC/deviance of the estimated model than the

According to this: http://www.redhat.com/support/wpapers/redhat/ext3/why.html ext3 is forward and backward compatible with ext2... Any user who wishes to un-journal a file system can do so easily... I am assuming un-journalling is the equivalent of converting it to ext2. How do you do this? I haven't been able to find anything. The reason I want to do this is so I can modify my

Hi All R Experts, I met with below mentioned statistics in paper "Stock Index Volatility Forecasting with High Frequency Data" by Eugenie Hol, Siem Jan Koopman http://ideas.repec.org/p/dgr/uvatin/20020068.html I would like to ask that what is "Box-Ljung portmantacau statistic based on N squared autocorrelation" ? Is it same as "Box-Ljung Statistics" of stats

Hi I found on google that there is an implementation of hot-deck imputation in SAS: http://ideas.repec.org/c/boc/bocode/s366901.html Is there anything similar in R? Many Thanks Eleni Rapsomaniki

Hello R users & R friends, I just want to ask you if density() can produce a confidence interval, indicating how "certain" the density() line follows the true frequency distribution based on the sample you feed into density(). I've heard of loess.predict(loess(y ~ x), se=TRUE) which gives you a SE estimate of the smoothed scatterplot - but density() kernel smoothing is not the

I only recently discovered options("available_packages_filters" = list(add = TRUE, "license/FOSS")) [cf. help("available.packages", package="utils") in R 2.10.0 or later] which goes nicely with my options("checkPackageLicense" = TRUE) [new in R 2.11]. But now I want to purge my library of packages I would not have installed had I known about

Hi I am trying to install "rms" package but while installing it shows following error package 'survival' 2.36-2 is loaded, but >= 2.36.3 is required by 'rms' what to do? i am using linux OS I have tried by updated r-base-core but it didnt work regards GRR [[alternative HTML version deleted]]

Hi Everybody, please can you help me how to call BayesX in R in order to see the graph already exist in BayesX Thanks ---------- Forwarded message ---------- From: <r-help-request@r-project.org> Date: Sat, Jan 22, 2011 at 5:00 AM Subject: R-help Digest, Vol 95, Issue 22 To: r-help@r-project.org Send R-help mailing list submissions to r-help@r-project.org To subscribe or

CRAN (and crantastic) updates this week New packages ------------ * Bergm (1.0) Alberto Caimo http://crantastic.org/packages/Bergm Functions implementing Bayesian estimation for exponential random graph models via exchange algorithm Updated packages ---------------- lmtest (0.9-26), logcondens (1.3.5), MTSKNN (0.0-4), pmml (1.2.21), r2lUniv (0.9.4), rattle (2.5.11), rgdal (0.6-23),

Hello Folks, We have 2 (identical hardware) busy public 2U rack servers, both have several hundred users & peak email (for example) deliveries at 300/5min during the daytime. Both servers run RAID1 over 2 EIDE disk drives. One server is fine - it never has problems - it is running RH7.3 and CPU loads are normal. The other server has constant disk activity & kjournald at CPU load of