hi Srinivas
I'm not sure I understand: your example result has 10-9 for mu1 but 1
for mu2, which is not equal to 11-9. I will assume you just got the
numbers wrong. (Also J100 should be J1, no?)
I am also wondering why you want what you say. It doesn't seem very
sensible to mix the results from t
Asserting `y isn't noramal' only eliminates one possibility. Is it
binary, a count, long-tailed, short-tailed, ...?
This is really a statistical consulting problem (and apparently about
SAS): see the R posting guide for where to seek advice on such problems.
As the guide says, had we had suffic
In non-parametric bootstrapping, you do draw a sample of the same size as
the original. So
D <- as.vector(WorldPhones[,1])
sample(D, replace=TRUE)
[1] 45939 64721 60423 79831 71799 45939 76036
sample(D, replace=TRUE)
[1] 45939 68484 79831 68484 76036 71799 79831
are two different bootstrap
Dear friends,
I have a dataset: response var--y, class var-group, and the third variable-x.
I want to test whether there is statistical significance bewteen
group for y with the controlled x. First, i want to use analysis of
covariance in SAS, but i found that y isn't noramal and can't become
nor
Have you tried RSiteSearch("spatial ecology")? I just got 47 hits
from that. Some of them might be relevant to your question.
If that fails, you might consider providing this group with the math
behind the automata models you are considering. I might expect them to
be exp
Hello,
I have a simple function in the form as follows:
> f<-function(x){sum(v^x)}
where v is a vector. I was trying to integrate f using the command
> I<-integrate(f,0,1)
However, this will not work and seems that the reason is to use
"integrate", the f must be a function that with input and
Try this which reads in the data as a data frame and appends a
chron variable datetime:
Lines <- " MM DD HH data1 data2
2002 12 01 01 0.002 0.003"
DF <- read.table(textConnection(Lines), header = TRUE)
library(chron)
DF$datetime <- with(DF, chron(paste(MM, DD, , sep = "/")) + HH/24)
DF
2006/6/7, Recep Aykaya <[EMAIL PROTECTED]>:
>
> hi.
> i'm a statistics student and studying bootstrap in R.
>
> i'm trying to draw bootstrap samples from a sample, using the following R
> code:
>
> > *boot = numeric(200)*
> *> {for (i in 1:200)*
>
> * boot [i] = (sample(data,replace=T))}*
>
>
>
>
R-help doesn't foward attached data files like this, but Brian
kindly forwarded it to me.
You need to restrict X so that it is full rank, it now has
rank 19 and column dimension 29 (with intercept). See
for example svd(cbind(1,x)).
I'll add some better checking for this, but it will basically
hi.
i'm a statistics student and studying bootstrap in R.
i'm trying to draw bootstrap samples from a sample, using the following R
code:
> *boot = numeric(200)*
*> {for (i in 1:200)*
* boot [i] = (sample(data,replace=T))}*
i obtain 200 samples but all of them are the same.
i want to obtain
I suggest you investigate with().
On Wed, 7 Jun 2006, Manuel Morales wrote:
> Dear list members,
>
> In some of my functions, I attach the data internally to allow subset
> commands or to specify a data frame. This works well except for cases
> where there is a "masking" conflict (which returns a
On Wed, 7 Jun 2006, roger koenker wrote:
> Since the "crash" occurs plotting the lm object it is unclear what
> this has to do with quantreg, but maybe you could explain
>
> 1. what you mean by crash,
> 2. something about x,y,
>
> This is best addressed to the maintainer of the packa
Dear list members,
In some of my functions, I attach the data internally to allow subset
commands or to specify a data frame. This works well except for cases
where there is a "masking" conflict (which returns a warning). I see
some alternative listed in ?attach, but I'm not sure which of them do
I attached the data file here. I restarted the PC but it still happens. It
says a memory address could not be written. I am not sure it is a problem of
R or quantreg but I plot without problems before I load quantreg.
Thank you.
Tian
On 6/7/06, Prof Brian Ripley <[EMAIL PROTECTED]> wrote:
Wit
Without y and x we cannot reproduce this.
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
On Wed, 7 Jun 2006, Mu Tian wrote:
> I forgot to mention my R version is 2.3.1 and quantreg is the most updated
> too.
It has a version number, which the posting guide tells y
Since the "crash" occurs plotting the lm object it is unclear what
this has to do with quantreg, but maybe you could explain
1. what you mean by crash,
2. something about x,y,
This is best addressed to the maintainer of the package rather than to
R-help, provided, of cou
I forgot to mention my R version is 2.3.1 and quantreg is the most updated
too.
On 6/7/06, Mu Tian <[EMAIL PROTECTED]> wrote:
>
> I was trying "quantreg" package,
>
> lm1 <- lm(y~x)
> rq1 <- rq(y~x)
> plot(summary(rq1)) #then got a warning says singular value, etc. but this
> line can be omited
>
I was trying "quantreg" package,
lm1 <- lm(y~x)
rq1 <- rq(y~x)
plot(summary(rq1)) #then got a warning says singular value, etc. but this
line can be omited
plot(lm1) #crash here
It happened every time on my PC, Windows XP Pro Serv. Pack 1, Pentium(4)
3.00G.
[[alternative HTML version del
Hello much more experienced R-users,
I have got a txt - file which contains data formated like this:
MM DD HH data1 data2 ... data31
2002 12 01 01 0.002 0.003 ... 312.0
The single columns are divided by at least one space.
Is their an easy and fast way to make R understand that the first
I am not sure I understand what you want to do, but maybe some of this
will be helpful. I first generate some data that should resemble yours:
dat<-expand.grid(Region=1:3, Species=1:4, Sex=c("M","F"))
dat<-do.call("rbind",lapply(1:10,function(x) dat))
dat$Bodysize<-rnorm(nrow(dat),10,2)
Now what
Hello all,
Here is the code that I am using for finding differentially expressed genes.
#Normalization
library(affy)
library(Biobase)
library(limma)
library(gcrma)
pd<-read.phenoData("file.txt",header=TRUE,row.names=1,as.is=TRUE,sep="\t")
Data <- ReadAffy(filenames=pData(pd)
[Pikounis, Bill [CNTUS]]
> view <- function(x) {
> warnopt <- options()$warn
> options(warn=-1)
> on.exit({sink(); options(warn=warnopt)})
> edit(x)
> invisible()
> }
I'm surprised by the necessity of "sink()". Presuming it is necessary
indeed, the above could be simplified a bit like
Nicolas Degallier wrote:
> Hi!
>
> I am trying to install in my R environment the rhdf5 package and
> library but it seems to have vanished from either the CRAN or
> BioConductors sites.
>
> Can you tell me where it would be possible to find it or any R
> library (or function) able to read
hello:
I have 3 data.frame objects.
First df object:
Of dim (149,31). Columns 2:31 are marked as T1..T14
and N1..N16.
Name T1T2N1 T3 N2 N3 N4 T4
mu1 1010910 9 9 8 10
mu2 1111911 9 9 9 11
...
muN 1212911 9
>Not a direct answer to your question, but if you use a logspline density
>estimate rather than a kernal density estimate then the logspline
>package will help you and it has built in functions for dlogspline,
>qlogspline, and plogspline that do the integrals for you.
>
>If you want to stick with t
Pedro wrote:
> I have made a simple kernel density estimation by
>
> x <- c(2,1,3,2,3,0,4,5,10,11,12,11,10)
> kde <- density(x,n=100)
>
> Now I would like to know the estimated probability that a
> new observation falls into the interval 0
> How can I integrate over the corresponding interval?
Spencer,
I agree that there exist biases within large organizations toward
commercial software packages, but I humbly disagree that this can be
reduced to geopolitical differences. I work for a US Government Agency
and while there is some truth to the feeling "[they] refuse to believe
that any
Not a direct answer to your question, but if you use a logspline density
estimate rather than a kernal density estimate then the logspline
package will help you and it has built in functions for dlogspline,
qlogspline, and plogspline that do the integrals for you.
If you want to stick with the KDE
Dear R-list,
I have made a simple kernel density estimation by
x <- c(2,1,3,2,3,0,4,5,10,11,12,11,10)
kde <- density(x,n=100)
Now I would like to know the estimated probability that a
new observation falls into the interval 0https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting
Also Consider Bonferroni Hochberg Holm type procedures or .
Dunn OJ. Multiple contrasts using rank sum tests. Technometrics
1964;6:241#/52.
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mai
On Thu, 2006-06-08 at 00:10 +0800, zhijie zhang wrote:
> Dear Rusers,
> As we all know , there are many methods to do multiple comparison in the
> parametric statistical analysis, But i can't find some in nonparametric
> statistical analysis.
> Could anybody give some suggestions?
Have you looke
Perfect. Thank-you!
"Deepayan Sarkar" <[EMAIL PROTECTED]>
06/07/2006 11:19 AM
To
"[EMAIL PROTECTED]" <[EMAIL PROTECTED]>
cc
r-help@stat.math.ethz.ch
Subject
Re: [R] Fw: Help needed using lattice for area plots lpolygon, xyplot.
On 6/7/06, [EMAIL PROTECTED]
<[EMAIL PROTECTED]> wrote:
On 6/7/06, [EMAIL PROTECTED]
<[EMAIL PROTECTED]> wrote:
> I am trying to learn how to use the graphics from the lattice package (
> and am very new to R).
>
> I am trying to replicate the example plot referenced below, by using the
> lattice xyplot & lpolygon to create panels. I get what appears t
Dear Rusers,
As we all know , there are many methods to do multiple comparison in the
parametric statistical analysis, But i can't find some in nonparametric
statistical analysis.
Could anybody give some suggestions?
[[alternative HTML version deleted]]
_
On Wed, 2006-06-07 at 10:30 -0500, Mark L Sessing wrote:
> Hello,
>
> I am learning how to use R, and I cannot figure out how to plot more
> than one data set on a single plot. Can you help me out?
>
> Cheers,
> Mark
It depends upon the type of plot (scatter, lines, bar, etc.) and whether
or n
I use a simple-minded, dirty wrapper on edit:
view <- function(x) {
warnopt <- options()$warn
options(warn=-1)
on.exit({sink(); options(warn=warnopt)})
edit(x)
invisible()
}
I say "dirty", because sometimes the output is re-directed from stdout
somewhere else as an unintended side effec
Without more inormation about what you are specifically trying
to do, see
?line
?points
Sarah
On 6/7/06, Mark L Sessing <[EMAIL PROTECTED]> wrote:
>
> Hello,
>
> I am learning how to use R, and I cannot figure out how to plot more
> than one data set on a single plot. Can you help me out?
>
> Ch
The best rebuttal I've heard recently to arguments like that is Linux
(www.linux.org): It's distributed under the same general public license
(GNU) license as R.
A perspective that I don't recall having seen on this list is that
the cost of producing and distributing softwa
Hello,
I am learning how to use R, and I cannot figure out how to plot more
than one data set on a single plot. Can you help me out?
Cheers,
Mark
--
Mark Sessing
CIMMS Research Fellow Meteorologist
NWS Warning Decision Training Branch
3200 Marshall Ave Ste. 202
Norman, OK 73072
Phone: 405-57
Hi Mike,
If you're not committed to using a panel function, another option is to
use the function lineplot.CI, available in the package sciplot at
http://mutualism.williams.edu/sciplot
E.g.
# Define and generate variables in "long" format
range <- vector()
voice <- vector()
string <- strsplit(as
I couldn't help but respond to this one, it's not often I see my own name.
Using data from the survival library:
library(survival)
lung[1:10, c('time', 'status')]
Surv(lung$time, lung$status)[1:10]
--Matt
Matt Austin
Statistician
Amgen, Inc
800 9AMGEN9 x77431
805-447-7431
-Original Mess
try something like:
surv.data <- data.frame(times = rexp(100, 1/10), events = rbinom(100,
1, 0.7))
surv.data[sample(nrow(surv.data), replace = TRUE), ]
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of L
I have generated some some survival times and censoring indicators.
Thus I have an ordered pair for each observation. How do I sample these
ordered paris? I only know how to sample from a vector? I would
appreciate any help I could get.
Thanks
Matt
_
If you want to change the data,_fix_ will be a good option.But if you
just want to browse the data ,then invisible(edit(data)) is better
one.
I remember this question showed up some time ago.
2006/6/7, Ulrich Keller <[EMAIL PROTECTED]>:
> fix(data)
>
> will invoke edit(data) and store changes yo
"Antonio, Fabio Di Narzo" <[EMAIL PROTECTED]> writes:
> 1. If you have time to change internal functions naming, you can rename
> internal functions by putting a leading '.'.
> Even without namespace, I have noticed there is no check for corresponding
> docs for such functions.
>
> 2. If you don't
fix(data)
will invoke edit(data) and store changes you make in data without
displaying anything.
stat stat schrieb:
> Dear all R users,
>
> I have a query on "Edit" function. Suppose I have a data frame named
> "data". I can use EDIT function to see the materials contained in data, by
>
I am trying to learn how to use the graphics from the lattice package (
and am very new to R).
I am trying to replicate the example plot referenced below, by using the
lattice xyplot & lpolygon to create panels. I get what appears to be the
correct shape of the filled region, but cannot get t
On Wed, 7 Jun 2006, Nicolas Degallier wrote:
> I am trying to install in my R environment the rhdf5 package and
> library but it seems to have vanished from either the CRAN or
> BioConductors sites.
>
> Can you tell me where it would be possible to find it or any R
> library (or function) able to
I have three completely different(!) but workable solutions. Thanks much
for your help.
=Randy=
R. Zelick email: [EMAIL PROTECTED]
Department of Biology voice: 503-725-3086
Portland State University fax: 503-725-3888
mailing:
P.O.
Hi!
I am trying to install in my R environment the rhdf5 package and
library but it seems to have vanished from either the CRAN or
BioConductors sites.
Can you tell me where it would be possible to find it or any R
library (or function) able to read hdf files?
Sincerely,
Nicolas Degallier
I am trying to learn how to use the graphics from the lattice package (
and am very new to R).
I am trying to replicate the example plot referenced below, by using the
lattice xyplot & lpolygon to create panels. I get what appears to be the
correct shape of the filled region, but cannot get t
On Wed, 7 Jun 2006, stat stat wrote:
> Dear all R users,
>
> I have a query on "Edit" function. Suppose I have a data frame named
> "data". I can use EDIT function to see the materials contained in data,
> by using the command:
>
> > edit(data)
>
> But when I close the window then again the m
See: ?itsSubset
On 6/7/06, Ivan Kalafatic <[EMAIL PROTECTED]> wrote:
> If I understood correctly in irregular time series (its) objects, values are
> indexed by time stamps in POSIX format.
> But if I try to select the value of my time series corresponding to specific
> time stamp in the following
[Presuming you mean monotone; many matrices are already monotonous.]
General solutions to this problem are discussed in
@ARTICLE{Friendly:02:corrgram,
author = {M. Friendly},
title = {Corrgrams: Exploratory displays for correlation matrices},
journal = {The American Statistician},
yea
Dear all R users,
I have a query on "Edit" function. Suppose I have a data frame named "data".
I can use EDIT function to see the materials contained in data, by using the
command:
> edit(data)
But when I close the window then again the materials contained in data is
displayed
If I understood correctly in irregular time series (its) objects, values are
indexed by time stamps in POSIX format.
But if I try to select the value of my time series corresponding to specific
time stamp in the following way:
x - its object
i <- as.POSIXct("2006-05-19 15:30:00")
x[i,] or x[i] or
Federico Calboli wrote:
> Hi All,
>
> I am using plot(x, type = 'l') for some plotting, but I would like rounded
> edges
> rather than jagged edges in the plot (purely for aestetic reasons).
>
> How could I achieve that?
>
Perhaps you want something like:
x<-rnorm(50)
plot(spline(1:50,x),type=
Prof Brian Ripley wrote:
> It I understand you aright, that is done by par(lend) but the default is
> "round". So maybe your graphics device (unstated) on your OS (unstated)
> does not support this.
graphics device = X11 (xserver-xorg)
OS = Debian GNU/Linux, Kernel 2.4.27-2-686-smp
>
> We ne
On 6/7/06, Romain Francois <[EMAIL PROTECTED]> wrote:
> Le 07.06.2006 11:20, [EMAIL PROTECTED] a écrit :
> > Spencer Graves a écrit :
> >
> >
> >> I agree it would be great to sort the variables in a correlation
> >> matrix to make it easier to read and see patterns. I don't know any
> >> fu
On Wed, 7 Jun 2006, Federico Calboli wrote:
> Hi All,
>
> I am using plot(x, type = 'l') for some plotting, but I would like
> rounded edges rather than jagged edges in the plot (purely for aestetic
> reasons).
>
> How could I achieve that?
It I understand you aright, that is done by par(lend)
Le 07.06.2006 11:20, [EMAIL PROTECTED] a écrit :
> Spencer Graves a écrit :
>
>
>> I agree it would be great to sort the variables in a correlation
>> matrix to make it easier to read and see patterns. I don't know any
>> functions for doing that. If it were my problem, I might "order"
Dimitris Rizopoulos wrote:
> probably you want to use the `lend' argument of ?par(); I hope it helps.
Does not seem to work in my case.
F
--
Federico C. F. Calboli
Department of Epidemiology and Public Health
Imperial College, St Mary's Campus
Norfolk Place, London W2 1PG
Tel +44 (0)20 7594 1
probably you want to use the `lend' argument of ?par(); I hope it
helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuven, Belgium
Tel: +32/(0)16/336899
Fax: +32/(0)16/337015
Web:
Hi All,
I am using plot(x, type = 'l') for some plotting, but I would like rounded
edges
rather than jagged edges in the plot (purely for aestetic reasons).
How could I achieve that?
Cheers,
Federico
--
Federico C. F. Calboli
Department of Epidemiology and Public Health
Imperial College, St
Spencer Graves a écrit :
> I agree it would be great to sort the variables in a correlation
> matrix to make it easier to read and see patterns. I don't know any
> functions for doing that. If it were my problem, I might "order" the
> variables by their first principal component. There
Thanks everybody - it's working
Rainer
Dimitris Rizopoulos wrote:
> try something like:
>
> lapply(f, summary)
> sapply(f, function(x) AIC(logLik(x)))
>
>
> Best,
> Dimitris
>
>
> Dimitris Rizopoulos
> Ph.D. Student
> Biostatistical Centre
> School of Public Health
> Catholic University
1. If you have time to change internal functions naming, you can rename
internal functions by putting a leading '.'.
Even without namespace, I have noticed there is no check for corresponding
docs for such functions.
2. If you don't want to rename all internal functions, the best way is
writing an
Dan Rabosky wrote:
> Hello.
>
> I am creating an R package that I'd like to submit to CRAN (OS Windows
> XP). How do I distinguish among 'public' functions, e.g., those that are
> intended to be called by users of the package and for which I am providing
> documentation & examples, and 'privat
It's confusing to have time values that are different but `match', so
first you want to reduce your times to hh:mm (e.g. with substr).
Then you can cbind() the data frames, and use tapply() to do the
averaging. E.g.
d1 <- read.table("foo1")
d2 <- read.table("foo2")
d <- cbind(d1, d2)
d$time <-
69 matches
Mail list logo