[R] glm: quasi models with logit link function and binary data
# Hello R Users, # # I would like to fit a glm model with quasi family and # logistical link function, but this does not seam to work # with binary data. # # Please don't suggest to use the quasibinomial family. This # works out, but when applied to the true data, the # variance function does not seams to be # appropriate. # # I couldn't see in the # theory why this does not work. # Is this a bug, or are there theoretical reasons? # One problem might be, that logit(0)=-Inf and logit(1)=Inf. # But I can't see how this disturbes the calculation of quasi-Likelihood. # # Thank you very much, # best, # # Björn set.seed(0) y - sample(c(0,1), size=100, replace=T) # the following models work: glm(y ~ 1) glm(y ~ 1, family=binomial(link=logit)) glm(y ~ 1, family=quasibinomial(link=logit)) # the next model doesn't work: glm(y ~ 1, family=quasi(link=logit)) __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Re: [R] read.spss Error reading system-file header
Hi Jake, I had and still have got the same problem. There were two columns in the data frame which contain just missing values. When I deleted those the impord seemed to work, but there is still a warning message. Warning message: C:\test.sav: Unrecognized record type 7, subtype 13 encountered in system file. Jake Wegelin wrote: Is there any documentation on what kind of SPSS file can and cannot be read by read.spss? Alternatively, how can one modify or clean an SPSS file to make it readable by read.spss? What properties must a *.sav file before read.spss can read it? The file in this example is 270KB, with 5 rows and 173 columns. I have no trouble reading larger files with read.spss, so it's not merely a size problem. I have no difficulty opening the file in SPSS. I also have no trouble getting read.spss to read a dummy SPSS file with 5 rows and 173 columns, where each entry was randomly sampled from c(letters, LETTERS). / library(foreign) / / junk-read.spss(indata/z2EXvideo.sav) / Error in read.spss(indata/z2EXvideo.sav) : Error reading system-file header. In addition: Warning message: indata/z2EXvideo.sav: File layout code has unexpected value 50331648. Value should be 2, in big-endian or little-endian format. Thanks for any information Jake Wegelin / version / __ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html