Re: [R] Bootstrapping Confidence Intervals for Medians

2007-01-06 Thread Prof Brian Ripley
On Sat, 6 Jan 2007, [EMAIL PROTECTED] wrote: I apologize for this post. I am new to R (two days) and I have tried and tried to calculated confidence intervals for medians. Can someone help me? Later, you say you want a confidence interval for a difference in medians, not the same thing.

[R] negative binomial family glm R and STATA

2007-01-06 Thread Patrick Giraudoux
Dear Lister, I am facing a strange problem fitting a GLM of the negative binomial family. Actually, I tried to estimate theta (the scale parameter) through glm.nb from MASS and could get convergence only relaxing the convergence tolerance to 1e-3. With warning messages:

Re: [R] eval(parse(text vs. get when accessing a function

2007-01-06 Thread jim holtman
I agree with what you are saying. That is the reason I have What is the problem you are trying to solve on my signature. The other way of saying that is Tell me what you want to do, not how you want to do it. On 1/5/07, Thomas Lumley [EMAIL PROTECTED] wrote: On Fri, 5 Jan 2007, Ramon

Re: [R] help with gls

2007-01-06 Thread Markus Jäntti
On Sat, 2007-01-06 at 01:21 -0500, Zhenqiang Lu wrote: Hello R-users, I am using gls function in R to fit a model with certain correlation structure. The medol as: fit.a-gls(y~1,data=test.data,correlation=corAR1(form=~1|aa),method=ML) mu-summary(fit.a)$coefficient With the toy data I

[R] listing all functions in R

2007-01-06 Thread Gavin Simpson
Dear List, I'm building an R syntax highlighting file for GeSHi [*] for a website I am currently putting together. The syntax file needs a list of keywords to highlight. How can I generate a list of all the functions in a base R installation? Ideally the list would be formatted like this:

Re: [R] eval(parse(text vs. get when accessing a function

2007-01-06 Thread Ramon Diaz-Uriarte
On 1/5/07, jim holtman [EMAIL PROTECTED] wrote: The other reason for considering which of the different approaches to use would be performance: f.1 - function(x) x+1 f.2 - function(x) x+2 system.time({ + for (i in 1:10){ + eval(parse(text=paste('f.', i%%2+1,

[R] ANCOVA

2007-01-06 Thread John Cardinale
Hello, I am analyzing a data set that has possible covariates. It is also multivariate. Are there any R function which can do analysis of covariance? Thanks. John Cardinale [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch

Re: [R] listing all functions in R

2007-01-06 Thread Prof Brian Ripley
Could you tell us what you mean by - 'function' (if() and + are functions in R, so do you want those?) - 'a base R installation'? What is 'base R' (standard + recommended packages?) And on what platform: the list is platform-specific? Here is a reasonable shot: findfuns - function(x) {

Re: [R] eval(parse(text vs. get when accessing a function

2007-01-06 Thread Gabor Grothendieck
Here are a few more timing alternatives from slowest to fastest. Of the ones that dynamically construct the name of the frunction (all except last example), get seems to be the fastest although its still 6x slower than the if statement. f.1 - function(x) x+1 f.2 - function(x) x+2 n - 10^4; i

Re: [R] eval(parse(text vs. get when accessing a function

2007-01-06 Thread Ramon Diaz-Uriarte
Dear Greg, On 1/5/07, Greg Snow [EMAIL PROTECTED] wrote: Ramon, I prefer to use the list method for this type of thing, here are a couple of reasons why (maybe you are more organized than me and would never do some of the stupid things that I have, so these don't apply to you, but you can

Re: [R] listing all functions in R

2007-01-06 Thread Gavin Simpson
On Sat, 2007-01-06 at 13:48 +, Prof Brian Ripley wrote: Could you tell us what you mean by Thank you for your reply, Prof. Ripley. - 'function' (if() and + are functions in R, so do you want those?) I was thinking about functions that are used like this: foo() So I don't need things

Re: [R] eval(parse(text vs. get when accessing a function

2007-01-06 Thread Ramon Diaz-Uriarte
On 1/6/07, Thomas Lumley [EMAIL PROTECTED] wrote: On Fri, 5 Jan 2007, Ramon Diaz-Uriarte wrote: I see, this is direct way of dealing with the problem. However, you first need to build the f list, and you might not know about that ahead of time. For instance, if I build a function so

Re: [R] eval(parse(text vs. get when accessing a function

2007-01-06 Thread Gabor Grothendieck
I guess the problems with eval(parse(...)) are: - speed - verbosity - security - length of string to be parsed get is faster, less verbose and safer than eval(parse(...)) in this case. The length of string problem seems not applicable here but I have enountered one situation where it was.

Re: [R] eval(parse(text vs. get when accessing a function

2007-01-06 Thread Brian Ripley
On Sat, 6 Jan 2007, Ramon Diaz-Uriarte wrote: Dear Greg, On 1/5/07, Greg Snow [EMAIL PROTECTED] wrote: Ramon, I prefer to use the list method for this type of thing, here are a couple of reasons why (maybe you are more organized than me and would never do some of the stupid things that

Re: [R] listing all functions in R

2007-01-06 Thread Duncan Murdoch
On 1/6/2007 9:25 AM, Gavin Simpson wrote: On Sat, 2007-01-06 at 13:48 +, Prof Brian Ripley wrote: Could you tell us what you mean by Thank you for your reply, Prof. Ripley. - 'function' (if() and + are functions in R, so do you want those?) I was thinking about functions that are

Re: [R] listing all functions in R

2007-01-06 Thread Gabor Grothendieck
The arguments to the functions can differ too even if they exist on multiple platforms. system() on Windows has the input= argument but not on UNIX. On 1/6/07, Duncan Murdoch [EMAIL PROTECTED] wrote: On 1/6/2007 9:25 AM, Gavin Simpson wrote: On Sat, 2007-01-06 at 13:48 +, Prof Brian

Re: [R] listing all functions in R

2007-01-06 Thread Gavin Simpson
On Sat, 2007-01-06 at 10:43 -0500, Duncan Murdoch wrote: On 1/6/2007 9:25 AM, Gavin Simpson wrote: On Sat, 2007-01-06 at 13:48 +, Prof Brian Ripley wrote: Could you tell us what you mean by Thank you for your reply, Prof. Ripley. - 'function' (if() and + are functions in R, so

Re: [R] listing all functions in R

2007-01-06 Thread Gavin Simpson
On Sat, 2007-01-06 at 10:58 -0500, Gabor Grothendieck wrote: The arguments to the functions can differ too even if they exist on multiple platforms. system() on Windows has the input= argument but not on UNIX. That's a good point Gabor, and one I hadn't considered as yet. As I'm only just

Re: [R] ANCOVA

2007-01-06 Thread Michael Kubovy
On Jan 6, 2007, at 8:34 AM, John Cardinale wrote: Are there any R function which can do analysis of covariance? ?lm RSiteSearch('ancova') _ Professor Michael Kubovy University of Virginia Department of Psychology USPS: P.O.Box 400400Charlottesville, VA

Re: [R] eval(parse(text vs. get when accessing a function

2007-01-06 Thread Ramon Diaz-Uriarte
On 1/6/07, Brian Ripley [EMAIL PROTECTED] wrote: On Sat, 6 Jan 2007, Ramon Diaz-Uriarte wrote: (...) cvFunct - function(whatever, genefiltertype, whateverelse) { internalGeneSelect - eval(parse(text = paste(geneSelect, genefiltertype, sep

Re: [R] ANCOVA

2007-01-06 Thread Ramon Diaz-Uriarte
On 1/6/07, Michael Kubovy [EMAIL PROTECTED] wrote: On Jan 6, 2007, at 8:34 AM, John Cardinale wrote: Are there any R function which can do analysis of covariance? ?lm RSiteSearch('ancova') Given the question, you'll probably need to find how to do an ancova with lm. Several documents in

[R] memory problem

2007-01-06 Thread Zoltan Kmetty
Hi! I had some memory problem with R - hope somebody could tell me a solution. I work with very large datasets, but R cannot allocate enough memoty to handle these datasets. I want work a matrix with row= 100 000 000 and column=10 A know this is 1 milliard cases, but i thought R could handle

Re: [R] ANCOVA

2007-01-06 Thread Richard M. Heiberger
Please look at the help file ?ancova in the HH package. Please get HH_1.17 which was up on CRAN yesterday. The source file has propagated to the mirrors. As of a few minutes ago, the Windows binary is on cran.at.r-project.org but not yet at the mirrors. Be sure to have history recording on in

[R] Using VGAM's vglm function for ordinal logistic regression

2007-01-06 Thread Inman, Brant A. M.D.
R-Experts: I am using the vglm function of the VGAM library to perform proportional odds ordinal logistic regression. The issue that I would like help with concerns the format in which the response variable must be provided for this function to work correctly. Consider the following example:

[R] has anyone implemented LARS with the positive lasso?

2007-01-06 Thread David Reiss
Hi, I am interested in a modification to LARS that allows for positive-only constraints in the variables (with details about how to implement this as described in section 3.4 of the Efron et al (2003) LARS paper). Before I dive into the lars package code myself, I was wondering if anyone knew of