similar to: fitting the null loglinear model with MASS::loglm??

In running the following example of a loglinear model for the Titanic data, I was surprised to see NaN reported for the Pearson chisq > loglin(Titanic, margin=list(1:3, 4)) 2 iterations: deviation 2.273737e-13 $lrt [1] 671.9622 $pearson [1] NaN $df [1] 15 $margin $margin[[1]] [1] "Class" "Sex" "Age" $margin[[2]] [1] "Survived" Tracing it back,

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

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 <-

Hi all, I've been running loglinear models for three-way tables: one of the variables having three levels, and the other two having two levels each. An example looks like below: > yes.no <- c("Yes","No") > switch <- c("On","Off") > att <- c("BB","AA","CC") > L <- gl(2,1,12,yes.no) > T <-

Kai, Apologies for the double message, it didn't go to the list last time. On Thu, Oct 22, 2015 at 9:59 AM, Kai Nitschke < kai.nitschke at uniklinik-freiburg.de> wrote: > > ?.Call? calls the C/C++ function ?C_LogLin?. But when I am running it line > by line I get > the following error on line 23/24: > Error: object 'C_LogLin' not found > > Hence, my

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)

Hi everyone, I have a question regarding a C function of the "stats" package in R. I tried to understand the ?loglin? basic function of the ?stats? package implemented in R. The implemented function itself runs without any problem (perhaps see sample). When I ran it line by line it stopped at the lines 23-24 of the loglin-function; (the following line): z <- .Call(C_LogLin,

Summary: loglin will fail if there is a global variable called c. I have fixed this for 0.64.2 and 0.65. On Sat, 22 May 1999, Cor en Aylin wrote: > Dear Brian, > > Using the mass library for R (the most recent port dated 16 may 1999) I Could you give version numbers in such reports please: at the moment there is an official version (5.3pl037) of MASS for R, plus a pre-release of MASS

I am using the loglin function of the base package to fit log-linear models. I am interested in obtaining the parameter values and their standard errors. Parameters are easily obtained, but I haven't found the way of obtaining their standad errors. Is this possible with the loglin function? If not, is there any other function to get them? Many thanks, -- Vicente Piorno Departamento de Ecolox?a

Hello, I am almost new in R, so perhaps my question will be silly. I try to use R for analyzing capture-recapture data in epidemiology. A cancer registry has different sources of patients. We know in each list, patients already known in all other list. The aim is to use capture-recapture models for estimating the number of patients unknow of all the sources. Because no order in sources, one

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

Hello, I'm planning to start using R. Before getting into it, I'd like to ask a couple of questions. Does R carry out loglinear model analysis? That is, will it provide the chi-squared goodness of fit test statistic for a given hierarchical loglinear model? Maybe even do a model selection procedure (like Brown's two-step procedure, or forward/backward selection)? Thanks

Dear all, sorry to bother you all with this but I've been trying to use the loglm in MASS package (v2.8.0) and cannot get any sensible output. I'm wondering am I doing something very foolish or missing something obvious. For example, I tried the documentation help(loglm) example - here's the code # Case 1: frequencies specified as an array. sapply(minn38,

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

Hi, I'm analyzing experimental results where two different events ("T1" and "T2") can occur or not during an experiment. I made my experiments with one factor ("Substrate") with two levels ("Sand" and "Clay"). I would like to know wether or not "Substrate" affects the occurrence probability of the two events. Moreover, for each

Dear netters, Sorry for cross-posting this question. I am sure R-Help is not a research methods discussion list, but we have many statisticians in the list and I would like to hear from them. Any function/package in R would be able to deal with the problem from this researcher? ---------- Forwarded message ---------- From: Heidi Bertels Date: Tue, Jun 5, 2012 at 4:31 PM Subject: How to best

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

I am fairly new to log-linear modelling, so as opposed to trying to fit modells, I am still trying to figure out how it actually works - hence I am looking at the interpretation of parameters. Now it seems most people skip this part and go directly to measuring model fit, so I am finding very few references to actual parameters, and am of course clear on the fact that their choice is irelevant for

I'm pleased to announce the release of the vcdExtra package, v. 0.5-0 from R-Forge to CRAN, on its way to a CRAN server near you. vcdExtra was originally designed to serve as a sandbox for introducing extensions of mosaic plots and other visualizations for categorical data, particularly those that apply to (poisson surrogate) loglinear models fitted using glm() and related, generalized

I have been tinkering with mosaicplot() and friends as a way of learning R. As part of this, I've written a pair.table() method for mosaic matrices, and would like to extend mosaicplot to work with loglin and logln (MASS) objects. I'm using R 1.4.0 on Win 98. I've been trying to figure out the formula interface, and think there's a bug, but not sure how to find it, yet alone fix