similar to: Mixing text and graphics in PDF

I was using R 1.5.1 today with an old script from S-Plus and discovered that R does not have reorder.factor. It's a very simple function, but very handy for re-arranging factors in lattice graphics (among other places). Would it make sense to add this function to the base code in order to increase compatability with S-Plus and add additional functionality? reorder.factor <-

I just installed R-1.4.0 on my Win2k machine from SetupR.exe. Everything seems to be fine with the base and recommended packages. However, when I install the lattice and grid packages they do not seem to work. Everything I try results in something like this: Instruction 0x00543f06 refernced memory 0x00000000. The memory cannot be "read". and I am getting kicked out of R. Has

Is there a function that provides the Fisher information matrix for a generalized linear model? I do not see how to access the off-diagonal matrix elements of the value returned by glm. (I'm particularly interested in logistic regression.) If not, what is a good way to use R to compute Hessians or other partial derivatives of log likelihoods? I would appreciate any guidance. David

Does anyone know how to use the equivalent of debug() on an S4 method? I would like R to enter the browser not for the generic function, but for the method of the class that I specify. Thanks, David _____________________________ David Bickel http://davidbickel.com Research Scientist Pioneer Hi-Bred International Bioinformatics & Exploratory Research 7250 NW 62nd Ave., PO Box 552 Johnston,

Given a numeric vector of observations, does R have any generic way to estimate the parameters of commonly used distributions (exponential, gamma, etc.) without numerically optimizing the likelihood function? Thanks, David _______________________________________ David R. Bickel http://davidbickel.com Research Scientist Pioneer Hi-Bred International Bioinformatics & Exploratory Research 7250

I encountered an unexpected difference in the way that mean and var work: > mean(numeric(0)) [1] NaN > var(numeric(0)) Error in var(numeric(0)) : `x' is empty Looking into this, I discovered in mean.default that sum(x)/n evaluates as > sum(numeric(0))/length(numeric(0)) [1] NaN The function var calls an .Internal and seems to evaluate this: >

Is there a way I can get R to automatically enter the browser inside a user-defined function on the generation of an error? Specifically, I'm trying to debug this: Error in as.double.default(sapply(lis, FUN)) : (list) object cannot be coerced to double In addition: There were 38 warnings (use warnings() to see them) > traceback() 8: as.double.default(sapply(lis, FUN)) 7:

I'm using R 1.6.1 on Windows 95 with cygwin bash as my shell. The 'shell' command tries to determine which shell the user prefers via the following steps: shell <- Sys.getenv("R_SHELL") if (!nchar(shell)) shell <- Sys.getenv("SHELL") if (!nchar(shell)) shell <- Sys.getenv("COMSPEC") For me, shell

Is there any R function that computes the convolution of the double exponential distribution? If not, is there a good way to integrate ((q+x)^n)*exp(-2x) over x from 0 to Inf for any value of q and for any positive integer n? I need to perform the integration within a function with q and n as arguments. The function integrate() is giving me this message: "evaluation of function gave a

Hi List, Am trying to customize a legend in trellis: Draws 2x5 lines in 5 colors and 2 linetypes. I would like to add two more items to the legend showing the key for the line types above the colored legend. Any suggestions welcome - thanks Herry ############################# #Following example code: library(gregmisc) trellis.device(bg="white") i1=0 i2=-1.89767506 i3=-1.17087085

Hey Timo, I think I've asked you about this before, and I'm sure its not a big priority, but this is really a pain when trying to check default settings... The output of dovecot -a is huge, and the fact that it is not sorted alphabetically (nor is dovecot -n, but thats not as big of a problem, although I'd like to see that sorted as well for consistency) makes it really difficult to

Dear R-users, When I tried to import a text file (tab delimited) which has 2000+ rows with the following command (With the importData in S, it works though), x <- read.table(textfile, sep= "\t", skip=5, stringAsFactors=F) I received the following warning message: Error in scan(file, what, nmax, sep, dec, quote, skip, nlines, na.strings,: number of items read is not a multiple of

Does anyone know how to interrupt R in RedHat? Neither control-c nor Esc is working. What I don't want to do is close the window or kill the entire R process. Thanks, David This communication is for use by the intended recipient and ...{{dropped}}

After reading the help page on set.seed, I am unsure about how to solve the following problem. I need to call function f a thousand times. The list of values returned by f, should be as random as possible. f calls g twice: f <- function(){g1 <- g(1); g2 <- g(2); c(g1; g2)} The function g in turn calls sample and returns a number, but also depends on its argument, so, starting from the

I am attempting to utilize BackupPC on a Fedora 14 server to backup a remote client. As I understand, it's primary mechanism for finding clients is performing a nmblookup <clientname> This works fine for computers connected to the local network. My issue is extending ?Samba's? search to encompass our other network - the point-to-point VPNs Using OpenVPN, we have a number of road

An NA in the rownames of a matrix (or dataframe) causes misalignment when the matrix is printed: R> x <- matrix(1:12, 3,4, dimnames=list(letters[1:3], LETTERS[1:4])) R> rownames(x)[2] <- NA R> x A B C D a 1 4 7 10 <NA> 2 5 8 11 c 3 6 9 12 The bug is in function Rstrlen, in src/main/printutils.c. MatrixRowLabel and MatrixColumnLabel (same file) rely on Rstrlen

Easy question that I can't find an answer for. I'm trying to subset a data frame and want to exclude the positive values, i.e. I want the NA values. My data: > summary(temp$tuna) Min. 1st Qu. Median Mean 3rd Qu. Max. NA's 1 2 3 3 4 5 1211 Querying for subset(temp, tuna %in% "NA", select.... subset(temp, tuna == NA,

Hi, I'm sorry to send email to everyone on this list, but I have a simple question which is bothering me and I can't seem to figure out the correct answer. I just downloaded R and I'm trying to reproduce some simple analysis I've done on other packages. In particular, if I have a vector of experimental values and a vector of their uncertainties, is there a way to plot the

Hi all, I'm building a plot of the values in tmeant (below) against positions 1 to 5, using matplot. tmeant looks like this: case1 case2 pos1 861.8466 818.5909 pos2 961.2841 976.3466 pos3 878.6080 1262.8523 pos4 950.8011 1129.6080 pos5 968.1080 1063.3920 I also have lower (object tl) and upper (object tu) bounds on the confidence intervals as follows: tl: pos1

Hello all, In the course of trying to implement the CDF of an InverseGammaDistribution, I have run across the need for an igamma() function. Several others have needed this function but the answers I have found so far are not totally clear to me. I'm writing for three reasons: 1) to present a small error in the docs 2) to clarify the approach we are expected to take 3) to request,for the