similar to: test for a single variance

Hi R Gurus! This is definitely off topic, but I thought I'd try: what is the way to put in url's into a Latex file, please? I know that you start with a \newcommand statement, with some kind of argument that will change, but I'm not sure how to set it up, please. thanks, Edna Bell

Dear R Gurus: Is there a problem with Rmpi on x86 with SUSE 10.1, please? I've tried everything and it still won't load. Has anyone else dealt with this please? Thanks, Edna Bell

Dear R Gurus: I read somewhere that functions are considered vectors. Is this true, please? thanks Edna Bell

Hi R Gurus! I have a package that I am putting together. I waited until the end to do the documentation and would like to generate a set of "empty" .Rd files. Is there a quick way to do this, please? Thanks, Edna Bell PS I can't seem to access the "Writing R Extensions" manual from 2.6.1 on Windows.

Dear R gurus I would like to take a monthly time series and convert it to a data frame without losing the tsp items, pleae I've tried as.data.frame and data.frame but I get the series without the time element. Any suggestions, please? tia Edna Bell

Hi! Suppose I have a factor: > fac1 <- factor(rep(c("dog","cat","tree"),c(2,3,4))) > fac1 [1] dog dog cat cat cat tree tree tree tree Levels: cat dog tree > length(fac1) [1] 9 Now I want to get rid of the dogs: > fac1 <- fac1[3:9] > fac1 [1] cat cat cat tree tree tree tree Levels: cat dog tree > How should I remove the dog from the

Hello R Gurus: I would like to take a character string and split at the $ sign. I thought that strsplit would do it, but here are the results: > vv [1] "whine$ts1" > vv [1] "whine$ts1" > strsplit(vv,"$") [[1]] [1] "whine$ts1" Does anyone have any suggestions, please? Thanks, Edna Bell

Dear R Gurus: I'm working slowly through "R Programming for Bioinformatics", which is really interesting! Anyway, my question now is: what determines if a function is a builtin vs. a closure, please? For instance: > typeof(sqrt) [1] "builtin" > typeof(mean) [1] "closure" > Thanks, Edna Bell

Hi R Gurus! I have 2 questions, please: a. I'm putting together a little baby package on SUSE Linux 10.1. I want to create a zip of the package for windows. How should I create that please? b. I don't know if the following is possible but here it is: I have R and a bunch of libraries set up very nicely. I would like to hand off this setup to someone else. Should I just copy

Dear R Gurus: I'm trying to find R on another Linux system. I'm using the find command, (surprise), but I only want to see the output where it exists, not all of the other stuff. Is there an option that I could select, please? Right now I have find / -name R TIA, Sincerely, Edna

Hi R Gurus! I'm trying to create a test package using the package.skeleton function. I wanted to add some compiled code too. In the src library, I put together a baby subroutine, compiled it and created a test.dll When I use the R cmd build, it works fine. But I get into trouble with the R CMD check section. /home/Desktop/R-2.5.1/bin # ./R CMD check mypkg * checking for working latex ...

Hi! What is the difference between using data(somedataset) vs. attach(somedataset), please? thanks, Edna Bell

I have put together a test class and methods > setClass("foo",representation(x="numeric")) [1] "foo" > setMethod("plot","foo",function(x,y,...)boxplot(x,...)) [1] "plot" > x <- rnorm(100) > class(x) <- "foo" > plot(x) > y <- new("foo",x=rnorm(50)) > class(y) [1] "foo" My

Hi again! Suppose I have the following: > xy <- round(rexp(20),1) > xy [1] 0.1 3.4 1.6 0.4 1.0 1.4 0.2 0.3 1.6 0.2 0.0 0.1 0.1 1.0 2.0 0.9 2.5 0.1 1.5 0.4 > table(xy) xy 0 0.1 0.2 0.3 0.4 0.9 1 1.4 1.5 1.6 2 2.5 3.4 1 4 2 1 2 1 2 1 1 2 1 1 1 > Is there a way to set things up to have 0 - 0.4 0.5 - 0.9 etc. please? I know there is the cut

Dear R Gurus: How do I find the functions which are primitives, please? Thanks, Edna Bell

Hi R Gurus: Here is some code that I was experimenting with, please: > f1 <- function(x) { + e1 <- new.env(parent=.GlobalEnv) + environment(e1) + print(environment()) + return(mean(x)) + } > f1(1:15) <environment: 0x02525444> [1] 8 > My question: why isn't the environment within the function set to e1, please? Thanks, Edna Bell

Hi still yet again! I have the following code: > try(log(rnorm(25)),silent=TRUE) [1] -0.26396185 NaN NaN -0.13078069 -2.44997193 -2.15603971 NaN 0.94917495 0.07244544 NaN [11] -1.06341127 -0.42293099 -0.53769569 0.95134763 0.93403340 NaN -0.10502078 NaN 0.30283262 NaN [21] -0.11696872 -3.84122332 NaN NaN -0.12808690

Hi R Gurus: I have a tiny data set: > xs [1] 0.7 2.8 0.1 1.9 0.0 1.4 0.2 2.3 0.3 0.2 > and generate the density function for it. I would like to replicate this. Is there a straightforward way to do this, please? According to the help file, the FFT is used. How does this compare to Silverman's version from his page 15 Equation 2.2a, please? Thanks, Edna Bell

What is 1. par.ests <- optimfit$par 2. fisher <- hessb(negloglik, par.ests, maxvalue=maxima); 3. varcov <- solve(fisher); 4. par.ses <- sqrt(diag(varcov)); Thanks a lot, fit.GEV <- function(maxima) { sigma0 <- sqrt((6. * var(maxima))/pi) mu0 <- mean(maxima) - 0.57722 * sigma0 xi0 <- 0.1 theta <- c(xi0, mu0, sigma0) #10/5/2007: removed assign() for maxima.nl

Dear R-List members, I want to check if a set of measurements follows better a gamma or a lognormal distribution (see data below). Using shapiro.test I can test for normality (shapiro.test(log (Lt)). To test for gamma (and normal) distribution I would use ks.test but I need to specify its shape and scale. How should I calculate these values in R? I tried > Lt.fit <- glm(Lt ~ 1,