[R] CircStats and Anova

2005-05-15 Thread Mueller, Adrienne
Hi, If I have two sets of directional data (in radians) and want to compare them with a multifactorial anova. Is it even legitimate to compare circular data with an anova? The books I've picked up from the library don't really say, but it looks unlikely. If it is allowable, is my having stored

[R] testing equality of covariance matrices

2005-05-15 Thread K. Steinmann
Dear R-mailers, I would like to test for equality of population covariance matrices. But I can't find a R tool to do so. I saw, that other people had the same question, but I could not find an answer to it, I would appreciate to know the missed link. Thank you, b.w. K. Steinmann

Re: [R] what means Best sample in clustering using Clara?

2005-05-15 Thread Uwe Ligges
Amir Safari wrote: Hi All, Can some body tell what means Best Sample in clustering using Clara? ?clara.object tells you: sample: labels or case numbers of the observations in the best sample, that is, the sample used by the clara algorithm for the final partition. Uwe Ligges

Re: [R] Job Opportunity: Statistical Guru CC 083

2005-05-15 Thread Uwe Ligges
BTW: I am looking for a wife that is very intelligent, earns at least 300,000$ a year, is less than 25 years old, is a real beauty, happily cooks, cleans, and works in the garden, and is never moody, of course. She would have the great opportunity to solve one of the most complex challenges

[R] Re: [CSL #265590] Schily ueber Deutschland

2005-05-15 Thread lab
Greetings. (This is an automated response. There is no need to reply.) Your message regarding: Schily ueber Deutschland has been received and assigned a request number of 265590. In order help us track the progress of this request, we ask that you note the following: 1) Please save this

[R] Save parameters from optim during iteration procedure

2005-05-15 Thread Camarda, Carlo Giovanni
Dear R-users, I am going to try to be as clearer as possible, showing also an example. 1) I have a function (in my real case it's much more complex) 2) I use optim to minimize 3) I want to use as method L-BFGS-B for several reasons 4) I know I could use trace=6 (in control) in order to see live

[R] Suppressing warning messages

2005-05-15 Thread Tolga Uzuner
How do I suppress the following ? Warning messages: 1: the condition has length 1 and only the first element will be used in: if (strike == forward) atmvol(forward, t, alpha, beta, rho, upsilon) else { 2: the condition has length 1 and only the first element will be used in: if (x(z) == 0) 1

Re: [R] Suppressing warning messages

2005-05-15 Thread Kjetil Brinchmann Halvorsen
Tolga Uzuner wrote: How do I suppress the following ? Warning messages: 1: the condition has length 1 and only the first element will be used in: if (strike == forward) atmvol(forward, t, alpha, beta, rho, upsilon) else { 2: the condition has length 1 and only the first element will be used in:

Re: [R] Suppressing warning messages

2005-05-15 Thread Gabor Grothendieck
On 5/15/05, Tolga Uzuner [EMAIL PROTECTED] wrote: How do I suppress the following ? Warning messages: 1: the condition has length 1 and only the first element will be used in: if (strike == forward) atmvol(forward, t, alpha, beta, rho, upsilon) else { 2: the condition has length 1 and

Re: [R] Suppressing warning messages

2005-05-15 Thread Uwe Ligges
Kjetil Brinchmann Halvorsen wrote: Tolga Uzuner wrote: How do I suppress the following ? Warning messages: 1: the condition has length 1 and only the first element will be used in: if (strike == forward) atmvol(forward, t, alpha, beta, rho, upsilon) else { 2: the condition has length 1 and only

RE: [R] Job Opportunity: Statistical Guru CC 083

2005-05-15 Thread Mike Waters
Wow, Cutting edge to bleeding heart - and all in one day! 8¬ Mike P.s. Uwe - start with an easy one, like a Unicorn. -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of Uwe Ligges Sent: 15 May 2005 12:00 To: Don Alexander Cc: r-help@stat.math.ethz.ch

Re: [R] Save parameters from optim during iteration procedure

2005-05-15 Thread Uwe Ligges
Camarda, Carlo Giovanni wrote: Dear R-users, I am going to try to be as clearer as possible, showing also an example. 1) I have a function (in my real case it's much more complex) 2) I use optim to minimize 3) I want to use as method L-BFGS-B for several reasons 4) I know I could use trace=6 (in

[R] Seeking friend for life (not)

2005-05-15 Thread Chris Evans
I'll take a risk following Uwe's wonderful response to that advert. I'm looking for someone, perhaps particularly a stats student, who might want to do a piece of work with me on using R to present and analyse routine data that psychotherapists might submit on a cgi-bin interface and perhaps

RE: [R] Job Opportunity: Statistical Guru CC 083

2005-05-15 Thread Ted Harding
On 15-May-05 Mike Waters wrote: Wow, Cutting edge to bleeding heart - and all in one day! 8¬ Mike P.s. Uwe - start with an easy one, like a Unicorn. Dear Uwe, Some further advice. I fear your vacancy notice has not been written to attract the best candidate. There are even important

[R] adjusted p-values with TukeyHSD?

2005-05-15 Thread Christoph Strehblow
hi list, i have to ask you again, having tried and searched for several days... i want to do a TukeyHSD after an Anova, and want to get the adjusted p-values after the Tukey Correction. i found the p.adjust function, but it can only correct for holm, hochberg, bonferroni, but not Tukey. Is it

Re: [R] CircStats and Anova

2005-05-15 Thread Spencer Graves
It's not clear to me what you are asking. Have you considered aov.circular in package circular? The following example comes from the help page for that function: x - c(rvonmises(50, 0, 1), rvonmises(100, pi/3, 10)) group - c(rep(0, 50), rep(1, 100)) aov.circular(x,

[R] Not sure if this is aggregate or some other task.

2005-05-15 Thread David L. Van Brunt, Ph.D.
I have data where where I've taken some measurements three times... twice in rapid succession so I could check test-retest reliability of a piece of equipment, and then a third measurement some time later. Not I'd like to do an analysis where I have two scores... the first being the mean of

[R] Finding Gradient from Contour map

2005-05-15 Thread Laura Quinn
I have a detailed map matrix in R (x,y,z co-ords), and am wishing to calculate the physical gradient at a number of points on the map wrt different axes and averaged over a number of distances (i.e. slope in N-S direction averaged over 10m, 50m, 100m respectively). Is there any function within R

Re: [R] How to make label in multi plot

2005-05-15 Thread Paul Murrell
Hi (cc'ed to Pierre Lapointe because this should answer the question about [R] Centered overall title with layout() as well) Muhammad Subianto wrote: Dear R-Help, As a reference about multi plot, http://finzi.psych.upenn.edu/R/Rhelp02a/archive/48725.html I want to know how can I make a label for

[R] Mann-Whitney Wilcoxon Rank Sum

2005-05-15 Thread Jim BRINDLE
Hello, I am hoping someone could shed some light into the Wilcoxon Rank Sum Test for me? In looking through Stats references, the Mann-Whitney U-test and the Wilcoxon Rank Sum Test are statistically equivalent. When using the following dataset: m -

[R] row.names need reordering

2005-05-15 Thread Anders Schwartz Corr
Hi, The row.names in my matrix seem to be out of order. I don't remember putting row.names in in the first place, I don't see what use they are, and they are out of order (perhaps because I sorted them at one point when the data was in data.frame format). Can I delete the rownames? or at least

[R] turning labels into a vector

2005-05-15 Thread Briggs, Meredith M
Hello 1/ 'priors' is a table looking like: W123 T678 S789 23 42 11 12 35 9 etc 2/ WBS - labels(priors) gives me a result of class list and length 1 looking like: W123 T678 S789 I want to read W123 into X[1] as W, T687 into X[2] as T and S789 into X[3]

Re: [R] Mann-Whitney Wilcoxon Rank Sum

2005-05-15 Thread Prof Brian Ripley
On Mon, 16 May 2005, Jim BRINDLE wrote: Hello, I am hoping someone could shed some light into the Wilcoxon Rank Sum Test for me? In looking through Stats references, the Mann-Whitney U-test and the Wilcoxon Rank Sum Test are statistically equivalent. Yes, but not numerically: they differ by a