François Pinard wrote:
> I wonder what happened, for R to hide the underlying Scheme so fully, at
> least at the level of the surface language (despite there are hints).
As far as I understood, the original version of Ihaka/Gentleman was
written in Scheme. But even if you look at the source ve
Hi
Jake Verschuyl wrote:
> Hi there,
>
>
>
> I have some bird flight height data that follows a gamma distribution. The
> data (x) goes from 0 to 700 meters (n=1055). The calculated parameters
> calculated from the fitdistr(x) are (shape = 5.1379, rate = 0.017541), and
> therefore the scale
I hope you don't mind that I offer also two solutions. No.1 is really
bad. No.2 should be on par with the other ones.
Best,
Roland
mydata <- matrix(rnorm(10*10), ncol=10)
threshold.value <- 1.5
mydata2 <- matrix(0, nrow=nrow(mydata), ncol=ncol(mydata))
mydata3 <- matrix(0, nrow=nrow(mydata),
Alberto Monteiro wrote:
> I have a time series x = f(t), where t is taken for each
> month. What is the best function to detect if _x_ has a seasonal
> variation? If there is such seasonal effect, what is the
> best function to estimate it?
>
From my own experience, I had the impression that ther
Hi,
Tiandao Li wrote:
> Hello,
>
> I have some files generated from microarray experiments. I used scan() to
> read the files, and assigned each file to a unique name with many rows and
> columns. Now I want to create a array (ArrayA) with unique names, and I can
> use ArrayA[1,2][[6]] to refer t
Dong GUO 郭东 wrote:
> Hi,
>
> I want to save a array (say, array[6,7,8]) write a cvs file. How can I do
> that??? can I write in one file?
For array[6,7,8], you don't need a csv(!) file since it is only a
scalar. If this is what you want, check
?write.table
But what you probably meant is how to w
Eric Doviak wrote:
>
> I need to find some way to overcome these constraints and work with large
> datasets. Does anyone have any suggestions?
I might be not the most authoritative person on this subject but I put
all my large datasets[1] into an SQLite database and extract/summarize
data from
Dear all,
I saw that the Free Software Foundation (www.fsf.org) released Version 3
of their General Public Licence on Friday.
Just out of curiosity, are there any plans to 'upgrade' the licensing
terms of R?
Maybe one can update the licence() function since the link given there
(http://www.gnu
Dear all,
I just saw for the first time the 'Matlab Programming Contest' at
http://www.mathworks.com/contest/overview.html
Is there something comparable for R? Is someone interested in organizing
such a thing?
Of course, one can try to solve the Matlab Problems Sets in R, but I
think it could
Hi Eric,
eric zoukekang wrote:
> Hello!
> I wonder if you might help me with informations about how to import data with
> a 2.4.1 R version without the menu "Import data".
> Best regards.
the best information about importing data in R is contained in the
manual "R Data Import/Export" which shou
Hi K,
Devarajan, Karthik wrote:
> Does anyone have any suggestions\thoughts
Hard to say where the problem(s) originate(s). Why don't you provide
commented, minimal, self-contained, reproducible code as is recommended
on all messages on this mailing list?
This would make it much easier to trac
Hi,
since you are very new to R, just a small advice: try to give a minimal
reproducible, self-contained example. This helps quite often to spot
some mistakes yourself. This is my own experience and has been
experienced by others, see, for example:
library(fortunes)
fortune("Ripleyed")
:-)
[EMAIL PROTECTED] wrote:
> Hi There,
>
> How to display a number in the formation of percentage?
>
> For example, a=0.25, how to display 25%?
if you write 'how to display' I guess you mean in a plot?
a <- 0.25
plot(0:10, 0:10, type="n")
text(x=5, y=5, labels=paste(round(a*100,0), "%", sep=""))
Hi Matthias,
please see below for some code example.
Matthias von Rad wrote:
> Hi,
> I am doing lots of scatterplots for my dissertation and to make the
> comparable, I would like to have equal x- amd - y axis. Can I specify
> their scale?
> Another question adresses overlay scatterplots. Having
Hi Neil,
[EMAIL PROTECTED] wrote:
>
> At this point, want to jump in avoiding all the
> Mathematical proofs and just apply R and the packages for what I want to
> do.
>
I'd still recommend Venables/Ripley: Modern Applied Statistics with S
(or often abbrev. MASS, which is also name of the packag
Hi,
[EMAIL PROTECTED] wrote:
> I am new to using R and would appreciate some advice on
> which books to start with to get up to speed on using R.
>
> My Background:
> 1-C# programmer.
> 2-Programmed directly using IMSL (Now Visual Numerics).
> 3- Used in past SPSS and Statistica.
>
> I put toget
李俊杰 wrote:
> Dear R-lister,
>
> One of my friends wanted to produce random number which is 64 bits. He did
> it with Fortune. I think R can do it also. But I don't know how to display a
> very big integer in the complete form but not scientific form. And what's
> the biggest integer R can display
Hi Jonathan,
Jonathan Morse wrote:
> I am new to Linux (not to R) and recently installed Mandriva Spring 2007 on
> my partitioned hard drive. My next objective is to install R in the Linux
> environment, unfortunately Mandriva is not one of the Linux distributions
> available for download...
Hi Patrick,
[EMAIL PROTECTED] wrote:
> Dear R-help list members,
>
> I am new to R, and having problems with plotting part of a time series. I
> have read in my data using read.table, and my 'time' column is recognised
> as a numeric variable. When I convert this to a time format, I am no longe
Hi,
assuming you are doing a Cox-PH-Model, you can check:
library(survival)
?coxph.object
There it says that the components 'residuals' refers to the martingale
residuals. I hope I recall it correctly but there exists a simple
relation between the martingale residuals ('mgr') and the unmodified
[EMAIL PROTECTED] wrote:
> Hi,
>
> Please see the email exchanges below. I am having trouble generating
output that is large enough
> for our needs, specifically when using the GaussRF function.
However, when I wrote Dr. Schlather
> (the author of the GaussRF function), he indicated that th
Markus voigt wrote:
> because it should start at 2. Is there a parameter where you can define
> the start on the x-axis?
Maybe this will help you?
datax1 <- 1:8
datax2 <- 2:9
datay1 <- runif(length(datax1))
datay2 <- runif(length(datax2))
plot(x=datax1, y=datay1, type="l", col="blue", xlab="X",
ss. So the '@'-character has some
ambiguous meaning in my setting. Sorry that I did not think of that
beforehand.
Thank you very much for your help,
Roland
Prof Brian Ripley wrote:
> On Wed, 16 May 2007, Roland Rau wrote:
>
>> Dear all,
>>
>> in the past I have bee
Dear all,
in the past I have been able to access websites with data directly. For
example the following code works nicely
mydata <-
read.table("http://www.lifetable.de/data/MPIDR/POL_2004.txt";, header=TRUE)
But what happens if I need a username and password (a different site)?
How do I do
mister_bluesman wrote:
> Thanks people! Sorry to sound thick but after I download the poltrix package,
> how do i install it?
>
> Thank you!
which platform are you on?
On Windows I would suggest you use the Menu of the GUI "Install packages
from local zip files" in the main category 'packages'
Hi,
summary(yourdata) is often a very good starting point for descriptive
data statistics. Or you can write your own little function which returns
what you actually like to see (the code below was written very quickly.
No care is taken for the presence of missing values or anything else).
exa
Well, I guess it makes quite a difference in survival analysis whether
you know that a person was alive/censored or experienced the event of
interest at a certain point of time/age. You could have tried it easily
for yourself by slightly modifying the example on the help page of
'survdiff'.
li
Hi Jessica
[EMAIL PROTECTED] wrote:
> __
> R-help@stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-co
!
Roland
Roland Rau wrote:
> Dear all,
>
> I have an ASCII file where records are separated by a blank. I would
> like to read those data; however, only the data in rows 1, 3, 5, 7, ...
> are important; the other lines (2,4,6,8,) contain no useful
> information for me.
&
And since Peter Dalgaard also just answered (without advertising his book):
if you (or your library) happen to have 'Introductory Statistics with R' by
Peter Dalgaard, have a look at section 12.4.
Best,
Roland
On 5/1/07, raymond chiruka <[EMAIL PROTECTED]> wrote:
>
> how do l programme the logran
Hi,
On 5/1/07, raymond chiruka <[EMAIL PROTECTED]> wrote:
>
> how do l programme the logrank test. l am trying to compare 2 survival
> curves
if you simply want to use the logrank test, have a look at the first example
of the function survdiff in the survival package. If you read the help page
t
Dear all,
I have an ASCII file where records are separated by a blank. I would like to
read those data; however, only the data in rows 1, 3, 5, 7, ... are
important; the other lines (2,4,6,8,) contain no useful information for
me.
So far I used awk/gawk to do it:
gawk '{if ((FNR % 2) != 0) {pr
On 4/17/07, Luca Braglia <[EMAIL PROTECTED]> wrote:
>
> I have found a strange "ifelse" behaviour (I think)
Don't you think it is rather consistent behavior?
> ifelse(T,1+1,1+2)
[1] 2
> ifelse(F,1+1,1+2)
[1] 3
> ifelse(T,"hello","goodbye")
[1] "hello"
> ifelse(F,"hello","goodbye")
[1] "goodbye"
Hi,
wouldn't it be enough for simple purposes (as stated in the original
message: "Greatest common divisor of two numbers") to use:
gcd <- function(a,b) ifelse (b==0, a, gcd(b, a %% b))
> gcd(12,4)
[1] 4
> gcd(4,12)
[1] 4
> gcd(123456789,987654321)
[1] 9
> gcd(987654321,123456789)
[1] 9
>
Best,
I am not sure what the question is but maybe this provides a starting point?
http://cran.r-project.org/bin/linux/ubuntu/README
HTH,
Roland
On 4/16/07, Erin Hodgess <[EMAIL PROTECTED]> wrote:
>
>
> __
> R-help@stat.math.ethz.ch mailing list
> https://st
Dear all,
on programming.reddit.com, I found a link to a PDF document which compares
on 17 pages commands/functions in Matlab, Python and R.
http://37mm.no/download/matlab-python-xref.pdf
I hope I am not re-posting which has been posted recently.
Best,
Roland
[[alternative HTML version
Hi Lorenzo,
On 4/5/07, Lorenzo Isella <[EMAIL PROTECTED]> wrote:
>
> I would like to have suggestions about where to collect info about:
> (1)Institutions (not only academia) using R
A starting point might be to look at the R-project homepage and look at the
members and donors list. This is, of
Hi,
On 3/7/07, j.joshua thomas <[EMAIL PROTECTED]> wrote:
>
> I have used RODBC to get the database i can view the tables in RGUI-2.4.1
> how can i query the records in R>
> i tried with sqlQuery
>
> need some help
>
> A bit of code showing what you did would help others to track down what
may hav
Hi,
On 2/21/07, Venkatesh Mantha <[EMAIL PROTECTED]> wrote:
>
> We are working on a project on forecast modeling and would like to know if
> there are any examples on how to READ/WRITE to a database (e.g.
> PostgreSQL)
did you have a look at the Manual "R Data Import/Export"? I have to admit
th
On 2/16/07, Petr Pikal <[EMAIL PROTECTED]> wrote:
>
> Hi
>
> slight modification of your function can be probably even quicker:
>
> fff<-function(x) sum(2^(which(rev(x))-1))
> :-)
> Petr
>
>
Yes, your function is slightly but consistently faster than my suggestion.
But my "tests" show still Bert Gu
not simply:
>
> sum(x * 2^(rev(seq_along(x)) - 1)) ?
>
>
> Bert Gunter
> Genentech Nonclinical Statistics
> South San Francisco, CA 94404
> 650-467-7374
>
>
> Bert Gunter
> Nonclinical Statistics
> 7-7374
>
> -Original Message-
> From: [EMAI
That was a nice quick distraction. Unfortunately, I am not the first to
answer. :-(
Anyway, I offer two solutions (which are different from the one of Marc
Schwartz); I wrote it quickly but I hope they are correct.
Enjoy and thanks,
Roland
a <- c(TRUE, FALSE, TRUE)
b <- c(TRUE, FALSE, TRUE, TRUE)
Thanks Duncan and Greg.
My current solution is to use the rgl-package.
Is there an easy way to obtain a screenshot in eps- or pdf-Format from such
an rgl-window?
I saw the rgl.snapshot function but it does not provide this format.
So far, I take a snapshot, save it as jpeg and convert it to eps vi
Dear all,
I would like to put two perspective plots into one plot. The help page
for ?persp shows how one can add points and lines but not another
perspective plot.
data(volcano)
z <- 2 * volcano# Exaggerate the relief
x <- 10 * (1:nrow(z)) # 10 meter spacing (S to N)
y <- 10 * (1:ncol(
> someone to start doing these things! It was written for Linux, but the
> only difference on Windows is that you need those tools I mentioned above!
>
> Oleg
>
> --
> Dr Oleg Sklyar | EBI-EMBL, Cambridge CB10 1SD, UK | +44-1223-494466
>
>
> Roland Rau wro
Hi,
On 2/8/07, andy1983 <[EMAIL PROTECTED]> wrote:
>
>
> That was a neat trick. However, it created a new problem.
>
> Before, it took way too long for a 10,000 columns to finish.
>
> Now, I test the memory limit. With 10,000 columns, I use up about 1.5 GBs.
>
> Assuming memory is not the issue, I
Hi,
On 2/8/07, Oleg Sklyar <[EMAIL PROTECTED]> wrote:
>
> On Windows you need:
> - download and install Cygwin (cygwin.com) with default options,
> supposedly you install into c:\cygwin. Add path to
> c:\cygwin\bin;c:\cygwin\lib to your system PATH
No. You don't need Cygwin. I don't ha
Hi,
On 2/8/07, Albrecht, Dr. Stefan (AZ Private Equity Partner) <
[EMAIL PROTECTED]> wrote:
>
>
> I would very much appreciate any comments on my above remarks. I know
> there has been some discussions of R vs. Matlab on R-help, but these
> could be somewhat out-dated, since both languages are evo
Hi,
On 2/2/07, Martin Olivier <[EMAIL PROTECTED]> wrote:
>
> Hi all,
>
> Suppose I have a vector x with numerical values.
> In y, I have a categorial variable : y<-c(1,1,..2,2,...30,30,30)
> x and y have the same length.
> I would like to compute the mean for x for the modality 1 to 30 in y.
> me
On 2/1/07, xiaopeng hu <[EMAIL PROTECTED]> wrote:
>
> i get the message :
> configure: WARNING: you cannot build info or html versions of the R
> manuals
>
> how to deal with it ?
did you read the manual 'R Installation and Administration'?
see section 2.2 there.
it is mentioned there that you ne
This is what I was hoping for!
Thanks,
Roland
On 1/29/07, Achim Zeileis <[EMAIL PROTECTED]> wrote:
>
>
>
> On Mon, 29 Jan 2007, Roland Rau wrote:
>
> > Dear all,
> >
> > given I have a data.frame in a format like this
> >
> > mydf <- data.f
Dear all,
given I have a data.frame in a format like this
mydf <- data.frame(age=rep(1:3,5),
year=c(rep(1996,3), rep(1997,3), rep(1998,3),
rep(1999,3), rep(2000,3)),
income=1:15)
mydf
Now I convert it to some 2D-frequency table like thi
Hi,
David Mas wrote:
> I believe that the way basehaz (in the survival package) compute the
> baseline hazard function is false.
>
> I come to question this function when it gives me hazard probabilities
> greater than 1.
>
> Looking at the code I think I've localised the error :
>
> hazard proba
53 matches
Mail list logo
