Chris,
I haven't seen anyone post a reply yet so thought I'd
throw in my thoughts. I'm no R expert!
When you talk about an exponential trend line are you
refering to:
1) y=ax^b
or
2) y=ae^(bx)
If 1) then take base10 logs of y and x and then fit them
with simple linear regression. Then
As I said, I am new to R after spending far too many years using SAS. I'm
slowly getting the hang of R and like it very much.
Thanks for your insights and help.
Murray M Cooper, Ph.D.
Richland Statistics
9800 N 24th St
Richland, MI, USA 49083
Mail: richs...@earthlink.net
Dear All,
I am new to R and slowly learning how to use the system.
The following code is an exercise I was trying.
The intent is to generate 10 random samples of size 5 from
a vector with integers 1:10 and 2 missing values. I then want
to generate a matrix, for each sample which shows the
Paul,
I suggest looking up observer agreement. The description of your study
sounds like a classical
categorical observer agreement problem. I can't give
a reference off the top of my head, but if you get
stuck, e-mail me and I'll try and find a ref to get you started.
Murray M Cooper, Ph.D.
You don't say what the error was, for the R factorial function,
but it is probably irrelevant for your question.
Factorials get to be big numbers rather quickly and unless you
are using a program that does arbitrary precission arithmetic
you will quickly exceed the precission limits, for storing
For science yes. For pleasure I'll still take a pint instead of 570ml!
Murray
- Original Message -
From: Rolf Turner r.tur...@auckland.ac.nz
To: Emmanuel Charpentier charp...@bacbuc.dyndns.org
Cc: r-h...@stat.math.ethz.ch
Sent: Friday, April 03, 2009 6:18 PM
Subject: Re: [R] [OT ?]
The usual approach is to count the co-occurence within so many words of each
other.
Typical is between 5 words before and 5 words after a given word.
So for each word in the document, you look for the occurence of all other
words
within -5 -4 -3 -2 -1 0 1 2 3 4 5 words. Depending on the
Just a word of caution. Having done a lot of work with 24 hour blood
pressure and ecg recordings, these series are seldom stationary which
presents problems with spectral analysis. I don't know what your ultimate
goal is, but in my work I found it often better to work with subsets, of the
I'm deeply disappointed! I keep checking the mail list to see if you guys
are posting answers to questions I haven't asked yet. It would save me a lot
of time!
Best,
Murray
- Original Message -
From: Rolf Turner r.tur...@auckland.ac.nz
To: Sherri Heck sh...@ucar.edu
Cc:
Charlotta,
I'm not sure what you mean when you say simple linear
regression. From your description you have two groups
of people, for which you recorded contaminant concentration.
Thus, I would think you would do something like a t-test to
compare the mean concentration level. Where does the
samples.
Murray M Cooper, Ph.D.
Richland Statistics
9800 N 24th St
Richland, MI, USA 49083
Mail: richs...@earthlink.net
- Original Message -
From: David Winsemius dwinsem...@comcast.net
To: Murray Cooper myrm...@earthlink.net
Cc: Charlotta Rylander z...@nilu.no; r-help@r-project.org
Sent
David,
If you really want to do a test on this data, I would suggest
a Fisher's Exact test, but you want to use hypergeometric
probabilities. You would probably want to try the CMH
test, if the function allows a single table and actually uses
hypergeometric probabilities.
My suggestion, would
A recent thread on summary statistics, got me thinking. (Note this may not
happen often.) A function that would do summaries as describe below (similar
to SAS PROC UNIVARIATE) might be a nice addition to the main R system. Is
there a process by which functions, from packages can eventually be
Nidhi,
You neglect to say what your OS is.
If you are using Windows, see the R for Windows FAQ faq 2.8. I haven't
looked, but I'm sure there is a similar item for other OS.
Murray M Cooper, Ph.D.
Richland Statistics
9800 N 24th St
Richland, MI, USA 49083
Mail: richs...@earthlink.net
-
David,
This may be a case of If all you have is a hammer, everything looks like a
nail.
If all you want to do is remove the last line if it contains a CONTROL-Z,
why
not use something like perl to process the files?
Murray M Cooper, Ph.D.
Richland Statistics
9800 N 24th St
Richland, MI, USA
Consider yourself lucky!
I'm sure there are many people who would
prefer not to see their name in the NYT. ;-)
Murray Coooper
- Original Message -
From: Duncan Murdoch murd...@stats.uwo.ca
To: Mark Difford mark_diff...@yahoo.co.uk
Cc: r-help@r-project.org
Sent: Thursday, February 05,
Categorical data analysis is definitely the way you want to go.
Which test you use depends on how you are going to use the results.
For quick and dirty I would suggest using Fisher's exact test on all
2x2 submatricies of counts. In this case, with 4 treatments you have
6 possible 2x2
I am new to R, so maybe I'm missing the point of your question. But why
wouldn't you just use sum(a,b)?
Murray M Cooper, Ph.D.
Richland Statistics
9800 N 24th St
Richland, MI, USA 49083
Mail: richs...@earthlink.net
- Original Message -
From: Fuchs Ira irafu...@gmail.com
To:
I was about to post a similar reply.
Stavros's reply was very eloquent and should be taken to heart!
Murray M Cooper, Ph.D.
Richland Statistics
9800 N 24th St
Richland, MI, USA 49083
Mail: richs...@earthlink.net
- Original Message -
From: Stavros Macrakis macra...@alum.mit.edu
To:
I am working on a methodology for qualifying R, for GLP and GCP.
If I quailfy only the base R install, with no contributed packages, it
seems relatively simple to qualify R. However, from time to time I
will want to use a contributed package. If I use a contributed package,
does it leave
but been unable to locate the
file which contains model.matrix.
Thanks for your help.
Murray Cooper
- Original Message -
From: Duncan Murdoch murd...@stats.uwo.ca
To: zhijie zhang rusers...@gmail.com
Cc: r-h...@stat.math.ethz.ch
Sent: Saturday, January 17, 2009 6:00 AM
Subject: Re: [R] Where
Do you really need the p-value or do you want to test at one
of the socially acceptable levels (i.e. .05 or .01). If all you want
is the test, use:
quantile(bootsample,c(0.025,0.975))
If the quantile range includes 0 then you decide there is no evidence
that the mean is different from zero, at
Mostly I am interested in using R for statistics.
I am also interested in being able to look at source code.
I hope to be able to write extensions.
I tried the suggestion below but was unable to access vcov.lm.
methods(vcov)
[1] vcov.Arima* vcov.glm* vcov.lm*vcov.mlm* vcov.nls*
23 matches
Mail list logo