[R] Fortran 95 in R ?

2005-05-18 Thread Joel Bremson
Is it possible to run Fortran 95 code from R? I don't think so, but hopefully someone can prove me wrong. Here is the test I tried: A little fortran 95 subroutine: subroutine allloc() real, dimension(:, :), allocatable :: a integer :: n n = 10 allocate(a(n,n+1)) end I then compiled: g95 -c

SV: [R] Sweave and paths

2005-05-18 Thread =?iso-8859-1?Q?S=F8ren_H=F8jsgaard?=
Don't think it is possible. To overcome the problem I wrote a small R-program which replaced foo with foo\bar in the Sweave file. Clumpsy, sure, but... Best Søren Fra: [EMAIL PROTECTED] på vegne af Bill Rising Sendt: ti 17-05-2005 20:31 Til: rhelp Emne: [R]

[R] dendrogram and dendrapply

2005-05-18 Thread Martin Olivier
Hi all, I think I have some problems to use correctly the function dendrapply. Let suppose dend is a dendrogram object. I would likde to know the cardinal number of leaves depending of each node of the tree dend. It is right that the command attr(dend,members) gives the total number of leaves,

RE: [R] Vuong test

2005-05-18 Thread =?iso-8859-1?Q?St=E9phanie_PAYET?=
Hi Reza, thanks for your answer. I have already used the code of Dr. Jackman for fitting the double hurdle model, but it doesn't allow to specify different regressors for the zero and the poisson (or negative binomial) parts, which is a bit limitative. Moreover, I don't have Stata, I just dispose

[R] Testing for warning inside functions

2005-05-18 Thread Peter Wolf
I am looking for a way to get a warning message immediately after an evaluation within a function. To get error messages you can use geterrmessage(). But I found no function that allows me to check for warnings. Five years ago this questions has been posted but I haven't found any answer. Thanks

Re: [R] Testing for warning inside functions

2005-05-18 Thread Prof Brian Ripley
On Wed, 18 May 2005, Peter Wolf wrote: I am looking for a way to get a warning message immediately after an evaluation within a function. To get error messages you can use geterrmessage(). Well, only in an error handler, as an error normally throws you out of the function. (That is the point of

[R] fast matrix update

2005-05-18 Thread vincent
Hello, I use extensively m[i0:i1,j0:j1] = ... to update matrices parts. This writing is very convenient, but, since I'm doing this calculus many many times, I would like to know if there is a way (a hope ?) to do the same operation faster ? Thanks for any advice or pointer. Vincent

Re: [R] Testing for warning inside functions

2005-05-18 Thread Peter Wolf
Prof Brian Ripley wrote: On Wed, 18 May 2005, Peter Wolf wrote: I am looking for a way to get a warning message immediately after an evaluation within a function. To get error messages you can use geterrmessage(). Well, only in an error handler, as an error normally throws you out of the

Re: [R] fast matrix update

2005-05-18 Thread Duncan Murdoch
vincent wrote: Hello, I use extensively m[i0:i1,j0:j1] = ... to update matrices parts. This writing is very convenient, but, since I'm doing this calculus many many times, I would like to know if there is a way (a hope ?) to do the same operation faster ? It tends to be faster to use vector

[R] Re: Finding the right number of clusters

2005-05-18 Thread Adriano S. Melo
Dear Philip, Perhaps this reference might be useful: V. D. P. Pillar. How sharp are classifications? Ecology 80:2508-2516, 1999. Adriano S. Melo -- Message: 9 Date: Tue, 17 May 2005 08:42:02 -0400 From: Philip Bermingham [EMAIL PROTECTED] Subject: [R] Finding the right

[R] Bry Boschan routines in R

2005-05-18 Thread Lapointe, Pierre
Hello, Does anyone of you know if someone has programmed the Bry Boschan routines in R? It's also known as the NBER method for identifying economic cycles peaks and troughs. Or do you know any method in R for indentifying peaks and troughs for times series. I tried turnpoints(), but I get too

RE: [R] R -SQL

2005-05-18 Thread Gesmann, Markus
Your code looks more like Visual Basic rather than R. What you want is: for(j in 1:length(criteria$Title)){ sqlstring - paste(select q.type,crit.Title, r.Value from criteria crit, reply r,question_reply qr, question q, question_criteria qc, form_question fq where qr.reply=r.ID and

[R] standardization

2005-05-18 Thread Philip Bermingham
SAS Enterprise Miner recommendeds to standardize using X / STDEV(X) versus [X mean(X)] / STDEV(X) Any thoughts on this? Pros Cons Philip __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] R -SQL

2005-05-18 Thread Prof Brian Ripley
You have not told us what package you are using (and I think this is a repeat of a post also missing that crucial information). However, the main problem appears to be that you think operates on character strings. You have not told us your intentions, which might have been paste(select ...

RE: [R] DEV2bitmap: jpeg with res=400 not enough for CORELDRAW poster A0

2005-05-18 Thread Prof Brian Ripley
On Wed, 18 May 2005, Jan Verbesselt wrote: Thanks for the input! Finally I used: dev2bitmap(name,type=pdfwrite,height=8,width=14,res=1200) == resolution and colors were fine on the A0 poster! (and files are really small). Why not use R's native PDF driver if PDF is acceptable? This is a

Re: [R] standardization

2005-05-18 Thread Peter Dalgaard
Philip Bermingham [EMAIL PROTECTED] writes: SAS Enterprise Miner recommendeds to standardize using X / STDEV(X) versus [X – mean(X)] / STDEV(X) Any thoughts on this? Pros Cons When??? This makes absolutely no sense out of context. -- O__ Peter Dalgaard Blegdamsvej 3

Re: [R] standardization

2005-05-18 Thread Barry Rowlingson
Peter Dalgaard wrote: SAS Enterprise Miner recommendeds to standardize using X / STDEV(X) versus [X mean(X)] / STDEV(X) This makes absolutely no sense out of context. To paraphrase Tanenbaum: The nice thing about standardization is that there's so many ways to do it. Baz [[ Free On-line

RE: [R] applying a function over an array

2005-05-18 Thread McGehee, Robert
This will do it. a - array(1:12, c(4, 3)) d - function(a,b,c) {a+b+c} apply(a, 1, function(x) do.call(d, sapply(x, list))) [1] 15 18 21 24 -Original Message- From: BJ [mailto:[EMAIL PROTECTED] Sent: Wednesday, May 18, 2005 9:10 AM To: r-help@stat.math.ethz.ch Subject: [R] applying a

Re: [R] standardization

2005-05-18 Thread Peter Dalgaard
Barry Rowlingson [EMAIL PROTECTED] writes: The nice thing about standards is that there are so many of them to choose from, Curiously enough, the same quote came up today on dk.edb.system.unix in the context of translations. -- O__ Peter Dalgaard Blegdamsvej 3 c/

Re: [R] standardization

2005-05-18 Thread TEMPL Matthias
My thoughts on this is: Do not trust what SAS say´s and least of all what the Enterprise Miner said. Robust Statisticians recommendends to standardize using e.g. (X - median(X)) / ( MAD(X) / 0.675 ) Best, Matthias SAS Enterprise Miner recommendeds to standardize using X / STDEV(X) versus

Re: [R] applying a function over an array

2005-05-18 Thread Gabor Grothendieck
On 5/18/05, BJ [EMAIL PROTECTED] wrote: Is there a way to apply a function with several arguements over an array? For instance if you had a function d-function(a,b,c) {a+b+c} and a 4,3 array, could you apply the function over each line of the array? tapply seems to only allow one argument, and

Re: [R] align

2005-05-18 Thread Gabor Grothendieck
Check out the 'locf' function in the 'its' package and the 'na.locf' function in the 'zoo' package. On 5/18/05, Omar Lakkis [EMAIL PROTECTED] wrote: Is there a function in R that is similar to Splus's align? The idea is, if I have a data.frame, or an its object that is like this:

Re: [R] applying a function over an array

2005-05-18 Thread Gabor Grothendieck
On 5/18/05, Gabor Grothendieck [EMAIL PROTECTED] wrote: On 5/18/05, BJ [EMAIL PROTECTED] wrote: Is there a way to apply a function with several arguements over an array? For instance if you had a function d-function(a,b,c) {a+b+c} and a 4,3 array, could you apply the function over each line

[R] 'fitdistr' and two views of the same data?

2005-05-18 Thread Matt Jadud
Hello, I have detailed (with pictures and whatnot) my question on my weblog at http://www.cs-ed.org/blogs/mjadud/archives/2005/05/a_question_abou.html The short version of the question is this: When I ask 'fitdistr' to try and fit my distribution as a weibull distribution, it comes up with

RE: [R] standardization

2005-05-18 Thread bogdan romocea
You asked another question about clustering, so I presume you want to standardize some variables before clustering. In SAS, PROC STDIZE offers 18 standardization methods. See http://support.sas.com/91doc/getDoc/statug.hlp/stdize_sect12.htm#stat_stdize_stdizesm for details. If you're really

Re: [R] align

2005-05-18 Thread Gabor Grothendieck
Its in 'its' but due to the use of namespaces in 'its' one must do this to access it directly: its:::most.recent On 5/18/05, Omar Lakkis [EMAIL PROTECTED] wrote: From the definition of locf locf function (x) { if (!inherits(x, its)) stop(function is only valid for objects of

Re: [R] How to convert array to c()

2005-05-18 Thread Muhammad Subianto
Dear all, Many thanks for your help. Regards, Muhammad Subianto On this day 5/18/2005 4:57 PM, Liaw, Andy wrote: Is this what you want? split(a, row(a)) $1 [1] 1 5 9 $2 [1] 2 6 10 $3 [1] 3 7 11 $4 [1] 4 8 12 Andy On this day 5/18/2005 5:15 PM, [EMAIL PROTECTED] wrote: Look

Re: [R] How to convert array to c()

2005-05-18 Thread Li, Jia
How about? p1-a[1,] p1 [1] 1 5 9 p2-a[2,] p2 [1] 2 6 10 p3-a[3,] p3 [1] 3 7 11 p4-a[4,] p4 [1] 4 8 12 Jia - Original Message - From: Muhammad Subianto [EMAIL PROTECTED] To: r-help@stat.math.ethz.ch Sent: Wednesday, May 18, 2005 10:18 AM Subject: [R] How to convert array to c()

[R] Why can't I download window binary zip packages

2005-05-18 Thread Lisa Wang
Hello there, I tried so many times to download windows binary zip package but it told me that I don't have the access to do so. It worked for me a few months ago. Please help me with it. Thank you Lisa Wang Msc. Princess Margaret Hospital Toronto , Canada tel: (416) 946 4501 ext.5201

Re: [R] Why can't I download window binary zip packages

2005-05-18 Thread Duncan Murdoch
Lisa Wang wrote: Hello there, I tried so many times to download windows binary zip package but it told me that I don't have the access to do so. It worked for me a few months ago. Please help me with it. You need to talk to one of your system administrators. R won't tell you that you don't

Re: [R] lattice plots question

2005-05-18 Thread Sundar Dorai-Raj
Laura Holt wrote: Dear R People: Is there any way to have the background of lattice plots be white instead of grey, please? This is not a criticism by any means...the lattice stuff is UNbelievable! Thanks, Laura Holt mailto: [EMAIL PROTECTED] R Version 2.1.0 Windows Hi, Laura, See

[R] Advice requested: installation failing for foreign_0.8-7, mgcv_1.2-4, rpart_3.1-23, VR_7.2-15

2005-05-18 Thread Michael Kubovy
I'm running R Version 2.1.0 (2005-04-18), ISBN 3-900051-07-0 on Machine Model: Power Mac G5 CPU Type: PowerPC 970 (2.2) Number Of CPUs: 2 CPU Speed:2 GHz L2 Cache (per CPU): 512 KB Memory: 2 GB with software: System Version:Mac OS X 10.3.9 (7W98)

[R] SAMM package for mixed models

2005-05-18 Thread kwright
First, a disclaimer. I am not affiliatied with the SAMM package. I am only a user of the package, but I have been contacted (off list) by people requesting information about SAMM and so I am posting this information here. SAMM is software for fitting mixed models. Versions are available for

[R] loop

2005-05-18 Thread Omar Lakkis
Rather than using a loop, how can I remove all consequentially repeated values as in this example? I am guessing using diff would help but not quite sure how. get s date f 1 1999-01-01 1 2 1999-01-02 1 3 1999-01-03 1 4 1999-01-04 2 5 1999-01-05 2 v - s[1,'f']; for (i in

[R] Re: text mining: ttda

2005-05-18 Thread Weiwei Shi
Can anyone suggest some good text mining reference or books? thanks, weiwei On 5/18/05, Weiwei Shi [EMAIL PROTECTED] wrote: Hi, I am working on a text mining project and i am interested in ttda package. however, I really cannot find the document for this package in English. Can anyone give

[R] source-only package, but still: Error: package 'simple' was built for powerpc-apple-darwin7.9.0

2005-05-18 Thread Paul Shannon
I have a number simple R functions written for the biologists I work with. These functions will evolve, and the documentation will get steadily better. I hope to put them in a package on a local webserver, along with lots of help files, so that the users can easily update their installation.

Re: [R] loop

2005-05-18 Thread Prof Brian Ripley
On Wed, 18 May 2005, Omar Lakkis wrote: Rather than using a loop, how can I remove all consequentially repeated values as in this example? I am guessing using diff would help but not quite sure how. get s date f 1 1999-01-01 1 2 1999-01-02 1 3 1999-01-03 1 4 1999-01-04 2 5

[R] Fitting linear model to the matrix

2005-05-18 Thread Huhrik Beybutov
I'm afraid I need your help with linear models in R, I apologize in advance is this question was beaten to death. I'm trying to fit Y( i, j ), where i is row, j is column. My model is log [ Y( i, j ) ] = log A( i ) * B( j )+ C( i ) ] + error. I have i x j observations, I need to find 2*i + j

[R] from list to dataframe

2005-05-18 Thread sms13+
I was wondering if someone can help me figure out the following: I have two patient datasets, ds1 and ds2. ds1 has fields patid, date, and lab1. ds2 has patid, date, and lab2. I want to find all the patids that have at least 2 dated records for each lab. I started by splitting each dataset

Re: [R] Fortran 95 in R ?

2005-05-18 Thread Bill Northcott
On 18/05/2005, at 8:10 PM, Joel Bremson wrote: Is it possible to run Fortran 95 code from R? I don't think so, but hopefully someone can prove me wrong. This is not really any issue with R. It is matter of what compilers you have installed. If you have Fortran 95 compiler and you use it to

Re: [R] source-only package, but still: Error: package 'simple' was built for powerpc-apple-darwin7.9.0

2005-05-18 Thread Prof Brian Ripley
On Wed, 18 May 2005, Paul Shannon wrote: I have a number simple R functions written for the biologists I work with. These functions will evolve, and the documentation will get steadily better. I hope to put them in a package on a local webserver, along with lots of help files, so that the users

[R] Bias to do with search engines

2005-05-18 Thread Patrick Connolly
On Wed, 18-May-2005 at 02:55PM -0700, [EMAIL PROTECTED] wrote: | | First, a disclaimer. I am not affiliatied with the SAMM package. I am | only a user of the package, but I have been contacted (off list) by people | requesting information about SAMM and so I am posting this information | here.

[R] Call R from Fortran

2005-05-18 Thread Samiran Sinha
Hello, I need to call a R function from Fortran 77 program. How will I do that exactly? I will grately appreciate any help. Sincerely, -- Samiran Sinha, Ph.D Assistant Professor Department of Statistics Texas AM University TAMU 3143 College Station, TX-77843 Phone Number: (979) 845 2966(O)

Re: [R] Line width through plot size reduction

2005-05-18 Thread Paul Murrell
Hi Jacques VESLOT wrote: Dear R-users, Someone, who uses R under Mac, wants to insert a couple of small plots (each with several lines) in an article, but he has to reduce plots' size significantly. He did it (in pdf or enc. ps) but, unfortunately, everything is reduced but lines' width.

[R] Calculation of Durbin-Watson p-value

2005-05-18 Thread Ramesh Kolluru
Sir,   I am unable to get the source code for Durbin-Watson test, as I want to calculate the p-value for Durbin Watson statistic using interpolation method. I sent this mail to r-help, but it was rejected, please suggest me some way. I will be highly greatful to you. Thanks in advance Ramesh