I didn't set the variable "treat" as ordered factor,
just ordinary factor. In my example, the 3 levels of
"treat" are of equal distance in the order of from a
to b to c. So my understanding is that a contrast in
the form of "a+c-2*b=0" was to test a linear trend of
the response variable among the 3
Hi,
Where do I find the link to install RDCOM? I need to use RExcel for my
project.
Thanks a lot!
Uma
This e-mail may contain confidential and/or privileged i...{{dropped:13}}
__
R-help@r-project.org mailing list
https://stat.ethz.ch
If I understand you correctly, use outer() -- the for loops suggested below
are not the "R way" (nor do they seem to fully address the "and so on" part
of your request):
ratios <- outer(z,z,"/")
This produces a matrix, the first column of which is z/z[1], the second of
which is z/z[2], and so for
"Charles C. Berry" <[EMAIL PROTECTED]> wrote in
news:[EMAIL PROTECTED]:
> On Wed, 16 Jan 2008, David Winsemius wrote:
>
>> I am a physician examining an NHANES dataset available at the
>> NCHS website:
http://www.cdc.gov/nchs/about/major/nhanes/nhanes2005-2006/demo_d.xpt
snip
>>
>> TC.ran <-
אשה במודיעין
גידלה חמור בקומה שלישית בבית משותף
סרטון מדהים
להעביר הלאה
http://www.go443.co.il/Go443/Video/Video.asp __
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.o
I don't quite follow what you are trying to do but the second contrast
has a few interpretations with the same meaning in your case
1) are the 2-1 and 3-2 differences equal
2) lack of fit of a linear trend
3) is there a quadratic response
If you declare your factor to be "ordered" then the de
James Reilly wrote:
>> Max Kuhn wrote:
>>> Factors have huge benefits over character data in SAS. For a series
>>> regulatory filings, I had miles of SAS code to compute KxK tables
>>> where all the cells must show up. For example, if one of the levels of
>>> one of the variables was never observed
Hi, I am running a simple one-way ANOVA with an
independent factot variable "treat" (3 levels: a, b
and c) and a response variable "y". I want to test a
linear relationship of the response among the 3 levels
of the variable "treat" (ordered a->b->c). I used
glht() from multcomp package. Later I fou
Johannes Graumann wrote:
I really do not know ho to else title this ... I want to draw something like
the attached png with R and would like to poll you on how to start ... make
an empty plot first and then start positioning the characterstring
by 'text' and then drawing the lines ...
Joh
Joha
hi Juan,
It is not so elegant, but work fine. I know that our colleagues can do it on a
simple line.
z<-c(526,723,110,1110,34,778,614,249,14)
v1<-NULL
v2<-NULL
for (i in 1:(length(z)-1))
{
for (j in i:length(z))
{
v1<-rbind(v1,z[i])
v2<-rbind(v2,z[j])
}
}
df<-data.frame(cbind(v1=v1,v2=v2))
nam
> Max Kuhn wrote:
>> Factors have huge benefits over character data in SAS. For a series
>> regulatory filings, I had miles of SAS code to compute KxK tables
>> where all the cells must show up. For example, if one of the levels of
>> one of the variables was never observed, the corresponding row
Thank you very much...
That was helpful..
On Jan 15, 2008 12:58 AM, Charles C. Berry <[EMAIL PROTECTED]> wrote:
> On Mon, 14 Jan 2008, Marko Milicic wrote:
>
> > Dear all,
> >
> > I'm trying to process HUGE datasets with R. It's very fast, but I would
> like
> > to optimize it a bit more, by focu
I am using gls to fit a linear model with spatially-autocorrelated errors.
My first step is to fit a simple model, and inspect a semivariogram of
residuals.
The following example gives this error:
Error in FUN(X[[1L]], ...) : unused argument(s) (method = "euclidean")
#Example
library(nlme)
x <- r
On Thursday 17 January 2008 10:55:03 am Zembower, Kevin wrote:
> I've been given the job of extracting some data from the United States
> 2000 census (files at
> http://www2.census.gov/census_2000/datasets/Summary_File_2/Maryland/all_
> Maryland.zip 52M). I'm only interested in Census Block Groups
I'm afraid I can't follow your examples, but you seem to me to be
mixing contingency table tests and goodness of fit tests in a somewhat
incoherent fashion.
Note that your ``x2()'' function does a contingency table test, and not
a goodness of fit test.
Note that in chisq.test(), if ``x'' is a ma
It looks fine to me. Try str(d) and check to be sure
that Votes is a numeric value or integer value.
I ran this code with no problem.
x <- "Name Votes
John 300
Sean222
Andy 467
Sinead 740
David 124
James 641
William 380 "
d <- read.table(textConnection(x),
Dear R-helpers,
A student in my Department reports that his R is in Japanese:
"I lived in Japan and have a US-purchased mac with the Japanese input/
reader enabled, which must be where the Japanese came from. When I
uninstalled R and reinstalled it last time (I think I completely
uninstalled
R Help on 'chisq.test' states that
"if 'simulate.p.value' is 'TRUE', the p-value is computed by Monte
Carlo simulation with 'B' replicates.
In the contingency table case this is done by random sampling from
the set of all contingency tables with given marginals, and works
o
I've been given the job of extracting some data from the United States
2000 census (files at
http://www2.census.gov/census_2000/datasets/Summary_File_2/Maryland/all_
Maryland.zip 52M). I'm only interested in Census Block Groups (CBGs)
located within Baltimore City, Maryland. Additionally, I just ha
I am using the auto.arima package to do some basic forecasting based
on CPU usage. I now have found a calendar that has various activities
that partially control the computer's usage and want to factor that in
(They are effectively dummy variables indicating a particular type of
activity that week)
On Jan 17, 2008, at 9:44 AM, Johannes Graumann wrote:
> I really do not know ho to else title this ... I want to draw
> something like
> the attached png with R and would like to poll you on how to
> start ... make
> an empty plot first and then start positioning the characterstring
> by 'tex
Thanks to all!
This is really helpful!
Sincerely,
Erin
On Jan 17, 2008 12:00 PM, Marc Schwartz <[EMAIL PROTECTED]> wrote:
> Erin Hodgess wrote:
> > Hi R People:
> >
> > I'm reading "Statistical Computing with R", by Maria Rizzo, and it's
> > really good.
> >
> > Anyhow, I have a question about
Clint Bowman wrote:
> So how does SAS compare with one of the specialty languages such as perl.
> I've found the combination of perl and R to work quite satisfactorily (as
> long as I don't confuse the syntax and functions available in each.)
Now that the topic has drifted off the subject of wha
Thanks Hadley. Witty! Profound! Concise! I think this is definitely a
Fortunes candidate.
-- Bert Gunter
Genentech Nonclinical Statistics
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On
Behalf Of hadley wickham
Sent: Thursday, January 17, 2008 9:56 AM
To: Wittner,
Wittner, Ben, Ph.D. wrote:
> Several people have mentioned large, messy data sets.
> I am curious as to in what way messy data sets are messy.
> (I am also curious about what SAS does that helps one deal with them, but
> perhaps that's asking too much.)
>
One aspect is that in the "SAS culture"
On Jan 17, 2008 9:53 AM, Thompson, David (MNR)
<[EMAIL PROTECTED]> wrote:
> Hello Hadley,
>
> I am trying to reproduce the following with ggplot:
> a <- seq(0, 360, 5)*pi/180 ; a
> ac <- sin(a + (45*pi/180)) + 1 ; ac
> plot(a, ac, type='b', xaxt = "n")
> axis(1, at=seq(0,6,1), label
Erin Hodgess wrote:
> Hi R People:
>
> I'm reading "Statistical Computing with R", by Maria Rizzo, and it's
> really good.
>
> Anyhow, I have a question about something in there.
>
>> u<- runif(5)
>> u
> [1] 0.1177041 0.4271790 0.4601597 0.2204846 0.4051473
>> #in the book
>> sum(as.integer(u> 0.4
"Happy families are all alike; every unhappy family is unhappy in its own way."
Leo Tolstoy
and every messy data is messy in its own way - it's easy to define the
characteristics of a clean dataset (rows are observations, columns are
variables, columns contain values of consistent types). If you
Erin Hodgess wrote:
> Hi R People:
>
> I'm reading "Statistical Computing with R", by Maria Rizzo, and it's
> really good.
>
> Anyhow, I have a question about something in there.
>
>
>> u <- runif(5)
>> u
>>
> [1] 0.1177041 0.4271790 0.4601597 0.2204846 0.4051473
>
>> #in the book
>> sum
Hi R People:
I'm reading "Statistical Computing with R", by Maria Rizzo, and it's
really good.
Anyhow, I have a question about something in there.
> u <- runif(5)
> u
[1] 0.1177041 0.4271790 0.4601597 0.2204846 0.4051473
> #in the book
> sum(as.integer(u > 0.4))
[1] 3
> #what I would do
> sum(u
So how does SAS compare with one of the specialty languages such as perl.
I've found the combination of perl and R to work quite satisfactorily (as
long as I don't confuse the syntax and functions available in each.)
Clint
Clint BowmanINTERNET: [EMAIL PROTECTED]
Air Disp
Thnk you very much! It now works correctly!
Much Appreciated,
John.
hoogeebear wrote:
>
> Hi,
>
> I am having trouble with an error I keep getting. I am just trying to
> create a simple pic chart from a small table. Hope someone can help. I am
> new to R.
>
> Table:
> Name Votes
> John
Would it be possible for Matthew (the original person) to tell us what
he ended up with for his final talk, please?
Thanks,
Erin
On Jan 17, 2008 10:45 AM, Wittner, Ben, Ph.D.
<[EMAIL PROTECTED]> wrote:
> Several people have mentioned large, messy data sets.
> I am curious as to in what way messy
For some unknown reason I stopped receiving any messages from the R-
help mailing list. See if this test gets through.
Thanks,
Gang
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http
> d <- read.table("C:\\rep.csv", head=TRUE, sep=",")
> > pie(d$Votes,
> + labels=d$Name,
> + main="Class Rep Results\n(Final Results)")
>
> Error:
> Error in pie(d$votes, labels = d$name, main = "Class Rep Results\n(Final
> Results)") :
> 'x' values must be positive.
The first input to the pie
Several people have mentioned large, messy data sets.
I am curious as to in what way messy data sets are messy.
(I am also curious about what SAS does that helps one deal with them, but
perhaps that's asking too much.)
Thanks.
-Ben
> -Original Message-
> From: [EMAIL PROTECTED] [mailto:[E
Hi,
I am having trouble with an error I keep getting. I am just trying to create
a simple pic chart from a small table. Hope someone can help. I am new to R.
Table:
Name Votes
John 300
Sean222
Andy 467
Sinead 740
David 124
James 641
William 380
Commands:
d <
The argument for SAS (and Stata) when working with large dataset comes
up fairly often. I have not had much experience in this area, but have
been pleasantly surprised using R in combination with an SQL interface,
in situations with modestly large, messy datasets. I certainly would
appreciate
See ?outer
outer(Z, Z, function(x,y) x/y)
Gabor
On Thu, Jan 17, 2008 at 01:24:33PM -0300, Juan Pablo Fededa wrote:
> Dear Contributors:
>
> I have the next vector:
>
> "Z"
>
> 526
> 723
> 110
> 1110
> 34
> 778
> 614
> 249
> 14
>
> I want to generate a vector containing the ratios of all the
Dear Contributors:
I have the next vector:
"Z"
526
723
110
1110
34
778
614
249
14
I want to generate a vector containing the ratios of all the values
versus all the values of the z vector. I mean a vector containing the
values of 526/723, 526/110, and so on, 723/723, 723/110, and so on,
and so
Thompson, David (MNR) wrote:
> Hello Hadley,
>
> I am trying to reproduce the following with ggplot:
> a <- seq(0, 360, 5)*pi/180 ; a
> ac <- sin(a + (45*pi/180)) + 1 ; ac
> plot(a, ac, type='b', xaxt = "n")
> axis(1, at=seq(0,6,1), labels=round(seq(0,6,1)*180/pi),1)
> abline(v=
A new major upgrade of the Epi package for Epidemiological data analysis
has been put on CRAN, it is now at version 1.0.7.
It contains an entirely new way of representing follow-up data on
multiple timescales and multiple states. See the function Lexis().
Plus a lot of other useful stuff for epid
Max Kuhn wrote:
> Factors have huge benefits over character data in SAS. For a series
> regulatory filings, I had miles of SAS code to compute KxK tables
> where all the cells must show up. For example, if one of the levels of
> one of the variables was never observed, the corresponding row or
> co
Splendid, thanks for your quick response.
[EMAIL PROTECTED] wrote: > I have a data frame column in which I would like to
replace some
> of the numbers dependent on their value.
>
> data frame = zz
>
> AveExpr t P.Value FC
> 7.481964 7.323950 1.778503e-04 2.218760
> 7.585783 12.233056 6.679776
[EMAIL PROTECTED] napsal dne 17.01.2008 15:08:40:
> There must be a better way but this will do it for
> you.
>
> x <- runif(100, 0, 1)
> y <- runif(100, 0, 1)
> z <- data.frame(x,y)
>
> plot(subset(z, z$y >=.5), col="red", ylim=c(min(z$y),
> max(z$y)), pch=16)
> points(subset(z, z$y <=.49), co
Hello Hadley,
I am trying to reproduce the following with ggplot:
a <- seq(0, 360, 5)*pi/180 ; a
ac <- sin(a + (45*pi/180)) + 1 ; ac
plot(a, ac, type='b', xaxt = "n")
axis(1, at=seq(0,6,1), labels=round(seq(0,6,1)*180/pi),1)
abline(v=c(45*pi/180, 225*pi/180))
I can get the bas
> The help says that the exact method is computationally demanding, but even
> after days the computing it won't finish. Also, if I include a frailty-term,
> the exact method gives me results in no time. Is my setup incorrect?
Assume that at some particular time point there are k deaths and n su
You might want to have a look at the recode function
in the car package. By the way I think you meant
26-35 not 25-25.
===
Example
xx <- data.frame(age=c(25, 33, 22, 19,21, 30, 32,
31),
edu=c(2,5 ,3, 1,3, 4, 4, 1))
library(car)
a
I really do not know ho to else title this ... I want to draw something like
the attached png with R and would like to poll you on how to start ... make
an empty plot first and then start positioning the characterstring
by 'text' and then drawing the lines ...
Joh
<>___
There must be a better way but this will do it for
you.
x <- runif(100, 0, 1)
y <- runif(100, 0, 1)
z <- data.frame(x,y)
plot(subset(z, z$y >=.5), col="red", ylim=c(min(z$y),
max(z$y)), pch=16)
points(subset(z, z$y <=.49), col="blue", pch=15)
--- dxc13 <[EMAIL PROTECTED]> wrote:
>
> useR's
>
> Previously I used SAS for 23 years and now R/S-Plus for 17. SAS is
> effective for large datasets (in my work > 500,000 subjects) but except
> for that, R is far superior to SAS for data management and manipulation.
> Just four of the reasons are that you can
>
> - merge data frames multiple w
Walter Paczkowski wrote:
> Good morning,
>
> I use SAS and R/S-Plus as my primary tools so I have a lot of experience with
> these programs. By far and away, SAS is superior for handling the "messy"
> datasets, but also the very large ones. I work at times with datasets in the
> hundreds of t
Rob Robinson wrote:
> I wonder if those who complain about SAS as a programming environment have
> discovered SAS/IML which provides a programming environment akin to Matlab
> which is more than capable (at least for those problems which can be treated
> with a matrix like approach). As someone who
Brian,
On 16 January 2008 at 11:26, Brian O'Gorman wrote:
| I'm having trouble with R CMD INSTALL Rmpi_0.5-5.tar.gz
| --configure-args=~/lam
|
| lam is is installed locally.
| lamboot -d (or lamboot-d and also recon) works. make -k check from the
| lamtest suite passes all tests.
| Is this is
David,
The value of odfWeave is not limited to newbie users. It is vastly
useful for researchers in fields that do not accept LaTeX for journal
paper submission (for example, sociology, demography).
Best,
Shige
On Jan 17, 2008 5:46 PM, David Hajage <[EMAIL PROTECTED]> wrote:
> yes, but we loose
> Thank you. yes it works after loading the package. can you please tell
me
> what should be the location of the file during uploading for aaMI
function.
You can use an absolute path to a file, e.g.
"c:/source/project/myrproject/myfile.r"
.. Or an relative path from the current working director
Good morning,
I use SAS and R/S-Plus as my primary tools so I have a lot of experience with
these programs. By far and away, SAS is superior for handling the "messy"
datasets, but also the very large ones. I work at times with datasets in the
hundreds of thousands (and on occasion, millions)
Thankyou for the replies.
assign() works.
for (i in 1:7)
assign(filesBox[i,1],read.table(paste(dir2, filesBox[i,1], sep=""), header =
FALSE))
--
View this message in context:
http://www.nabble.com/using-an-element-of-an-array-as-a-new-object-tp14884435p14917349.html
Sent from the R help mail
Thank you Professor...
Shubha Karanth | Amba Research
Ph +91 80 3980 8031 | Mob +91 94 4886 4510
Bangalore * Colombo * London * New York * San José * Singapore *
www.ambaresearch.com
-Original Message-
From: Prof Brian Ripley [mailto:[EMAIL PROTECTED]
Sent: Thursday, January 17, 2008 2
Mango Solutions are pleased to announce the above course in London as
part of our schedule for Q1 2008.
---
Introduction to R and R Programming - 6-7th March 2008
--
>If I understand your question:
>
>x <- rnorm(100)
>plot(x)
>legend("topright", capture.output(t(t(summary(x)
Thank you for your help, but I'm afraid that is not what I meant.
As an example of what i am trying to say, imagine I had grouped a whole
bunch of people into 4 age ranges and then
Hi Navish,
did you run
require(aaMI)
?
Cheers
Andrew
On Thu, Jan 17, 2008 at 02:17:12AM -0800, navish wrote:
>
> hi
> i am new to R language. I want to use aaMI package which calculates the
> amino acid mutual interaction for a given protein sequence. I had installed
> the package but when
Dear Richie.
Thank you. yes it works after loading the package. can you please tell me
what should be the location of the file during uploading for aaMI function.
Richard Cotton wrote:
>
>> i am new to R language. I want to use aaMI package which calculates the
>> amino acid mutual interaction
Ricardo Perrone <[EMAIL PROTECTED]> wrote:
>
>Hi,
>
>How to join two large vectors ordered, where one
>has the variable's levels and another has the
>frequencies, in way similar to that showing by table
>function in R console? and considering this two
>vectors how to use summary function to pro
> i am new to R language. I want to use aaMI package which calculates the
> amino acid mutual interaction for a given protein sequence. I had
installed
> the package but when i run the program it gives me the error could not
find
> function "aaMI". can anyone tell me what might be the problem..
hi
i am new to R language. I want to use aaMI package which calculates the
amino acid mutual interaction for a given protein sequence. I had installed
the package but when i run the program it gives me the error could not find
function "aaMI". can anyone tell me what might be the problem..
--
V
yes, but we loose the possibility to use a 'simple' text editor. For
example, I like to use emacs+ess to edit and evaluate R code, and to write
my report in the same editor.
I like the idea that the input format could be writen with a simple text
editor, and the output format be chosen after. I kno
I wonder if those who complain about SAS as a programming environment have
discovered SAS/IML which provides a programming environment akin to Matlab
which is more than capable (at least for those problems which can be treated
with a matrix like approach). As someone who uses both SAS and R - grap
Thanks for your help! Works like a charm now - I can even append to an
expression abject as if it was plain 'c()' ...
Joh
Prof Brian Ripley wrote:
> On Thu, 17 Jan 2008, Johannes Graumann wrote:
>
>> Hi all,
>>
>> I asked something like this earlier but decided that a proper minimal
>> example
On Thu, 17 Jan 2008, Johannes Graumann wrote:
> Hi all,
>
> I asked something like this earlier but decided that a proper minimal
> example might be helpfull ;0)
>
> Why does this work with regards to the expression (substitution):
>
> require(stats)
> plot(cars)
> text(5,120,labels=substitute(i^{
Please re-check your time-series books. The acf at lag 1 is _not_ the
correlation between x and lag(x). For one thing, the variance of x is
computed from the whole series, and not from the series with either the
first or last value removed -- there is also the question of the divisor.
See MAS
Hi all,
I asked something like this earlier but decided that a proper minimal
example might be helpfull ;0)
Why does this work with regards to the expression (substitution):
require(stats)
plot(cars)
text(5,120,labels=substitute(i^{z+phantom()}*"("*a*"
AMU)",list(i="yx",z=2,a=0)))
text(c(5,5),c(
Hi R,
I have doubt.
>x= c(4,5,6,3,2,4,5)
>acf(x,plot=F,lag.max=1)
Autocorrelations of series 'x', by lag
0 1
1.000 0.182
But if I actually calculate the autocorrelation at lag1 I get,
>cor(x[-1],x[-length(x)])
[1] 0.1921538
Even in excel I get 0.1921538 val
73 matches
Mail list logo
