guRus and useRs,
As I instruct my students: "With R what is difficult in anything
else, is easy, usually one line of code; but, what is easy in
anything else, is *&^%$%#$... (translation: next to impossible) in R."
I just ran into such a case that has, to put it mildly, driven me to
use spr
?seq
anil kumar rohilla wrote:
> Hi List,
> I am new to this Rsoftware, i want to make a sereis for example which
> is having values like this, s<- 0,0.1,0.2,0.3,0.4,1
>
> i tryed this statement
> s<-0:0.1:1
> but this giving an error megssage.
> but by default increment 1 it is
On 14 Jun 2006, anil kumar rohilla wrote:
> Hi List,
> I am new to this Rsoftware, i want to make a sereis for example which
> is having values like this, s<- 0,0.1,0.2,0.3,0.4,1
>
> i tryed this statement
> s<-0:0.1:1
> but this giving an error megssage.
> but by default incremen
Hi List,
I am new to this Rsoftware, i want to make a sereis for example which is
having values like this, s<- 0,0.1,0.2,0.3,0.4,1
i tryed this statement
s<-0:0.1:1
but this giving an error megssage.
but by default increment 1 it is taking ,so what to do ,,
i want to use thi
On Tuesday 13 June 2006 13:00, Erin Hodgess wrote:
> Dear R People:
>
> I am using R for teaching purposes in a large classroom.
>
> Each computer has its own copy of R.
>
> However, every once in a while, about half of us will get
> thrown out of R, for no apparent reason.
>
> By the way, it has h
On Wed, 14 Jun 2006, Gad Abraham wrote:
> Prof Brian Ripley wrote:
>> On Tue, 13 Jun 2006, Gad Abraham wrote:
>>
>>> Prof Brian Ripley wrote:
On Tue, 13 Jun 2006, Gad Abraham wrote:
> Hi,
>
> I'm trying to plot several lag.plots on a page, however the second plot
> rep
Prof Brian Ripley wrote:
> On Tue, 13 Jun 2006, Gad Abraham wrote:
>
>> Prof Brian Ripley wrote:
>>> On Tue, 13 Jun 2006, Gad Abraham wrote:
>>>
Hi,
I'm trying to plot several lag.plots on a page, however the second plot
replaces the first one (although it only takes up the upp
The restriction is that objects are kept in memory
so if you have sufficient memory and your OS lets you
access it then you should be ok. S-Plus is a commercial
package similar to R but stores its objects in files and can handle
larger data sets if you run into trouble.
Given that R is free and o
Hi there,
I'm very new to R and am only in the beginning stages of investigating
it for possible use. A document by John Maindonald at the r-project
website entitled "Using R for Data Analysis and Graphics: Introduction,
Code and Commentary" contains the following paragraph, "The R system may
stru
[sorry for my previous empty follow-up -- my fingers got messed up.]
Spencer Graves <[EMAIL PROTECTED]> wrote:
> Did you try the following:
>
>> xto-xfrom
> Time in days:
> [1] 1.50 1.75 2.00 2.25
I wanted to find the total amount of time between a series of date/time's
corresponding to particu
On Tue, 13 Jun 2006 18:49:09 -0700,
Spencer Graves <[EMAIL PROTECTED]> wrote:
> Did you try the following:
>
>> xto-xfrom
> Time in days:
> [1] 1.50 1.75 2.00 2.25
>
Also, have you seen Gabor Grothendieck and Thomas Petzoldt. "R
help desk: Date and time classes in R". R News, 4
Did you try the following:
> xto-xfrom
Time in days:
[1] 1.50 1.75 2.00 2.25
Also, have you seen Gabor Grothendieck and Thomas Petzoldt. "R help
desk: Date and time classes in R". R News, 4(1):29-32, June 2004.,
downloadable from www.r-project.org -> "Documentation: Newsle
Prof Brian Ripley wrote:
> On Tue, 13 Jun 2006, Gad Abraham wrote:
>
>> Prof Brian Ripley wrote:
>>> On Tue, 13 Jun 2006, Gad Abraham wrote:
>>>
Hi,
I'm trying to plot several lag.plots on a page, however the second plot
replaces the first one (although it only takes up the upp
Nathan Johnson <[EMAIL PROTECTED]> wrote:
> I constructed a data frame with columns corresponding to the table
> fields, and tried with and without a null column for the auto-
> increment internal id column. I also tried renaming all the column
> names to match those of the table fields. Al
On 6/13/06, H. Paul Benton <[EMAIL PROTECTED]> wrote:
> Thank you,
> That was the problem. I know about using the extra brackets but in the R
> help manual it doesn't use brackets so I thought maybe R is different.
The problem is not that you must use brackets -- you don't. The problem is that
yo
Thank you,
That was the problem. I know about using the extra brackets but in the R
help manual it doesn't use brackets so I thought maybe R is different.
However, yes indeed, it was assigning the 'lgAB to 4'. I'm just starting to
really learn R.
Cheers,
Paul
PS thanks to the Patrick Burns who
RSiteSearch('Rosner')
?RSiteSearch
or search directly from CRAN.
Incidentally, I'll repeat what I've said before. Don't do outlier tests.
They're dangerous. Use robust methods instead.
-- Bert Gunter
Genentech Non-Clinical Statistics
South San Francisco, CA
"The business of the statistician
Hi,
Put all you .dat files in one directory ('c:/data' for example) and do this:
finaltable <-NULL
folder <-dir('c:/data')
for (i in folder){
table1 <-read.table(i)
mean1 <-mean(table1$V1)
var1 <-var(table1$V1)
rowname <-paste(i)
finaltable <-rbind(finaltable,c(rowname,mean1,var1))
}
Neuro
>Fr
Not sure what your increments are, but something like this might work:
result <- NULL
for (i in 4:5){
x <- read.table(paste('0.', i, '.dat', sep=''))
E <- data$V1[1000:length(data$V1)]
result <- rbind(result, c(paste('0.', i, sep=''), mean(E), var(E)
}
# you can then output 'result' to a
My second request of the day, sorry to be such a bother.
Can you tell me whether Rosner's test for outliers is available in
any of the R packages and, if so, which one? I've tried but I can't
find it.
Thank you very much,
Bob Powell
__
R-help@stat
On 13-Jun-06 Erin Hodgess wrote:
> Dear R People:
>
> I am using R for teaching purposes in a large classroom.
>
> Each computer has its own copy of R.
>
> However, every once in a while, about half of us will get
> thrown out of R, for no apparent reason.
>
> By the way, it has happened in ot
On 6/13/06, Rolf Turner <[EMAIL PROTECTED]> wrote:
> H. Paul Benton wrote:
>
> > Umm sorry to bother everyone again but I'm having trouble with my if
> > statement. I come from a perl background so that's probably my problem! :)
> > So here is my code:
> >
> > if (any(lgAB>4) | any(lgAB<-4)){
> >
H. Paul Benton wrote:
> Umm sorry to bother everyone again but I'm having trouble with my if
> statement. I come from a perl background so that's probably my problem! :)
> So here is my code:
>
> if (any(lgAB>4) | any(lgAB<-4)){
> freq_AB<-hist(lgAB, type="o", plot=F)
> else
> f
I am trying to use the ncdf package on a Mac OS X machine (10.4.6), but
I get the following error message:
> library(ncdf)
Error in dyn.load(x, as.logical(local), as.logical(now)) :
unable to load shared library
'/Library/Frameworks/R.framework/Resources/library/ncdf/libs/ppc/ncdf.so
':
Dear R People:
I am using R for teaching purposes in a large classroom.
Each computer has its own copy of R.
However, every once in a while, about half of us will get
thrown out of R, for no apparent reason.
By the way, it has happened in other classrooms as well.
Has anyone else run into th
I have many files (0.4.dat, 0.5.dat, ...) of which I would like to calculate
mean value and variance and save the output in a new file where each line
shouldlook like: "0.4 mean(0.4.dat) var(0.4.dat)" and so on. Right now I got a
a simple script that makes me unhappy:
1. I run it by "R --no-sav
... and again I wonder which courses are "near". This leads at once to
the question: "which metric is in use?". Probably some football related
metric: FIFA WM takes place in Dortmund and commercials say something
like "the world is our guest" ...
Now, let's escape from football to Austria and Vi
Hi, this is my first help request so please bear with me.
I've been running some normality tests using the nortest package. For
some of my datasets the Cramer-von Mises normality test generates an
extremely high probability (e.g., 1.637e+31) and indicates normality
when the other tests do no
I am using R 2.3.0 and trying to run the following
code.
BinCnt=round(ShorterVectorLength/10)
BinMatrix=matrix(nrow=BinCnt, ncol=2)
Increment=(VectorRange[2]-VectorRange[1])/BinCnt
BinMatrix[1,1]=VectorRange[1]
BinMatrix[1,2]=VectorRange[1]+Increment
Sarah Goslee wrote:
> You need more brackets:
>
> if(blah) {
> do something
> } else {
> do something different
> }
>
> Sarah
>
> PS Using underscores in variable names is not encouraged, and
> can cause you problems in certain contexts.
>
Hi, Sarah,
Why do you say this? And in what sit
Since several people have commented to me, on and off list:
On 6/13/06, Sarah Goslee <[EMAIL PROTECTED]> wrote:
> >
> > PS Using underscores in variable names is not encouraged, and
> > can cause you problems in certain contexts.
I realize that it is now legal in R to use underscores, since they
On 6/13/06, Sarah Goslee <[EMAIL PROTECTED]> wrote:
> You need more brackets:
>
> if(blah) {
> do something
> } else {
> do something different
> }
>
> Sarah
>
> PS Using underscores in variable names is not encouraged, and
> can cause you problems in certain contexts.
The only context I can thi
replace 'else' with '} else {'
Gabor
On Tue, Jun 13, 2006 at 11:18:00AM -0700, H. Paul Benton wrote:
> Umm sorry to bother everyone again but I'm having trouble with my if
> statement. I come from a perl background so that's probably my problem! :)
> So here is my code:
>
> if (any(lgAB>4) | any
You need more brackets:
if(blah) {
do something
} else {
do something different
}
Sarah
PS Using underscores in variable names is not encouraged, and
can cause you problems in certain contexts.
On 6/13/06, H. Paul Benton <[EMAIL PROTECTED]> wrote:
>
> Umm sorry to bother everyone again but
Umm sorry to bother everyone again but I'm having trouble with my if
statement. I come from a perl background so that's probably my problem! :)
So here is my code:
if (any(lgAB>4) | any(lgAB<-4)){
freq_AB<-hist(lgAB, type="o", plot=F)
else
freq_AB<-hist(lgAB, breaks=br,type
XLSolutions Corporation (www.xlsolutions-corp.com) is proud to announce:
(1) Introduction to R/S+ programming: Microarrays Analysis and Bioconductor
*** San Francisco / July 17-18, 2006 ***
*** Chicago / July 24-25, 2006 ***
*** Baltim
Hello,
I am a new user and I am looking for the description of the output of
rqss function (Additive Quantile Regression Smoothing). It is supposed to
be in rqss.object but I could not find any reference to rqss.object
anywhere.
thanks a lot.
Julia
[[alternative HTML version deleted
On Tue, 13 Jun 2006, Nathan Johnson wrote:
> I'm new to R and have been resisting posting to the list thus far,
> but I think I've exhausted my work around options.
Have you asked the package maintainer: the posting guide suggests you do
so first? He is by far the most likely source of help.
>
Hi
I'm new to R and have been resisting posting to the list thus far,
but I think I've exhausted my work around options.
I'm having a few problems with RMySQL and/or it's underlying packages.
I'm doing a very simple vsn transformation on large datasets, using
RMySQL to retrieve the data for
This page:
http://depts.washington.edu/dsc2005/
says the proceedings for DSC 2005 will be published online and that
the papers were due in final form by September 2005.
Anyone know if the proceedings have been published yet or when (if?)
they will be published?
Thanks.
Kevin Wright
Hello,
I have a data frame :
> linkptpn[1:6,]
pedigree id fa mo sex status rs3789604 rs3811021 rs1217413 ss38346942
1 1 1 0 0 1 2
2 1 2 0 0 2 2
3 1 3 1 2 2 3 3/3 2/2 3/34/4
4 1 4 1 2 2 3 3/3
On Tue, 13 Jun 2006, Niels Vestergaard Jensen wrote:
> I have been using levelplot but have had trouble calling it inside
> functions - something seems to go wrong when it's not called directly from
> the R command prompt. Simplest reproducible example:
>
> $ R --vanilla
>> library(lattice)
>> lev
This is 7.22 of the R FAQ:
http://cran.r-project.org/doc/FAQ/R-FAQ.html#Why-do-lattice_002ftrellis-graphics-not-work_003f
On 6/13/06, Niels Vestergaard Jensen <[EMAIL PROTECTED]> wrote:
> I have been using levelplot but have had trouble calling it inside
> functions - something seems to go wrong
I have been using levelplot but have had trouble calling it inside
functions - something seems to go wrong when it's not called directly from
the R command prompt. Simplest reproducible example:
$ R --vanilla
> library(lattice)
> levelplot(matrix(1:4,2,2))
- This gives a nice plot in soothing pas
Albert Vilella wrote:
> Hi all,
>
> I would like to know if it is possible to sort the columns and rows in
> a matrix given a specified order of their colnames and rownames. For
> example:
>
> I would like the original matrix:
>
> RNO FRU ANG CEL DAR PTR
> RNO 3.45 1.35 2.16 2.25 1
Dear R users,
Suppose that I generate the following :
library(survival)
set.seed(123)
n=200
x=rbinom(n,size=1,prob=0.5)
v=rep(rgamma(n/2,shape=1,scale=1),2)
w=rweibull(n,shape=1,scale=1)
b=-log(2)
t=exp(-x*b-log(v)+log(w))
c=rep(1,n)
id=rep(seq(1:(n/2)),2)
group=rep(1:(n/10),10)
fit2=coxp
Hi all,
I would like to know if it is possible to sort the columns and rows in
a matrix given a specified order of their colnames and rownames. For
example:
I would like the original matrix:
RNO FRU ANG CEL DAR PTR
RNO 3.45 1.35 2.16 2.25 1.43 1.20
FRU 1.31 Inf 2.22 2.36 1.34 1.4
Just a quick point which may be easy to correct. Whilst typing the
wrong thing into R 2.2.1, I noticed the following error messages,
which seem to have some stray quotation marks and commas in the list
of available families. Perhaps they have been corrected in the latest
version (sorry, I d
Hi all,
I would like to ask if it is possible to use bubbleplot for a 20x20
matrix, instead of a dataframe with factors in columns.
The idea would be to get a tabular representation with bubbles like in
Rnews_2006_2 article, which look very nice.
Thanks in advance,
Albert.
Dear all R-users,
I wanted to fit a Garch(1,1) model to a dataset by:
>garch1 = garch(na.omit(dat))
But I got a warning message while executing, which is:
>Warning message:
>NaNs produced in: sqrt(pred$e)
The garch parameters that I got are:
> garch1
Call:
garch(x = na.omit(dat))
Coefficie
On Tue, 13 Jun 2006, Jouanin Celine wrote:
I want to do a nonparametric regression. Im using the function loess.
The variable are the year from 1968 to 1977 and the dependant variable is
a proportion P. The dependant variable have missing value (NA).
The script is :
year
curve(coef$A * dnorm(x, coef$mu, coef$sig), -4, 4)
---
Jacques VESLOT
CNRS UMR 8090
I.B.L (2ème étage)
1 rue du Professeur Calmette
B.P. 245
59019 Lille Cedex
Tel : 33 (0)3.20.87.10.44
Fax : 33 (0)3.20.87.10.31
http://www-good.ibl.
I want to do a nonparametric regression. Im using the function loess.
The variable are the year from 1968 to 1977 and the dependant variable
is a proportion P. The dependant variable have missing value (NA).
The script is :
year <- 1969:2002
length(year)
Eric Hu wrote:
> Hi I am trying to plot two data set in the same picture window without
> overlapping with each other. I am using the format plot(x1,y1,x2,y2)
> but get the following error message:
>
>> plot(as.numeric(r0[,4]),as.numeric(r0[,7]),as.numeric(r0[,4]),as.numeric(r0[,7][ind[,1]]))
> Er
Hello
thanks for the point, probably i saw DBMS as a need due to my ignorance about R.
If i can process all files like you said, i would not use DBMS as i
prefer to keep it simple and easy to manage and run, less software
dependencies the better.
thanks
gabi
On 6/12/06, bogdan romocea <[EMAIL
On Tue, 13 Jun 2006, Gad Abraham wrote:
> Prof Brian Ripley wrote:
>> On Tue, 13 Jun 2006, Gad Abraham wrote:
>>
>>> Hi,
>>>
>>> I'm trying to plot several lag.plots on a page, however the second plot
>>> replaces the first one (although it only takes up the upper half as it
>>> should):
>>>
>>
Prof Brian Ripley wrote:
> On Tue, 13 Jun 2006, Gad Abraham wrote:
>
>> Hi,
>>
>> I'm trying to plot several lag.plots on a page, however the second plot
>> replaces the first one (although it only takes up the upper half as it
>> should):
>>
>> par(mfrow=c(2,1))
>> a<-sin(1:100)
>> b<-cos(1:100)
On Tue, 13 Jun 2006, Gad Abraham wrote:
> Hi,
>
> I'm trying to plot several lag.plots on a page, however the second plot
> replaces the first one (although it only takes up the upper half as it
> should):
>
> par(mfrow=c(2,1))
> a<-sin(1:100)
> b<-cos(1:100)
> lag.plot(a)
> lag.plot(b)
>
> What's
58 matches
Mail list logo
