similar to: loglm fitting

Last week Giovanni Parrinello posted a message asking why various packages were loaded when he loaded an .Rdata file. Brian Ripley replied saying he thought it was because the saved workspace contained a reference to the namespace of ipred. (Correspondence copied below). This begs the question: how did the reference to the namespace of ipred come to be in the .Rdata file? Brian did say it is

Last week Giovanni Parrinello posted a message asking why various packages were loaded when he loaded an .Rdata file. Brian Ripley replied saying he thought it was because the saved workspace contained a reference to the namespace of ipred. (Correspondence copied below). This begs the question: how did the reference to the namespace of ipred come to be in the .Rdata file? Brian did say it is

Hi list, I would like to normalize a matrix (two actually for comparison) using iterative proportional fitting. Using ipf() would be the easiest way to do this, however I can't get my head around the use of the function. More specifically, the margins settings... for a matrix: mat <- matrix(c(65,4,22,24,6,81,5,8,0,11,85,19,4,7,3,90),4,4) using fit <-

I think SAS Macros has capability to call R, and execute it without it being in the picture anywhere. So you can use SAS Macros in a file called R.sas In this you can create a macro called %Describe that can call R , load Hmisc ,run the describe function Note you will need repeated use of %put in this %describe for the mapping to take place Use %INCLUDE to include that file in all SAS

Hello, Eversince I upgraded Samba to 2.0 version I've gotten so many calls from the users regarding Excel 'save as HTML' problem. What seems to be happening is that the users cannot save their Excel spreadsheet as html file when they were able to with the previous Samba version we had installed (Samba 1.9.18p10). According to MS 'knowledge database', this is a known MS Excel

[This email is either empty or too large to be displayed at this time]

Dear R-help community, I am working with multidimensional contingency tables and I am having trouble getting loglm to run on all dimensions without typing out each dimension. I have generated random data and provided output for the results I want below: d1.c1 <- rnorm(20, .10, .02) d1.c2 <- rnorm(20, .10, .02) d2.c1 <- rnorm(20, .09, .02) d2.c2 <- rnorm(20, .09, .02) d3.c1 <-

The null loglinear model is an intercept-only model for log frequency, log(f) = \mu For a one-way table the test of the null model is the same as the chisq.test. This can be fit using loglin(), but I don't think there is any way to specify this using MASS::loglm > t1<- margin.table(Titanic,1) > t1 Class 1st 2nd 3rd Crew 325 285 706 885 > loglin(t1, NULL) 0

loglm() is a port of an original written for S-PLUS. The fact that it carries only a reference to the data frame is neither intentional nor unintentional, but an unnoticed side-effect. I can see advantages both ways. (I'm not so sure, either, that what you say is standard behaviour for model fitting functions really is so universal.) It would not be too hard to come up with a version that

Bug report: there is a problem with loglm (or some underlying code) when used with updating within a function. For example, this code works as expected, when used outside a function: ------------------------------------------------- library(MASS)

I'd like to fit and explore a collection of hierarchical loglinear models that might range from the independence model, ~ 1 + 2 + 3 + 4 to the saturated model, ~ 1 * 2 * 3 * 4 I can use add1 starting with a baseline model or drop1 starting with the saturated model, but I can't see how to get the model formulas or terms in each model as a *list* that I can work with further. Consider

I have numerous replicated goodness-of-fit experiments (observed compared to expected counts in categories) and these replicates are nested within a factor. The expected counts in each cell are external (from a scientific model being tested). The calculations I need within each level of the nesting factor are a heterogeneity G test, with the total G and the pooled G across replicates. Then I

Dear R-helpers, Im trying to fit a Log-linear model on a dataset with bird counts from 60 sites over 14 years for 12 months each (factors for the models). One of the aims is to predict the missing values in this dataset with model predictions. Ive first tried to work with GLM's, that worked fine except for models with one or more interaction-terms. The GLMs run, run .. for hours. So I

Hi all, simple Q: how do I extract the upper and lower CI for predicted probabilities directly for a glm - I'm sure there's a one line to do it but I can't find it. the predicted values I get with the predict (.. "response") Thanks Gerard ********************************************************************************** The information transmitted is intended only for

Dear All, 1. Does anyone have experience of running R on a server inside a Citrix shell - I'd like to get R onto the server and would be greatful for any tips or direction on the matter. 2. This may seem like a silly question so forgive my ignornace. Most of the data I currently work with is held on a number of Lotus Notes (LN) Databases (well it's called a DB here but it's really a

Dear all, The standard call to ARIMA in the base package such as arima(y,c(5,0,0),include.mean=FALSE) gives a full 5th order lag polynomial model with for example coeffs Coefficients: ar1 ar2 ar3 ar4 ar5 0.4715 0.067 -0.1772 0.0256 -0.2550 s.e. 0.1421 0.158 0.1569 0.1602 0.1469 Is it possible (I doubt it but am

Hi all, I have a multiple queing situation I'd like to simulate to get some idea of the distributions - waiting times and allocations etc. Does R has a package available for this - many years ago there used to be a language called "simscript" for discrete event simulation and I was wondering if R has an equivalent (or hopefully with graphics, something better!). Apologies if there

Hi everyone, newbee query! I've installed R 2.8.0 and tried to run this simple glm - x is no of cars in a given year, y is the number voted in an election that year while n is the population 18+: votes <- data.frame(x = c(0.62,0.77,0.71,0.74,0.77,0.86,1.13,1.44), + y=c(502,542,711,653,771,806,934,1123), n=

Dear all, I've replicated the cheese tasting example on p175 of GLM's by McCullagh and Nelder. This is a 4 treatment (rows) by 9 ordinal response (cols) table. Here's my simple code: #### cheese library(MASS) options(contrasts = c("contr.treatment", "contr.poly")) y = c(0,0, 1, 7, 8,8,19, 8,1, 6,9,12,11, 7,6, 1, 0,0, 1,1, 6, 8,23,7,

Hi all, apologies if this is obvious - but I can't see it and would appreciate some quick help! the matrix mhouse is 26x3 and I'm computing odds ratios. The simple code below "should" compute the odds vector for every pair (325) i.e. 26C2 in cols 1 and 2. On the first i=1 outer loop the inner j loop runs from 2 to 26 ok and then I get the error (Error: subscript out of bounds)