similar to: partial association model

If I do not make a mistake,the partial association model is an extension of log-linear model.I read a papers which gives an example of it.(Sloane and Morgan,1996,An Introduction to Categorical Data Analysis,Annual Review of Sociology.22:351-375) Can R fit such partial association model? ps:Another somewhat off-topic question.What's the motivations to use log-linear model?Or why use

Is it generally known/has it been previously discussed here that the $aic() component in GLM-family objects (e.g. results of binomial(), poisson(), etc.) does not as implemented actually return the AIC, but rather -2*log-likelihood + 2*(model_has_scale_parameter) ? Can anyone in this forum gauge how a documentation patch would be received? This behaviour does not seem to be documented in ?family

Hello, I am now running a multiple linear regression program, but I do not know the difference between the command step and stepAIC. Thanks. Maggie [[alternative HTML version deleted]]

Hi I have two processes which take with a certain probability (p1 and p2) x number of years to complete (age1 and age2). As soon as thge first process is completed, the second one begins. I want to calculate the time it takes for the both processes to be completed. I have the following script which gives me the answer, butI think there must be a more elegant way of doing the calculations

At a loss. I've removed and reinstalled all the ports that use openssl.so.10 now, and still having the same issue. At this point, I'm looking at just upgrading the whole system to freebsd 12.0 and seeing what happens. Not sure why the core files don't seem to cooperate with gdb, butI assume that is some freebsd configuration issue. If anyone has any other ideas though I am all ears.

Hi everyone, I tried to solve this problem but I could not find the solution. I have about 105 matrices of equal size in the memory of**R, I need to do is extract from these matrices some known positions and create a new matrix with these columns. Show you an example with only three matrices (but in my case I have hundreds of them). A=matrix(1:15,5) [,1] [,2] [,3] [1,] 1 6 11 [2,] 2 7 12

Hello I want to compile some R-packages under windows (rmutil and revent for instance from Pr Lindsey) I followed the instrictions given in the help. The compilation succeeded, but I can't run the programs. I receive the message: "error in dyn.oad(file): unable to load shared library "C:\R\Rwo650/library/rmutil" LoadLibrary failure. I found the given path quite strange!!

Dear all, I am running a logistic regression and this is the output: glm(formula = educationUniv ~ brncntr, family = binomial) Deviance Residuals: Min 1Q Median 3Q Max # ???? ????? ?? ???????? -0.8825 -0.7684 -0.7684 1.5044 1.6516 Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) -1.06869 0.01155 -92.487 <2e-16 *** brncntrNo

Time, not sure who owner is, but I've tried as both admin and myself, to no avail. ?However, I told you that mesg is y, buti tried to send a mesg using ?' write ' ?and it says I have messages disabled. ????? ? I did that by opening two?terminals, on I logged in as ' su ', and sent the message from there. ?Is that peculiar? Might account for why I am notgetting anything from

Hello, Are there any packages/functions available for testing the assumptions underlying assumptions for a good GLM model? Like linktest in STATA and smilar. If not, could somebody please describe their work process when they check the validity of a logit/probit model? Regards, Jay

Chris, Finally got time to fully play with your suggestions. Had to re-code a bit of stuff but the concept works - just like a bought one :~). Many thanks for the solution, Kind Regards, Eric. For those curious - here''s what it was about; On Monday 02 January 2006 02:10, Eric Sloane tried to type something like: > Hi, > I''m trying to figure the most efficient way to

Hi Listers, I''m wondering about the possibility of creating a tabbed partial (or set of partials) dependant on a parent-children structure. For example; I have a parent JOB - job has many sons (Ooops wrong story). JOB contains CONTACTS, NOTES, STATISTICS, WORKFLOW, etc.. I want to show the general JOB data at the top of my screen, then below that load a set of tabs containing the

How the hell do you close em off after the action is completed? Ta, Eric

I am trying to fit a rank-frequency distribution with 3 unknowns (a, b and k) to a set of data. This is my data set: y <- c(37047647,27083970,23944887,22536157,20133224, 20088720,18774883,18415648,17103717,13580739,12350767, 8682289,7496355,7248810,7022120,6396495,6262477,6005496, 5065887,4594147,2853307,2745322,454572,448397,275136,268771) and this is the fit I'm trying to do: nlsfit

Hello everybody! I'm trying some datamining, but i'm having some problems with arule package, in the end of processing R "had to be closed". I already tryied to reinstall the 2.10 version, change the computer and realocated more virtual memory. Does anyone had this problem to? I had a hiphoteses that i have to prepare the data, somehow i don't know. Thanks for helping!!!

Hello When I run scatter.smooth(jitter(weight), jitter(height2), span = .25, evaluation = 50, pch = '.') I get the type of graph I thought I would get, but also a warning..... k-d tree limited by memory. ncmax= 528 I always get concerned when there are warnings I don't understand. What's a k-d tree? Is this something to be concerned about? Thanks Peter Peter L. Flom,

I see something new and unexpected here. > update.packages() trying URL `http://cran.r-project.org/src/contrib/PACKAGES' Content type `text/plain; charset=iso-8859-1' length 163467 bytes opened URL .......... .......... .......... .......... .......... .......... .......... .......... .......... .......... .......... .......... .......... .......... .......... ......... downloaded

Hi, I have the following code fragment; @comm_type = Communication.find_by_sql(''select id from communication_types where upper(description) = "EMAIL"'') Which results in >> @comm_type => [#<Communication:0x39a0f60 @attributes={"id"=>"4"}>] But when I try to test an ActiveREcord against this it doesn''t work;

If you have an interest in providing rich web interfaces for your applications beyond the abilities of DHTML, including audio and video, you may want to join me and 225 other developers on the Laszlo on Rails project. Our Google group is at http://groups.google.com/group/laszlo-on-rails. Be a part of something that will be featured in a talk at RailsConf -- join us!

This works: > model.matrix(~I(pos>3),data=data.frame(pos=c(1:5))) (Intercept) I(pos > 3)TRUE 1 1 0 2 1 0 3 1 0 4 1 1 5 1 1 attr(,"assign") [1] 0 1 attr(,"contrasts") attr(,"contrasts")$"I(pos > 3)" [1] "contr.treatment"