[R] (no subject)

2003-11-27 Thread Len Vir
The results of your email command are provided below. Attached is your original message. - Unprocessed: original message. - Results: Ignoring non-text/plain MIME parts - Unprocessed: I registered to your mailing list, but i don't get any mail from you. I checke

Re: [R] dvips function gives "Documents not found" error

2003-11-27 Thread Dirk Eddelbuettel
On Thu, Nov 27, 2003 at 04:21:04PM -0800, David Kelly wrote: > Frank Harrell wrote: > > > Did you make the changes I suggested (surrounding two items by > > dQuote())? > > > > Have you tried Linux? > > No, I didn't make those changes because I've just been working with a > binary dis

Re: [R] dvips function gives "Documents not found" error

2003-11-27 Thread David Kelly
Frank Harrell wrote: > Did you make the changes I suggested (surrounding two items by > dQuote())? > > Have you tried Linux? No, I didn't make those changes because I've just been working with a binary distribution. I may go ahead and pull down sources and try; I was trying to avoid

Re: [R] FDA and ICH Compliance of R

2003-11-27 Thread Frank E Harrell Jr
On Thu, 27 Nov 2003 09:59:57 -0500 (EST) "Gabor Grothendieck" <[EMAIL PROTECTED]> wrote: > > From: Frank E Harrell Jr <[EMAIL PROTECTED]> > > per year in SAS licenses and have to hire armies of non-intellectually > > challenged SAS programmers to do the work of significantly fewer > > programmers

[R] Re: FDA and ICH Compliance of R

2003-11-27 Thread Paul Sorenson
> "Antonia Drugica" <[EMAIL PROTECTED]> writes: > > > Does anybody know if R is FDA or ICH (or EMEA...) > compliant? AFAIK S-Plus > > is but that means nothing... > > As Thomas pointed out, that does mean nothing -- there was a group of > folks discussing what might be done to help, earlier this

Re: [R] dvips function gives "Documents not found" error

2003-11-27 Thread Frank E Harrell Jr
On Thu, 27 Nov 2003 10:04:51 -0800 David Kelly <[EMAIL PROTECTED]> wrote: > This is primarily an FYI to the Hmisc author, though any brilliant > suggestions from him or anyone else are always welcome. Re this problem > > reported earlier in this thread: > >getting an error from Hmisc when

[R] cclust - cindex - binary data

2003-11-27 Thread Bruno Giordano
Hi, I'm trying to debug a function I wrote to calculate the cindex for a hierarchical tree. For this it is useful to compare my calculations with those in output from the clustindex function, in the cclust library. There's no way, however, to have the cindex value for a given output of the cclust f

[R] Getting rid of loops?

2003-11-27 Thread Xiao-Jun Ma
I wrote a function to calculate cosine distances between rows of a matrix. It uses two loops and is slow. Any suggestions to speed this up? Thanks in advance. theta.dist <- function(x){ res <- matrix(NA, nrow(x), nrow(x)) for (i in 1:nrow(x)){ for(j in 1:nrow(x)){ if (i > j)

Re: [R] tcltk - tkcreate question

2003-11-27 Thread Thomas Stabla
On 27 Nov 2003, Peter Dalgaard wrote: > Thomas Stabla <[EMAIL PROTECTED]> writes: > > > > > I have successfully translated the code until the line > > > >$base.fm.cv create window 0 0 -anchor nw -window $hull -width 10 -height 500 > > -tag window > > > > which i don't fully understand because

Re: [R] tcltk - tkcreate question

2003-11-27 Thread Peter Dalgaard
Thomas Stabla <[EMAIL PROTECTED]> writes: > > I have successfully translated the code until the line > >$base.fm.cv create window 0 0 -anchor nw -window $hull -width 10 -height 500 -tag > window > > which i don't fully understand because i started with tcltk just this > week. > > I tried

[R] tcltk - tkcreate question

2003-11-27 Thread Thomas Stabla
Hello, i'm trying to translate following tcltk source code, which I found in newsgroup comp.lang.tcl, written by Tom Wilkason, into R Code. proc scrolled_Canvas {base} { frame $base.fm -borderwidth 2 -relief sunken canvas $base.fm.cv -yscrollcommand "$base.fm.cv_vertscrollbar set" scrol

Re: [R] lme v. aov?

2003-11-27 Thread Prof Brian Ripley
On Thu, 27 Nov 2003, Spencer Graves wrote: > Do you want to make inference about the specific subjects in your > study? If yes, the subjects are a fixed effect. If instead you want to > make inference about the societal processes that will generate the > subjects you will get in the fut

Re: [R] lme v. aov?

2003-11-27 Thread John Christie
Its not so much that I wasn't getting the difference between fixed and random effects. Although, I do like the way you put the comment below. For my purposes subject is a random effect. It was more on correct notation in lme with repeated measures designs (my a and b are repeated while the m

Re: [R] MASS fitdistr()

2003-11-27 Thread Prof Brian Ripley
On Thu, 27 Nov 2003, Oscar Linares wrote: > Dear R experts, > > I am trying to use the R MASS library fitdistr() to fit the following > list: Well, that is not a list! > k21stsList<-c(0.76697,0.57642,0.75938,0.82616,0.93706,0.77377,0.58923,0.37157,0.60796,1.00070,0.97529,0.62858,0.63504,0.68697

Re: [R] dvips function gives "Documents not found" error

2003-11-27 Thread David Kelly
This is primarily an FYI to the Hmisc author, though any brilliant suggestions from him or anyone else are always welcome. Re this problem reported earlier in this thread: getting an error from Hmisc when trying dvips(latex(describe(mtcars)),file="/kellytest/kelly.ps") that says Error i

[R] MASS fitdistr()

2003-11-27 Thread Oscar Linares
Dear R experts, I am trying to use the R MASS library fitdistr() to fit the following list: k21stsList<-c(0.76697,0.57642,0.75938,0.82616,0.93706,0.77377,0.58923,0.37157,0.60796,1.00070,0.97529,0.62858,0.63504,0.68697,0.61714,0.75227,1.16390,0.66702,0.83578) as follows, library(MASS) fitdistr(k

Re: [R] significance in difference of proportions: What problem are you solving?

2003-11-27 Thread Spencer Graves
Hi, Torsten: Thanks for the reference to library(exactRankTests). That seems like a reasonable alternative to "prop.test" with small samples. However, aren't "exact tests" and the related bootstrap methodology what Deming called "enumerative techniques", more relating to describin

RE: [R] significance in difference of proportions

2003-11-27 Thread Ted Harding
On 27-Nov-03 [EMAIL PROTECTED] wrote: > I've 2 samples A (111 items) and B (10 items) drawn from the same > unknown population. Witihn A I find 9 "positives" and in B 0 > positives. I'd like to know if the 2 samples A and B are different, > ie is there a way to find out whether the number of "posit

Re: [R] significance in difference of proportions

2003-11-27 Thread Torsten Hothorn
> Hello, > > I'm looking for some guidance with the following problem: > > I've 2 samples A (111 items) and B (10 items) drawn from the same unknown > population. Witihn A I find 9 "positives" and in B 0 positives. I'd like to > know if the 2 samples A and B are different, ie is there a way to fin

Re: [R] significance in difference of proportions

2003-11-27 Thread Jonathan Baron
On 11/27/03 17:04, [EMAIL PROTECTED] wrote: >Hello, > >I'm looking for some guidance with the following problem: > >I've 2 samples A (111 items) and B (10 items) drawn from the same unknown >population. Witihn A I find 9 "positives" and in B 0 positives. I'd like to >know if the 2 samples A and B a

Re: [R] lagsarlm - using mixed explanatory variables (spdep package)

2003-11-27 Thread Roger Bivand
On Thu, 27 Nov 2003, Roy Sanderson wrote: > Hello Usually it is easier to send package questions directly to the package maintainer, because they may not interest the whole list. > > I'm very new to R (which is excellent), so apologies if this has already > been raised. In the spdep package,

Re: [R] lme v. aov?

2003-11-27 Thread Spencer Graves
Do you want to make inference about the specific subjects in your study? If yes, the subjects are a fixed effect. If instead you want to make inference about the societal processes that will generate the subjects you will get in the future, that is a random effect. The function "lme" ha

Re: [R] help me!

2003-11-27 Thread Spencer Graves
English may have greater variation between regional dialects than does Spanish, The Scots roll their r's, and I suspect it is more like the Spanish rolls (perro) than taps (pero), though I have not had enough contact with the Scots to judge. The English r may be closer to the r in Portug

[R] lme v. aov?

2003-11-27 Thread John Christie
I am trying to understand better an analysis mean RT in various conditions in a within subjects design with the overall mean RT / subject as one of the factors. LME seems to be the right way to do this. using something like m<- lme(rt~ a *b *subjectRT, random= ~1|subject) and then anova(m,type

RE: [R] R 1.8.1 on SUSE 9.0

2003-11-27 Thread J.R. Lockwood
> > I think the problem is that "by default" g77 is not installed. However you > should still be able to find the rpm on the CDROM. > > HTH, > Andy Thanks to all for your replies. Indeed the package gcc-g77 was on the install disk, and I was able to install the R rpm with no problems. J.R. L

[R] significance in difference of proportions

2003-11-27 Thread Arne.Muller
Hello, I'm looking for some guidance with the following problem: I've 2 samples A (111 items) and B (10 items) drawn from the same unknown population. Witihn A I find 9 "positives" and in B 0 positives. I'd like to know if the 2 samples A and B are different, ie is there a way to find out whether

Re: [R] problems with R graph.

2003-11-27 Thread Prof Brian Ripley
If indeed you used the bitmap() function, please read its help page. The differences are very likely due to the installations of ghostscript and nothing to do with R. There never was an `R 1.6', but the current version of R is 1.8.1, so it looks as if an update is well overdue. On Thu, 27 Nov 20

[R] problems with R graph.

2003-11-27 Thread aurelie . defferrard
Hello, I have some problems to generate graph with R... I am working on two different platform : - Compaq Alpha Server (Unix True 64 5.1) + R 1.6 - Sparc Server (Sun Solaris 8) + R 1.6 I use different functions like the bitmap function, the legend function and the barplot function. The gra

Re: [R] FDA and ICH Compliance of R

2003-11-27 Thread Gabor Grothendieck
From: Frank E Harrell Jr <[EMAIL PROTECTED]> > per year in SAS licenses and have to hire armies of non-intellectually > challenged SAS programmers to do the work of significantly fewer > programmers that use modern statistical computing tools like R and S-Plus, > it is surprising that SAS is still

Re: [R] OT: apt-get and R in Debian

2003-11-27 Thread Stefano Calza
> Can you give us more details, please? Which version of Debian > (stable, testing or unstable) are you running and what are the > relevant parts of your /etc/apt/sources.list file? (Just send us the > whole file if you are not sure what parts are relevant.) Ok. Find attach my source.list + the

[R] help me!

2003-11-27 Thread Patricia
Would you help me to answer this question, please? /r/ is one of the most difficult English sounds to acquire and imitate. Describe in what ways it is different from our Spanish rolls (perro) and taps (pero). How does the pronunciation of this sound vary according to the context? Find at l

Re: [R] FDA and ICH Compliance of R

2003-11-27 Thread Frank E Harrell Jr
On Thu, 27 Nov 2003 07:48:02 +0100 "Antonia Drugica" <[EMAIL PROTECTED]> wrote: > I'm quite new to this medical stuff. But my associates told me that we > are not free in choice of Statistical Software because the FDA has high > standards concerning this topic. But if they would prefer a specific

Re: [R] FDA and ICH Compliance of R

2003-11-27 Thread Thomas Lumley
On Thu, 27 Nov 2003, Antonia Drugica wrote: > I'm quite new to this medical stuff. But my associates told me that we are > not free in choice of Statistical Software because the FDA has high > standards concerning this topic. But if they would prefer a specific > package (like SAS) that could mean

Re: [R] OT: apt-get and R in Debian

2003-11-27 Thread Dirk Eddelbuettel
On Thu, Nov 27, 2003 at 01:11:48PM +0100, Stefano Calza wrote: > Sorry for the OffTopic, but I always have a problem using apt-get to > update my debian siystem and R. > Anytime it updates the packages (right now I installed a self-compiled > version of 1.8.1), even if they are exactly the same.

[R] stl and NA

2003-11-27 Thread "Unternährer Thomas, uth"
Hi, I try to figure out what the stl-function exactly do. I was reading the paper by Cleveland et al. (1990) and tested some features of stl (the ability to decompose time series with missing values and the robustness feature). I tried the following: > data(co2) > co2.na <- co2 > is.na(co2.na[c

Re: [R] would like to know how to simulated a GARCH(1,2)

2003-11-27 Thread Adrian Trapletti
Follow the example in tseries, we can simulated a GARCH(0,2), n <- 1100 a <- c(0.1, 0.5, 0.2) # ARCH(2) coefficients e <- rnorm(n) x <- double(n) x[1:2] <- rnorm(2, sd = sqrt(a[1]/(1.0-a[2]-a[3]))) for(i in 3:n) # Generate ARCH(2) process { x[i] <- e[i]*sqrt(a[1]+a[2]*x[i-1]^2+a[3]*x[i-2]

[R] lagsarlm - using mixed explanatory variables (spdep package)

2003-11-27 Thread Roy Sanderson
Hello I'm very new to R (which is excellent), so apologies if this has already been raised. In the spdep package, I'm trying to undertake an autoregressive mixed model using the lagsarlm function. This is working fine, but there does not appear to be a method of including an explanatory variable

Re: [R] Blocked Mail Notification

2003-11-27 Thread Peter Dalgaard
[EMAIL PROTECTED] writes: > * eManager Notification ** > > Recipient, Content filter has detected a sensitive e-mail. > > Source mailbox: "[EMAIL PROTECTED]" > Destination mailbox(es): "[EMAIL PROTECTED]" > > *** End of message *** > > Re

RE: RE: [R] Correlation test in time series

2003-11-27 Thread Adrian Trapletti
Thanks for you help, And how to test covariance = zero in time series , cov(r_t, r_t-1)=0 and r_t are homoscedastik and dependent ? How about: ?acf ?pacf in package 'ts' Box.test from package 'ts' best Adrian -- Dr. Adrian Trapletti Trapletti Statistical Computing Wildsbergstrasse 31,

[R] OT: apt-get and R in Debian

2003-11-27 Thread Stefano Calza
Hi everybody, Sorry for the OffTopic, but I always have a problem using apt-get to update my debian siystem and R. Anytime it updates the packages (right now I installed a self-compiled version of 1.8.1), even if they are exactly the same. Anybody can help me? TIA, Stefano __

Re: [R] FDA and ICH Compliance of R

2003-11-27 Thread Spencer Graves
Have you tried "www.r-project.org" -> search -> "R site search"? This issue has been discussed in the past. With luck, you may find something there that might help. hope this helps. spencer graves Antonia Drugica wrote: I'm quite new to this medical stuff. But my associates told

Re: [R] would like to know how to simulated a GARCH(1,2)

2003-11-27 Thread Patrick Burns
Prelude for those not in the know: GARCH models the variance of a times series conditional on past information (often only the series itself). It is a reasonably good model of the variance of the returns of market-priced assets, which display big jumps upwards in variance followed by gradual deca

RE: [R] multiple peaks in data frame

2003-11-27 Thread Rieckermann Joerg
Hello! why not try the R site search for "peaks": http://finzi.psych.upenn.edu/search.html which gives you (among other things): http://finzi.psych.upenn.edu/R/Rhelp02a/archive/7593.html Hope this helps, J. Joerg Rieckermann Environmental Engineering Swiss Federal Institute for Environmental

[R] would like to know how to simulated a GARCH(1,2)

2003-11-27 Thread M. M. Palhoto N. Rodrigues
Follow the example in tseries, we can simulated a GARCH(0,2), n <- 1100 a <- c(0.1, 0.5, 0.2) # ARCH(2) coefficients e <- rnorm(n) x <- double(n) x[1:2] <- rnorm(2, sd = sqrt(a[1]/(1.0-a[2]-a[3]))) for(i in 3:n) # Generate ARCH(2) process { x[i] <- e[i]*sqrt(a[1]+a[2]*x[i-1]^2+a[3]*x[i-2]^

[R] Workshop on Statistics in Functional Genomics 2004

2003-11-27 Thread Christina Kuenzli
Apologies in advance if you receive multiple copies of this email. This is to announce and invite your participation in a workshop on Statistics in Functional Genomics, to be held from 27 June - 2 July 2004 at Ascona in the Italian-speaking part of Switzerland. The purpose of the workshop is to br