[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

[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]^2) }

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

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

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 me

[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: 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

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 *** Received: from

[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

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]^2) }

[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

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.

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] 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

[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

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

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 the

[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

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

[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] 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. Lockwood

[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] 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

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

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, I'm

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 are

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 find out

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 describing a

[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)

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 in

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!

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

[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

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 to

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 i started with

[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

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 trying

[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 year, but

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] 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 distribution. I