similar to: New version of "catspec" package

A multinomial logit model can be specified as a conditional logit model after restructuring the data. Doing so gives flexibility in imposing restrictions on the dependent variable. One application is to specify a loglinear model for square tables, e.g. quasi-symmetry or quasi-independence, as a multinomial logit model with covariates. Further details on this technique and examples with several

I've put two packages for R on my home page at http://www.xs4all.nl/~jhckx/R/. The "pcnt" package is for multiway percentage tables. I've posted a first effort called "ctab" on this group and a request for enhancing "ftable" with percentages on the wishlist. The "mcl" package is for estimating multinomial logistic models using conditional logistic

R has amazing capabilities, but percentage tables are a weak spot IMHO. There's prop.table but that's rather unwieldly, especially for multiway tables. CrossTable by Marc Schwartz in the gregmisc library makes percentage tables a breeze but is limited to two-way tables. So I decided to try my own hand at writing an R-function that would make it easy to produce nicely formatted percentage

I've posted a sample file for estimating loglinear models for square tables (mobility models) at http://www.xs4all.nl/~jhckx/mcl/R/ Comments and suggestions are welcome. John Hendrickx

Hi,  I tried fitting loglinear model using the glm(catspec). The data used is FHtab. . An independence model was fitted. Here summary() and fitmacro( ) give different values for AIC.   I understand that fitmacro( ) takes the likelilhood ratio L2(deviance) to calculate AIC and uses the formula AIC= L2- d.f(deviance)*2 and this AIC is used in nested models. (Am I right?)   The value given by

Hi,  I tried fitting loglinear model using the glm(catspec). The data used is FHtab. . An independence model was fitted. Here summary() and fitmacro( ) give different values for AIC.   I understand that fitmacro( ) takes the likelilhood ratio L2(deviance) to calculate AIC and uses the formula AIC= L2- d.f(deviance)*2 and this AIC is used for comparison of nested models. (Am I right?)   The value

Full_Name: John Hendrickx Version: 1.6.0 OS: Windows 98 Submission from: (NULL) (137.224.174.216) CrossTable in the "gregmisc" package fails when the fisher.exact test produces an error (I suspect this is because the number of cases is too large). This can be fixed using "FTt <- try(fisher.test(t, alternative = "two.sided"))" or by making the test optional.

Hi,Dear all R experts, I am trying to do the 2-way contingency table analysis by fitting the loglinear models. However, I found my table has several empty cells which are theoretically missing values.I have no idea of how to solve them coz we cannot compute the simulated p-value with zero marginals.Does someone have some suggestions? Please help me out, thanks a lot! Cheers, Yan

I am attempting to produce crosstabulations between two variables for subgroups defined by a third factor variable. I'm using by() and CrossTable() in package gmodels. I get the printing of the tables first and then a printing of each level of the INDICES. For example: library(gmodels) by(warpbreaks, warpbreaks$tension, function(x){CrossTable(x$wool, x$breaks > 30,

Dear list, I am trying to create a three-way table with percent occurrence instead of raw frequencies. However, I cannot get the results I expected: I have the following table: > ftable(table( mannerDF$agem, mannerDF$target, mannerDF$manner )) <snip> 50 bak 0 0 0 0 1 0 pak 0 0 0 0 3 0 sak

Dear All, I want to create a table for several variables. As example. I have a dataframe with following data: Gender transport driving 1 0 1 0 1 0 1 0 1 Now I want to create a table in the following form: Gender

The HTB qdisc has a compile time option, HTB_HYSTERESIS, that trades accuracy of traffic classification for CPU time. These patches change hysteresis to be a runtime option under the control of "tc". The effects of HYSTERESIS on HTB''s accuracy are significant (see chapter 7, section 7.3.1, pp 69-70 in Jesper Brouer''s thesis: http://www.adsl-optimizer.dk/thesis/ ),

Dear list. New to R, I'm looking for a way of using crosstab to output low-dimensional (higher than 2) contingency tables (frequencies, per-cents by rows, % by columns, mean, quantiles....) I'm looking for something of the following sort dataframe: singers, categorical variates: voice category (soprano,mezzo-soprano, ...) , voice type( drammatic, spinto, lirico-spinto, lirico,

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

Hello, I have a dataset with a continuous independent variable (fish length, range: 30-150 mm) and a binary response (foraging success, 0 or 1). I want to discretize fish length into 5 mm bins and give the proportion of individuals who successfully foraged in each each size bin. I have used the cut function to discretize the length values into my desired bins, but I can't figure out how to

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

Full_Name: John Hendrickx Version: 1.6.2 OS: Windows XP Submission from: (NULL) (80.126.78.108) (This is not a bug report but a request to add a feature to future versions of R. Hope this is an appropriate place). I'd like to suggest adding an option to ftable to allow percentages. It would be easy to do and backwards compatible. Percentage tables are useful in educational contexts, whereas

The HTB qdisc has a compile time option, HTB_HYSTERESIS, that trades accuracy of traffic classification for CPU time. These patches change hysteresis to be a runtime option under the control of "tc". The effects of HYSTERESIS on HTB''s accuracy are significant (see chapter 7, section 7.3.1, pp 69-70 in Jesper Brouer''s thesis: http://www.adsl-optimizer.dk/thesis/ ),

-- apologies if this is a dup - i got a bounce saying the message was unprocessed. Is there a straightforward way to get a table with percents in the cells rather than counts? I've looked at table, ftable, xtabs, and ctab, which did the conversion but returned the results in a single row without labels. any suggestions are appreciated. thank you.

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