similar to: tests performed by anova

Here's a more general question following up on the specific question I asked earlier: Can anybody recommend an R command other than mle.aic() (from the wle package) that will give back a ranked list of submodels? It seems like a pretty basic piece of functionality, but the closest I've been able to find is stepAIC(), which as far as I can tell only gives back the best submodel, not a

Folks, I'm using stepAIC(MASS) to do some automated, exploratory, model selection for binomial and Poisson glm models in R 1.3. Because I wanted to experiment with the small-sample correction AICc, I dug around in the code for the functions glm.fit stepAIC dropterm.glm addterm.glm extractAIC.glm and came across something I just don't understand. stepAIC() passes dropterm.glm() a

Dear R-users, I would like to announce the release of the new package JMbayes available from CRAN (http://CRAN.R-project.org/package=JMbayes). This package fits shared parameter models for the joint modeling of normal longitudinal responses and event times under a Bayesian approach using JAGS, WinBUGS or OpenBUGS. The package has a single model-fitting function called jointModelBayes(),

Dear R-users, I would like to announce the release of the new package JMbayes available from CRAN (http://CRAN.R-project.org/package=JMbayes). This package fits shared parameter models for the joint modeling of normal longitudinal responses and event times under a Bayesian approach using JAGS, WinBUGS or OpenBUGS. The package has a single model-fitting function called jointModelBayes(),

CentOS Errata and Bugfix Advisory 2014:1270 Upstream details at : https://rhn.redhat.com/errata/RHBA-2014-1270.html The following updated files have been uploaded and are currently syncing to the mirrors: ( sha256sum Filename ) x86_64: 6637f912283e727dbb61803200b26ad1974e4ed1c254d8cb315d79be3cb052f6 redhat-lsb-4.1-27.el7.i686.rpm

Note that I sent this to r-devel, yesterday. However, it didn't appear on the mailing list. So, I'm resending it. Today, I plotted the following: > filled.contour (,,z, color.palette=terrain.colors) It looked OK, in R. However, when I created a PDF document, the plot (and other similar plots) had grid (and other) lines in it, that shouldn't be there. Note that this problem is more

All, I'm using lmer for some repeated measures data and have specified the contrasts for a time factor such that say time 3 is the base. This works fine. However, when I next use the subset argument to remove the last two time values, the output indicates that the specified contrast is not maintained (see below). I can solve this by creating a new dataframe for the subset of interest

I know this a newbie question, but I've only just started using AIC for model comparison and after a bunch of different keyword searches I've failed to find a page laying out what the differences are between the AIC scores assigned by AIC() and mle.aic() using default settings. I started by using mle.aic() to find the best submodels, but then I wanted to also be able to make comparisons

Has anyone compiled Vim 7.3 on CentOS 6.2 64bit ? I got the following error after *configure --enable-perlinterp && make*: ./vim.h:2153:21: error: EXTERN.h: No such file or directory ./vim.h:2154:19: error: perl.h: No such file or directory ./vim.h:2155:19: error: XSUB.h: No such file or directory I've installed perl and perl-devel, but with no luck. Any ideas ? Thanks.

Hello, I've noticed that all contrast functions, like contr.treatment, contr.poly, etc., take a logical argument called 'contrasts'. The default is TRUE, in which case they do their normal thing of returning a n x n-1 matrix whose columns are linearly-independent of the intercept. If contrasts=FALSE, they instead return an n x n matrix with full rank (usually the identity matrix,

Send CentOS-announce mailing list submissions to centos-announce at centos.org To subscribe or unsubscribe via the World Wide Web, visit http://lists.centos.org/mailman/listinfo/centos-announce or, via email, send a message with subject or body 'help' to centos-announce-request at centos.org You can reach the person managing the list at centos-announce-owner at centos.org When

Dear R-experts, Here below a toy example with some error messages, especially at the end of the code (Tuning the parameters). Your help to correct my R code would be highly appreciated. ####################################### #libraries library(lattice) library(ggplot2) library(caret) library(randomForest) ?? #Data

Hi Sacha, On second thought, perhaps this is more the direction that you want ... X2 = cbind(X_train,y_train) colnames(X2)[3] = "y" regr2<-randomForest(y~x1+x2, data=X2,maxnodes=10, ntree=10) regr regr2 #Make prediction predictions= predict(regr, X_test) predictions2= predict(regr2, X_test) HTH, Eric On Tue, May 9, 2023 at 6:40?AM Eric Berger <ericjberger at gmail.com>

Dear R users, I am looking for a simple way to define a contrast in a linear model. I have a data frame with two factors (f1, f2) and a dependent y. x <- data.frame(y =rnorm(30), f1=gl(2, 15), f2=gl(3, 5, 30)) Now I want to specify the following contrast: "f1= 1 or 2 and f2=1" vs. "f1= 1 or 2 and f2=3" The best I can come up with is the following: x$new[x$f1==1 &

I would like to try a likelihood ratio test in place of waldtest. Ideally I'd like to provide two glm models, the second a submodel of the first, in the style of lrt (http://www.pik-potsdam.de/~hrust/tools/farismahelp/lrt.html). [lrt takes farimsa objects] Does anyone know of such a likelihood ratio test? Chris Elsaesser, PhD Principal Scientist, Machine Learning SPADAC Inc. 7921

Hello, unfortunately I have I big problem I can't solve. I have to analyse if a gene is tissue specific. For example for the gene xyz I have following expression values: Heart Liver Brain 8.998497 10.013561 12.277407 9.743556 10.137574 11.033957 For every tissue I have two values from two different experiments. Now I want to test if Heart is significant higher

Dear R People: Hello! I'm putting together another RcmdrPlugin package and need to add a FORTRAN subroutine to speed things up a bit. I've never added compiled code to a package. Anyhow, I put my code into a /src directory. Here are the statements: * using log directory 'c:/R/R-2.11.1/bin/RcmdrPlugin.push.Rcheck' * using R version 2.11.1 (2010-05-31) * using session charset:

Emit a deprecation warning when a deprecated function is used, so users have a way to know that they are using one. --- generator/perl.ml | 12 +++++++++++- 1 file changed, 11 insertions(+), 1 deletion(-) diff --git a/generator/perl.ml b/generator/perl.ml index efb31077c..af19650a8 100644 --- a/generator/perl.ml +++ b/generator/perl.ml @@ -335,7 +335,7 @@ PREINIT: List.iter ( fun {

Hi, I know that is not possible make a stepwise procedure using REML in R, I can use ML for this. For nested design it may be very dangerous due the difference in variance structure, mainly in a splitplot design. ML make significative variables that REML dont make. I read an article that is made a stepwise procedure using GENSTAT. from article: "Terms were dropped from a model in a

Dear R-helpers, I have run the code below which I expected to make an object called dd1, but that object does not exist. So, in summary, my problem is that my function is meant to make an object (dd1), and it does indeed make that object (I know that the last line of the function prints it out) but then, after the function has run, the object has disappeared. It's late on a Friday so I may