Skip to content

loglm fitting

2 messages · Gerard M. Keogh, Sunil Suchindran

#
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, function(x) length(levels(x)))
         ## hs phs fol sex f
         ##  3   4   7   2 0
         minn38a <- array(0, c(3,4,7,2), lapply(minn38[, -5], levels))
         minn38a[data.matrix(minn38[,-5])] <- minn38$fol
         fm <- loglm(~1 + 2 + 3 + 4, minn38a)  # numerals as names.
         deviance(fm)
         > [1] 0

The deviance is zero.
I tried other examples as well (Laura Thompson) and get pretty much
nonsense.

E.g.

      library(MASS)
      options(contrasts=c("contr.treatment", "contr.poly"))
      tbl.8.3 = data.frame(mara
      =factor(c("Yes","No"),levels=c("Yes","No")),
                            cig
      =factor(c("Yes","No"),levels=c("No","Yes")),
                            alc
      =factor(c("Yes","No"),levels=c("Yes","No")),
                            count=c(911,538,44,456,3,43,2,279) )
      fit.acm = loglm( count ~ alc*cig*mara,data=tbl.8.3, param=T, fit=T)


Any advice greatly appreciated.

Gerard


**********************************************************************************
The information transmitted is intended only for the person or entity to which it is addressed and may contain confidential and/or privileged material. Any review, retransmission, dissemination or other use of, or taking of any action in reliance upon, this information by persons or entities other than the intended recipient is prohibited. If you received this in error, please contact the sender and delete the material from any computer.  It is the policy of the Department of Justice, Equality and Law Reform and the Agencies and Offices using its IT services to disallow the sending of offensive material.
Should you consider that the material contained in this message is offensive you should contact the sender immediately and also mailminder[at]justice.ie.

Is le haghaidh an duine n? an eintitis ar a bhfuil s? d?rithe, agus le haghaidh an duine n? an eintitis sin amh?in, a bhearta?tear an fhaisn?is a tarchuireadh agus f?adfaidh s? go bhfuil ?bhar faoi r?n agus/n? faoi phribhl?id inti. Toirmisctear aon athbhreithni?, atarchur n? leathadh a dh?anamh ar an bhfaisn?is seo, aon ?s?id eile a bhaint aisti n? aon ghn?omh a dh?anamh ar a hiontaoibh, ag daoine n? ag eintitis seachas an faighteoir beartaithe. M? fuair t? ? seo tr? dhearmad, t?igh i dteagmh?il leis an seolt?ir, le do thoil, agus scrios an t-?bhar as aon r?omhaire. Is ? beartas na Roinne Dl? agus Cirt, Comhionannais agus Athch?irithe Dl?, agus na nOif?g? agus na nGn?omhaireachta? a ?s?ideann seirbh?s? TF na Roinne, seoladh ?bhair chol?il a dh?chead?.
M?s rud ? go measann t? gur ?bhar col?il at? san ?bhar at? sa teachtaireacht seo is ceart duit dul i dteagmh?il leis an seolt?ir l?ithreach agus le mailminder[ag]justice.ie chomh maith. 
***********************************************************************************