similar to: buglet in lm.influence() -- easy fix but... (PR#278)

Greetings! Recently when I tried to use lm.influence I get the following error: Error in .Fortran("lminfl", model$qr$qr, n, n, k, as.integer(do.coef), : Fortran symbol name "lminfl" not in DLL for package "base" This occurs on both Linux and Windows platforms (details below). Searching the mail lists and other sources indicates that the fortran code for

Hi all I'm getting a NaN returned on using dffits, as explained below. To me, there seems no obvious (or non-obvious reason for that matter) reason why a NaN appears. Before I start digging further, can anyone see why dffits might be failing? Is there a problem with the data? Consider: # Load data dep <-

>>>>> "Erin" == Erin Hodgess <hodgess@uhddx01.dt.uh.edu> writes: Erin> Suppose I create a function, x.e, via the fix command: (ESS users ESS don't ever do this, but still ...) >> x.e <- fix(x.e) Erin> When I type in x.e, the function code appears, along Erin> with the name/number of the environment. Ok, in the above it was

Dear statistitians / R users, I was told to analyse the effects of the four factors using lm or update for aov. The following is the result from lm. As you can see that in the last few lines enclosed coefficients. 1)Does this mean I can write a formula, like yield =0.26 + 0.03*H + 1.48*T + 0.04*L + 0.004*C 2) in the two levels design, is there any difference between I use (-1, 1) to

I am performing a locally weighted regression model using housing data, where I only include observations within a certain distance of the house in question. For cross-validation of the bandwidth I am collecting elements of the "hat matrix" (where y hat=hat matrix *y). I was convinced I could grab the diagonal elements for the hat matrix using lm.influence()$hat. In

Hello. I want to get the proper influence function for the glm coefficients in R. This is supposed to be inv(information)*(y-yhat)*x. So I am wondering what is the exact mathematical formula for the output that the functions: dfbeta() OR lm.influence()$coefficients return for a glm model. I am confused because: 1. Their columns don't sum to zero as influences should. 2. They

Hello, I am looking for an explanation and/or fix for a discrepancy in the behaviour of the R lm.influence() function [ version R 1.5.0 (2002-04-29) ] and the same function in Splus [ Splus version 5.1 release 1, running on SGI IRIX 6.2]. The discrepancy is of concern because I am migrating some Splus scripts to R and need to ensure consistency of results. Specifically, when I fit a glm()

Hi, I am trying to plot an interaction.plot with different color for each level of a factor. It has an erratic behavior. For example, it works for the first interaction.plot below, with the example from the ALDA book, but not with the other plots, from the NPK dataset: # from http://www.ats.ucla.edu/stat/r/examples/alda/ch2.htm tolerance <-

Hi all, Perhaps this in an inappropriate post, but I've found a bug in proj I'd like to track down a bit further before making a formal bug report. The example below shows the problem, if you change the rownames proj fails. The problem seems to be that there is a mismatch in the rownames in the qr objects constructed by aov and the rownames that proj is expecting them to have.

Dear Madame, Dear Sir, I am able to obtain the coefficients from a 'summary' of 'lm', but NOT from a 'summary' of 'aov'. The following example shows my steps. ## Initialize rm(list = ls()) # remove (almost) everything in the working environment utils::data(npk, package="MASS") # get data model <- yield ~ block + N*P*K ## Using lm npk.lm <-

Hi, This might be a minor bug, or it could be that I'm just mis-using the tools (again). If you modify the example for model.tables to introduce an NA, N <- c(0,1,0,1,1,1,0,0,0,1,1,0,1,1,0,0,1,0,1,0,1,1,0,0) P <- c(1,1,0,0,0,1,0,1,1,1,0,0,0,1,0,1,1,0,0,1,0,1,1,0) K <- c(1,0,0,1,0,1,1,0,0,1,0,1,0,1,1,0,0,0,1,1,1,0,1,0) yield <-

I'm baffled. When I run make check after installing from source, I get a Error 2. From my understanding of how these things work, it would appear to be coming from this (as at the end of base-Ex.Rout.fail: > has.VR <- require(MASS, quietly = TRUE) Attaching package 'MASS': The following object(s) are masked from package:base : confint confint.lm nclass.FD nclass.scott

I contacted John Fox about this first, because parts of the file are attributed to him. He says that he didn't write rstandard.glm(), and suggests asking r-devel. As it stands, rstandard.glm() has summary(model)$dispersion outside the sqrt(), while in rstandard.lm(), the sd is already sqrt()ed. This seems to follow stdres() in VR/MASS/R/stdres.R. Of course for the c("poisson",

Dear listmembers, I would like to define a class with a slot that takes either an object of class aov or NULL. I have been reading "S4 Classes in 15 pages more or less" and "Lecture: S4 classes and methods" #First I tried with list and NULL setClass(listOrNULL") setIs("list", "listOrNULL") setIs("NULL", "listOrNULL") #doesn't

Full_Name: Naoki Takebayashi Version: 0.65.0+R-release.diff (Oct 6, 1999) OS: Linux/Alpha Submission from: (NULL) (129.79.224.171) This will fix the "problem 2 (crash in fft)" in Bug ID #277 On Linux/Alpha, make check failed because R could not handle the following example in base-Ex.R ##___ Examples ___: # The Old Faithful geyser data data(faithful) : : ## Missing values: x <-

I have determined why there are many "WARNING: Duplicate symbol:" messages. *** My apologies for the length *** There is a lot of detail - but I hope the detail will help R - and others - setup correct options for shared libraries. *** As I press send, I have not stopped testing (my final trial here might not even work) - but! *** the message is that shared libraries do not need to have

I would appreciate it if someone could explain the results of the example from the aov() help file. The output given below is different from book Venables and Ripley - MASS The results In R1.5.1 (Under windows) is as follows: > N <- c(0,1,0,1,1,1,0,0,0,1,1,0,1,1,0,0,1,0,1,0,1,1,0,0) > P <- c(1,1,0,0,0,1,0,1,1,1,0,0,0,1,0,1,1,0,0,1,0,1,1,0) > K <-

The following R commands were written: >help.search("factorial") >data(npk) >npk >coef(npk.aov) In the output of coef command, please explain me the interpretation of coefficients of block1 to block 6 in this 2^3 confounded factorial experiment. Thanks. [[alternative HTML version deleted]]

Dear Sirs, I am trying to compile R v1.2.1 under Digital Unix 4.0D. The configure-script runs without any problems but during 'make' I receive the following error-message: [......] g77 -mieee -g -O2 -c dtrco.f -o dtrco.o g77 -mieee -g -O2 -c dtrsl.f -o dtrsl.o g77 -mieee -g -O2 -c eigen.f -o eigen.o g77 -mieee -g -O2 -c lminfl.f -o lminfl.o ar cr libappl.a Rsock.o approx.o bakslv.o

Not too hard to fix, the following buglet (new in 1.3.0), but I'm about to leave for one week (of vacation).. Here is an example (for Linux) : lynne{sfs}214> cd /tmp lynne{sfs}215> mkdir blue lynne{sfs}216> cd blue /tmp/blue lynne{sfs}217> echo 'data(package = "mva")'|R-1.3.0 --vanilla R : Copyright 2001, The R Development Core Team Version 1.3.0