similar to: loglinear analysis

Is there a way to automate fitting and assessing loglinear models for several nominal variables . . . something akin to step or drop1 or add1 for linear or logistic regression? Thanks. --Chris -- Christopher W. Ryan, MD SUNY Upstate Medical University Clinical Campus at Binghamton 40 Arch Street, Johnson City, NY 13790 cryanatbinghamtondotedu "If you want to build a ship, don't drum

Hi, is it possible to do this operation faster? I am going over 35k data entries and this takes quite some time. for(cnt in 2:length(sdata$date)) { if(sdata$value[cnt] < sdata$value[cnt - 1]) { sdata$ddtd[cnt] <- sdata$ddtd[cnt - 1] + sdata$value[cnt - 1] - sdata$value[cnt] } else sdata$ddtd[cnt] <- 0 } return(sdata) Thank you, Benjamin

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

Hi all: Here's a question about result of loglinear analysis. There're 2 factors:area and nation.The raw data is in the attachment. I fit the saturated model of loglinear with the command: glm_sat<-glm(fre~area*nation, family=poisson, data=data_Analysis) After that,I extract the coefficients: result_sat<-summary(glm_sat) result_coe<-result_sat$coefficients I find that all the

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

Hi R-masters! I have a problem and need your help. I have 9 discrete variables with 2 levels each. In exploratory analisys I generate one matrix with chi-square for tables with 2 ariables each with this script setwd("F:/") dados<-read.csv("log.csv")[,2:10] dados.x<-matrix(NA,ncol=9,nrow=9) for(i in 1:8){ for(j in (i+1):9){ tab<-table(dados[,i],dados[,j])

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

google for "centos error performing checksum" refers to createrepo but of course "yum install createrepo" also fails and the instructions appear to be overkill. Is there a "simple" soultion? here's the result of yum update: # yum -v update Loading "fastestmirror" plugin Config time: 0.118 Yum Version: 3.2.22 Setting up Package Sacks Loading mirror

Hi all, I have some some difficulties to use plot and its options. This is what i would like to do. I want to put three plot in a same window so i use par(mfrow=c(1,3)) and i obtain three columns for my three graphics. But, i would like that the two first take 90% of the window and the last 10%. Could anyone give me a soultion?? Best regards, Olivier --

Hi! I had some memory problem with R - hope somebody could tell me a solution. I work with very large datasets, but R cannot allocate enough memoty to handle these datasets. I want work a matrix with row= 100 000 000 and column=10 A know this is 1 milliard cases, but i thought R could handle it (other commercial software like spss could do), but R wrote out everytime: not enough memory.. any

Hi all, I''am new to Prototype and after a 2 hours of research I didn''t find a solution for the following question: I have the following list: <ul id="navi"> <li><a href="/about/">Home</a></li> <li><a href="/about/" class="dropdown" onmouseover="showSubnavi(this)">About</a>

Hello. I am trying to store the output from a loop into a matrix. Failed so far, thanks for the help. What I want to store into a new matrix: (just run the code once): temps<-rnorm(400,14,0.05) ttind<-NULL for(ti in 1:(length(temps)-9)) { if(temps[ti]-temps[ti+9] >= 0.1 && max(temps[ti]-temps[ti+1:9]) > -0.05) ttind<-c(ttind,ti) } for(ti in 1:length(ttind)) {

Hi, I'm not sure why this would be sent directly to me so I'm forwarding it to the CentOS Documentation mailing list as that would seem a more appropriate place. http://wiki.centos.org/Contribute#head-42b3d8e26400a106851a61aebe5c2cca54dd79e5 shake chen wrote: > hello > > I am a team member of iRedMail. http://code.google.com/p/iredmail/w/list > > the iredmail is open

Hi, I have seen many posts regarding this but have not found the soultion. So please cna someone shed some light on this Thanks - padma I get this error when I do try to create a domU guest (xen 3.0.0, 2.6.12.6). The vif is set to vif = [ '' mac=aa:00:00:00:11, bridge=xenbr0'' ] and a "brctl show" shows xenbr0 8000,fefffffff interfaces vif0.0 peth0 Xend.log

Hi everybody, As the Asterisk community is getting larger and larger, I was wondering that the features which are provided in Asterisk and are programmed by the open source community under GPL, or GUIs like FreePBX which also come loaded with wonderful features and uses same Asterisk, are they anywhere violating any patent laws? Most of the features work the same way as Nortel, Avaya and other

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

Hi All, I'm trying to implement a quasi symmetry model for data on twin pairs. A crosstabulation of twin 1 by twin 2 (assumed symmetrical) stratified by another variable. There is a good paper on this by Phil (?) McCloud and Darroch in Biometrika (1995) which explains the method, but I've not done this before so am not clear how to code these models. Any help would be greatly

I have a dataset with 25 subjects and 25 items. For each subject-item combination, there's a 0/1 score for two parts, A and B. I'm thinking of this as a set of 2 x 2 tables, 25 x 25 of them. I'd like to fit a log-linear model to this data to test the independence of the A and B scores. If I ignore the subject and item parts, the following works just fine: glm(count ~ A * B,

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

I'm actively working on Windows with LLVM/clang (not directly contributing yet). First of all, it does compile on Windows (only tested in Visual Studio 2010). But, there are still many problems: (1) x86-64 won't compile as far as I know. CMake generates 64-bit soultion, but one project is missing. But, 32-bit solution compiles safely. (2) Although safe compiliation, there are many