similar to: print.glm() signif digits (PR#7765)

Hi -- I am having the following problem with table() when applied to vectors of type (mode) integer. When I use the table() command, I can *only obtain an entry in the table for NA values by using exclude=NULL*. Just issuing exclude=NaN will not do it. See below, where x is double at first, and then coerced to integer and notice the difference. Is this a bug or is there something that I do not

After some futile searches, I decided to ask the list to see if any of the sages out there would have an answer: I have a function I wrote a few years ago in S, which calls glim numerous times. I'd like to port it to R, but glm works differently from glim, which takes as part of its input an X design matrix. I probably could write a function to convert glim to glm, but hope this

Dear Experts, In glim() of S Plus, one of the returned values is "var", the estimated variance matrix of coefficients. However, in glm() of R (there is no glim() in R), "var" is not one of the returned values. Anyone know what could I get the varience matrix of coefficients in glm() in R? As a novice in R and S+, I'd appreciate your help Sincerely, Charlie Liu

Dear All, I wonder if anyone can advise me as to whether there is a consensus as to how the effect size should be calculated from GLIM models in R for any specified significant main effect or interaction. In investigating the causes of variation in infection in wild animals, we have fitted 4-way GLIM models in R with negative binomial errors. These are then simplified using the STEP procedure,

Hi all I try to make a weibull survival analysis on R. I know make this on GLIM, and now I try to make the GLIM exercice GLEX8 on R to learning and compare the test. The variables are: time censor group bodymass In GLIM I make: $calc %s=1 $ to fit weibull rather than exponential $input %pcl weibull $ $macro model group*bodymass $endmac$ $use weibull t w %s $ Then, GLIM estimate an alpha for the

Dear fellow R users, I would like to draw a "sunflowerplot" because I have data (decade by month) that plots multiple times on the same x-y co-ordinates. Further I would like to colour each of the points/sunflower leaves on the plot according to the group they belong to (i.e. which type of event each represents within that decade and month). I thought that this would be relatively

I noticed that boxplot computes a 95% CI for the median by using median +/- 1.58*IQR./sqrt(n) Where does the 1.58 constant come from? -- Rick White Statistical Consultant U.B.C. Vancouver B.C. Canada rick at stat.ubc.ca -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send

Hi all, In my continuous transition of GLIM to R I try to make a glm with binomial errors. The data file have 3 vectors: h -> the factor that is ajusted (have 3 levels) d -> number of animais alive (the response) n -> total number of animals To test proportion of alive, make d/n. In GLIM: $yvar d$ $error binomial n$ $fit +h$ scale deviance = 25.730 (change = -9.138) at cycle 4

I came to R from GLIM and its glm. My data sets (ecological community data) are severely over-dispersed, and so I was delighted to find out that the MASS library has glm.nb which is an advancement from the GLIM macros I had used (N.E.Breslow, Applied Statistics 33, 38--44; 1984). However, I need to use offset, but that failed. I am not (yet --- hopefully) fluent enough in R to be able to

Hello, Has anyone written R functions for applying self-controlled case series methods (http://statistics.open.ac.uk/sccs/). In fact only thing needed is to modify glm function to allow absorption of effect. Eg. in Poisson model individual effect is used as factor, but it is considered as nuisance term where parameter estimates are not needed. Could anyone point how absorbing individual

Asterisk as a phone system makes perfect sense in a condo. You can get all the DID's you want and eliminate costs for the owners. You can offer standard FXO for people who don't care and IP sets for people who want to "upgrade" to feature sets. Your door openner is a piece of cake. 1. Create an option in your dialplan only in the "from-access-door" context that reads

Hi all I have another GLIM question. I have been using R as well as Genstat (version 6) in order to fit GLIM models to the data (displayed below). The same models are fitted but the answers supplied by the two packages are not the same. Why? Can anyone help? A discription of the data and the type of model/s fitted can be found below. Regards Allan The

I've downloaded R for windows (9.0.1) and it is great! I've converted all my lecture notes for my GLM course to run on R (they are available on my web page below). I must admit I particularly like the default contrast options, which are identical to GLIM. Also I like the gl function - very useful! I have a couple of questions/bugs: 1. predict.glm doesn't work, but predict.lm does -

I want to do a logistic regression without an intercept term. This option is absent from glm, though present in some of the inner functions glm uses. I gather glm is the standard way to do logistic regression in R. Hoping it would be passed in, I said > r <- glm(brain.cancer~epilepsy+other.cancer, c3, > family=binomial(link="logit"), intercept=FALSE) which produced

Hello, I would like to know if there is any package that allow us to fit Generalized Linear Models via Maximum Likelihood and Linear Models using Generalized Least Squarse in R as the functions glim and gls, respectively, from S-Plus. Also, anybody know if there is any package that fit Log-Linear Models using Generalized Least Squares? Any help will be very useful. Thanks, -- Frederico

I appeal to those entrusted with the keeping of the R flame, the S flame, and the S-Plus flame to relate a bit of history. How did S, S-Plus, and R get their names? Going from S to S-Plus appears clear, a commercial company purchased rights to S, developed a product that they wanted to indicate was related to S, but was more fully developed. This leaves me with only a rumor how S got its original

Dear Friends, My problem is related to how to measure probabilities from a probit model by changing one independent variable keeping the others constant. A simple toy example is like this Range for my variables is defined as follows y=0 or 1, x1 = -10 to 10, x2=-40 to 100, x3 = -5 to 5 Model output <- glim(y ~ x1+x2+x3 -1, family=binomial(link="probit")) outcoef <-

Hello All: I hope I can get someone interested in this problem: Agresti in "Analysis of Categorical Data," p. 289, applies a "row and column effects model" to analyze a two-dimensional cross-classification of ordinal data. He got his results in either SAS or GLIM. Is there a way to replicate his results with R? He claims the RC model fits well with G^2(RC) = 3.57 with df =

I have fitted the faults.data to glm.nb and to the function negbin from the package aod. The output of both is the following: summary(glm.nb(n~ll, data=faults)) Call: glm.nb(formula = n ~ ll, data = faults, init.theta = 8.667407437, link = log) Deviance Residuals: Min 1Q Median 3Q Max -2.0470 -0.7815 -0.1723 0.4275 2.0896 Coefficients:

Hello R-helpers, I was wondering if anyone is porting S lib packages for design of experiments: conf.design, glim, graff(all from statlib)? I was able to change the format of S packages into R (including docs). In conf.design, by W. Venables, sort.list function is missing from R. I just put sort.list <- sort, but I got more error messages. Before doing more work, would like to know if this