Dear R experts,
I was wondering if there are any R functions that give the tail area
of a sum of chisquare distributions of the type:
a_1 X_1 + a_2 X_2
where a_1 and a_2 are constants and X_1 and X_2 are independent chi-square
variables with different degrees of freedom.
Thanks,
Klaus
Hi All
I have read this document - An Introduction to the .C Interface to R,
which primarily tells how to interface C language with R.
Is there and more elaborative and online documentation regarding this
interface.
Any pointers appreciated
-thanks
-gaurav
=
Dear All,
I am Ricky Jacob, a project Student from India who is working on
Forest Inventories.
Input data:
Plot(area = .1 ha) data having the following information:
1) Basal Area
2)Tree Density
3)Volume
So I am applying this information to the corresponding pixels in the
satellite imagery of the st
John Fox wrote:
> Dear Horace,
>
> The Bonferonni p-value is obtained from the "unadjusted" p-value by
> multiplying the latter by the number of observations, and provides a
> conservative (although usually quite accurate) outlier test. When the
> adjusted p-value exceeds 1 you can take that as an
Thank you Gabor, your solution works great, even for veeery long expressions.
//Daniel
On 3/28/07, Gabor Grothendieck <[EMAIL PROTECTED]> wrote:
> The solution in my post has the advantage of not using
> eval or character conversions (except in setting up
> L where you want such character convers
hi,
I want to use EM for MLE estimation.WHAT would be the relevant package or
code??
Regards,
Subhajit.
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
Hi,
On Wednesday 28 March 2007 07:32:08 pm Paul Murrell wrote:
> You could call the grid function getNames() and if that is zero-length,
> then the device is either blank or only contains traditional graphics.
Thanks Dr. Murrell. I was looking for something simple like this.
>
> > A related ques
Dear Horace,
The Bonferonni p-value is obtained from the "unadjusted" p-value by
multiplying the latter by the number of observations, and provides a
conservative (although usually quite accurate) outlier test. When the
adjusted p-value exceeds 1 you can take that as an indication that there are
n
On Wed, 28 Mar 2007, Frank E Harrell Jr wrote:
> Eric Elguero wrote:
>> Hi everybody,
>>
>> recently I had to teach a course on Cox model, of which I am
>> not a specialist, to an audience of medical epidemiologists.
>> Not a good idea you might say.. anyway, someone in the
>> audience was very ho
I did read your second message about the problem symptoms disappearing
but I thought that I might make a couple of suggestions about your
code anyway.
There are a number of helper functions declared in Rinternals.h such
as ScalarReal, which is equivalent to your mkans. (Also
ScalarInteger, Scalar
Thanks Christos for this tip on the book. I needed something like that
get started on the topic.
AA.
- Original Message -
From: "Christos Hatzis" <[EMAIL PROTECTED]>
To: "'John Kane'" <[EMAIL PROTECTED]>; "'Chuck Cleland'"
<[EMAIL PROTECTED]>
Cc: "'R R-help'"
Sent: Wednesday, March 28, 2
Andy Jacobson said the following on 3/28/2007 4:01 PM:
> Hi,
>
> I'd like to try using "options(error=recover) in my ~/.Rprofile, but
> it appears that the function "recover" is not defined during R
> startup when the .Rprofile is processed. "recover" is defined after
> I get an R prompt,
Hi
Prasenjit Kapat wrote:
> OK, here is the issue:
> I have a graphcis device (say X11) which I can set as active by
> dev.set(n). Is there ANY way to identify/infer whether this X11 device
> was obtained from a lattice function (like xyplot, levelplot...) or a
> standard function (like plot, sca
On Wed, 2007-03-28 at 08:38 -0700, Brian Dolan wrote:
> Hello,
>
> I would like to take the string
>
> "mylist = list(a = 5, b = 7)"
>
> and evaluate it as a list. I have attempted to use parse and several
> other functions with no success.
>
> Thanks for your time.
Not sure what combination
Hi,
I'd like to try using "options(error=recover) in my ~/.Rprofile, but
it appears that the function "recover" is not defined during R
startup when the .Rprofile is processed. "recover" is defined after
I get an R prompt, however. Can anyone shed light on why this is,
and whether a work
cmd = "mylist = list(a = 5, b = 7)"
(eval(parse(text=cmd)))
b
On Mar 28, 2007, at 11:38 AM, Brian Dolan wrote:
> Hello,
>
> I would like to take the string
>
> "mylist = list(a = 5, b = 7)"
>
> and evaluate it as a list. I have attempted to use parse and
> several other functions with no suc
Hi folks,
I use the outlier.test in package car to test a lm model and the bonferroni p
value returned is shown as NA. When the object is typed it indicates the p
value is greater than 1. I'm not sure how to interpret it.
Thanks in advance.
Horace W. Tso
> outlier.test(mod)$test
max|rstuden
Hello,
I would like to take the string
"mylist = list(a = 5, b = 7)"
and evaluate it as a list. I have attempted to use parse and several other
functions with no success.
Thanks for your time.
-brian dolan
~~~
may all your sequences converge
[[alternative HTML version deleted]]
Marc Schwartz <[EMAIL PROTECTED]> writes:
> On Wed, 2007-03-28 at 19:55 -0200, Alberto Monteiro wrote:
>> Prof Brian Ripley wrote:
>> >
>> > We have already seen three solutions.
>> >
>> > I don't like to see the use of c() for its side effects. In this
>> > case Marc's as.vector seems to me to
Prof Brian Ripley wrote:
>
> We have already seen three solutions.
>
> I don't like to see the use of c() for its side effects. In this
> case Marc's as.vector seems to me to be self-explanatory, and that
> is a virtue in programming that is too often undervalued.
>
I agree; but for our enligh
Hi.
The dprep library has at least three different methods for outlier
detection: baysout, mahaout, robout.
I wanted to test them on a very simple data set:
vrmat<-cbind((1:22),c(8,14,14,17,21,20,27,23,25,33,31,32,30,36,37,40,42,44,52,61,81,265))
As you can see by eyeballing this, the last poi
On Wed, 2007-03-28 at 19:55 -0200, Alberto Monteiro wrote:
> Prof Brian Ripley wrote:
> >
> > We have already seen three solutions.
> >
> > I don't like to see the use of c() for its side effects. In this
> > case Marc's as.vector seems to me to be self-explanatory, and that
> > is a virtue in
Hi again!
Please try something like this:
> zone1.df
zone sex len
11 male 15
21 fema 20
31 fema 17
42 fema 19
52 male 18
> str(zone1.df)
'data.frame': 5 obs. of 3 variables:
$ zone: Factor w/ 2 levels "1","2": 1 1 1 2 2
$ sex : Factor w/ 2 levels "fema","male": 2
Thankyou Achim.
Along these lines, how might I extract subsets from that set by day of week?
eg: create a new series that includes Mondays only.
Thanks,
Alf
- Original Message -
From: "Achim Zeileis" <[EMAIL PROTECTED]>
To: "Alfonso Sammassimo" <[EMAIL PROTECTED]>
Cc:
Sent: Thursday, Ma
Hello:
Is EM clustering the same as latent class analysis?
Thank you.
Zhongmiao Wang
Senior Analyst
RMG Connect
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE d
Thanks to everyone for the information and to Chuck
for the links. I must have misspelt something because
I missed most of them in the RsiteSearch.
Anyway back to selling R :)
--- Chuck Cleland <[EMAIL PROTECTED]> wrote:
> John Kane wrote:
> > A colleague was asking me if R does multi-level
> >
[EMAIL PROTECTED] wrote:
> Hi!
> I’m trying to calculate de C index (concordance probability) through the
> somers2 function (library Hmisc). I’m interesting on including the
> sampling effort as a weight factor for the evaluation of model predictions
> with real data. I’ve some questions about tha
Is this something like what you want?
ab <- rep(1,4)
bb <- rep(2,4)
cc <- rep(3,4)
mydata <- data.frame(ab,bb,cc)
unlist(mydata)
unlist(data.frame(t(mydata)))
--- A Ezhil <[EMAIL PROTECTED]> wrote:
> Hi,
>
> I have a matrix HR(9x27). I would like to make a
> single vector with elements: t(HR
have you tried replacing
.duplicate()
in your my.causality() by
vars:::.duplicate()
?
b
On Mar 28, 2007, at 3:11 PM, Leeds, Mark ((IED)) wrote:
> I am using the vars package and it calls a function causality() which
> then calls something
> called .duplicate. I had to modify the causality fu
Hi All,
I'm trying to use the geoRglm package to run a poisson spatial glm
on a dataset with several covariates. When I run without covariates I
have no problems.
control1.data.geo <- mcmc.control(S.scale=0.2, thin = 1)
model1.data.geo <- list(cov.pars = c(1,1), beta=c(1), family="poisson"
OK, here is the issue:
I have a graphcis device (say X11) which I can set as active by
dev.set(n). Is there ANY way to identify/infer whether this X11 device
was obtained from a lattice function (like xyplot, levelplot...) or a
standard function (like plot, scatterplot,...).
A related question is
I am using the vars package and it calls a function causality() which
then calls something
called .duplicate. I had to modify the causality function slightly for
my purposes and
I called it my.causality() but now the .duplicate function is no longer
known to the my.causality function.
I'm fairly c
Thanks Christos
Some low life managed to get it out of the library!
It looks worth reading, I must say.
--- Christos Hatzis <[EMAIL PROTECTED]> wrote:
> Also, take a look at Andrew Gelman's recent book
> with Jennifer Hill on
> multilevel modeling.
> It is written as a practical how-to text book
Also, take a look at Andrew Gelman's recent book with Jennifer Hill on
multilevel modeling.
It is written as a practical how-to text book using R as the primary
platform. I have found it very easy to read, extremely useful and very good
value for its price (I am not related to the authors or the p
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
No one can help you.
Ranjan
On Wed, 28 Mar 2007 11:47:03 -0600 (MDT) [EMAIL PROTECTED] wrote:
> Hi all
>
> I was wondering if someone could help me.
We have already seen three solutions.
I don't like to see the use of c() for its side effects. In this case
Marc's as.vector seems to me to be self-explanatory, and that is a virtue
in programming that is too often undervalued.
On Wed, 28 Mar 2007, Marc Schwartz wrote:
> On Wed, 2007-03-28 at
Hi all
I was wondering if someone could help me.
I have to estimate some parameters, so I am using the function nlm. Inside
this function I have to integrate, hence
I am using the function adapt.
I don't understand why it is giving the following warnings:
At the beginning:
Warning: a final empt
That is not the only way to apply logistic regression to this problem
(although it is a common error in the analysis of cancer studies).
One can discretize time and apply logistic regression to survival over
each short time period (jointly): doing so comes pretty close to what the
Cox proporti
A Ezhil wrote:
>
> I have a matrix HR(9x27). I would like to make a
> single vector with elements: t(HR[,1]) followed by
> t(HR[,2]) and then t(HR[,3] ... etc. Is there any neat
> way of converting this matrix into a vector rather
> doing something like c(t(HR[,1]), t(HR[,2]), t(HR[,3])
> ..)?
>
On Wed, 2007-03-28 at 10:27 -0700, A Ezhil wrote:
> Hi,
>
> I have a matrix HR(9x27). I would like to make a
> single vector with elements: t(HR[,1]) followed by
> t(HR[,2]) and then t(HR[,3] ... etc. Is there any neat
> way of converting this matrix into a vector rather
> doing something like c(t
A matrix is just a vector with a dim attribute, so if you remove
the dim attribute it'll be a vector. The elements of a matrix
are stored columnwise so you'll get just what you want, if i get
your question right.
> g <- matrix(1:6, nc=2, nr=3)
> g
[,1] [,2]
[1,]14
[2,]25
[3,
c(HR) will do it.
Ravi.
---
Ravi Varadhan, Ph.D.
Assistant Professor, The Center on Aging and Health
Division of Geriatric Medicine and Gerontology
Johns Hopkins University
Ph: (410) 502-2619
Fax: (410) 614-96
Hi,
I have a matrix HR(9x27). I would like to make a
single vector with elements: t(HR[,1]) followed by
t(HR[,2]) and then t(HR[,3] ... etc. Is there any neat
way of converting this matrix into a vector rather
doing something like c(t(HR[,1]), t(HR[,2]), t(HR[,3])
..)?
Thanks in Advance.
Kind reg
How I test for differences among time points without assuming
parallelism, With a profile analysis model, Using manova R function.
For example, with Potthoff and Roy data (Potthoff and Roy 1964) how
I can test
H: ABC=0,
when
A<-diag(2)
and
C<-rbind(diag(3),rep(-1,3))
(Davis, 2005 p.p 81 ).
John Kane wrote:
> A colleague was asking me if R does multi-level
> modelling as opposed to multiple regression. Since I
> have no knowledge of multi-level modelling (except 5
> minutes googling ) I thought that I would as here.
>
> Does are offer any multi-level modeling packages? It
> looked
On the same point, transforming time-to-event data to binary outcomes so
that contingency-table analysis (odds ratios etc) or logistic regression can
be applied will result in loss of information that could lead to misleading
conclusions.
For example, assuming that there is a good-prognosis group
Yes, typically one uses lme4 or nlme to do multilevel modeling in R.
--
Andrew J Perrin - andrew_perrin (at) unc.edu - http://perrin.socsci.unc.edu
Assistant Professor of Sociology; Book Review Editor, _Social Forces_
University o
A colleague was asking me if R does multi-level
modelling as opposed to multiple regression. Since I
have no knowledge of multi-level modelling (except 5
minutes googling ) I thought that I would as here.
Does are offer any multi-level modeling packages? It
looked like arm might be one but I wa
Hi!
Im trying to calculate de C index (concordance probability) through the
somers2 function (library Hmisc). Im interesting on including the
sampling effort as a weight factor for the evaluation of model predictions
with real data. Ive some questions about that: first of all Im not
really sure
You can (and I have) fit survival data with logistic regression. Agresti (1990,
pp 189--196) has an introductory discussion.
The issue is whether the occurrence of the event is of interest or whether the
time-to-event is of interest. If the study lasts 180 days (as in my case)
logistic regress
The answer is correct, and is the way it should be. Cluster indicators are only
nominal: there is no ordering in the ids, and hence the means are in different
order.
Either way, the reason for this (and you could get totally different answers
also) is because of how R initializes kmeans (random
On Wed, 2007-03-28 at 10:49 -0500, Douglas Bates wrote:
> On 3/28/07, Marc Schwartz <[EMAIL PROTECTED]> wrote:
> > On Wed, 2007-03-28 at 02:42 -0700, francogrex wrote:
> > > Does anyone know of an R package that is equivalent of S+SeqTrial for
> > > analysis of clinical trials using group sequentia
[EMAIL PROTECTED] wrote:
> Mich beschäftig folgende Fragestellung. Ich kenne die Verteilung
> (lognormal) zusätzlich weiss ich das 99%, das 90% und das 1% Quantil.
> Gibt es in R eine Möglichkeit die Lognormalverteilung zu finden, das
> heisst den korrespondierenden logmean und logsd?
>
(Nich
On 3/28/07, Marc Schwartz <[EMAIL PROTECTED]> wrote:
> On Wed, 2007-03-28 at 02:42 -0700, francogrex wrote:
> > Does anyone know of an R package that is equivalent of S+SeqTrial for
> > analysis of clinical trials using group sequential methods? Thanks.
>
> I don't know that there are fully functio
Dear R-Helpers,
I performed kmeans clustering on the following data set(y):
YEAR PRODUCTS
1 10
2 42
3 25
4 42
5 40
6 45
7 44
8 47
9 42
with this code:
cluster<-kmeans(y[,c("YE
Eric Elguero wrote:
> Hi everybody,
>
> recently I had to teach a course on Cox model, of which I am
> not a specialist, to an audience of medical epidemiologists.
> Not a good idea you might say.. anyway, someone in the
> audience was very hostile. At some point, he sayed that
> Cox model was use
On Wed, 28 Mar 2007, Alfonso Sammassimo wrote:
> Is there a way of aggregating 'zoo' daily data according to day of week? eg
> all Thursdays
Sure, the easiest way will probably differ depending on the time stamp
class. One example might be this:
## small example with Date index
z <- read.z
Hi everybody,
recently I had to teach a course on Cox model, of which I am
not a specialist, to an audience of medical epidemiologists.
Not a good idea you might say.. anyway, someone in the
audience was very hostile. At some point, he sayed that
Cox model was useless, since all you have to do is
This mailing list is in english
If you know that for the quantiles, say:
1%: 3
90%: 5.2
99%: 6.1
then you could quickly look for the minimum of the following function
foo <- function(x){
(qlnorm(0.01, meanlog = x[1], sdlog = x[2]) - 3)^2 +
(qlnorm(0.9, meanlog = x[1], sdlog = x
This could work, but not with big matrix!
year <- c(1990,1995,2000,2005)
Prod <- c(2478,3192,NA,1594)
matrix <- data.frame(cbind(year,Prod))
for (i in 1:dim(matrix)[1])
{
if (is.na(matrix[i,2])) {matrix[i,2] <- 0}
}
__
R-help@st
On Wed, 28 Mar 2007, Jeff Lusk wrote:
> After reading the help files (?Startup) and using
> RSiteSearch("defaultPackages"), I have been trying to add several packages
> to my default startup list using the following code:
>
> local({
> old<-getOption("defaultPackages")
> options(defaultPackages=c(
On Wed, 28 Mar 2007, [EMAIL PROTECTED] wrote:
>
> Dear all,
>
> I have some C code using the zlib uncompress. It compiles fine under Win XP
> (with SP2) with R up to 2.3.1 and generates four files (Makedeps, *.d,
> *.dll, and *.o):
No one every said that zlib interface would be exported from R, a
Since R-2.4.0 less symbols are exported from R.dll.
Uwe Ligges
[EMAIL PROTECTED] wrote:
> Dear all,
>
> I have some C code using the zlib uncompress. It compiles fine under Win XP
> (with SP2) with R up to 2.3.1 and generates four files (Makedeps, *.d,
> *.dll, and *.o):
>
> C:\Data\R-packages
See ?is.na and use its result for indexing.
Uwe Ligges
Sergio Della Franca wrote:
> Dearl R-Helpers,
>
>
> I have the following data set:
>
> YEAR PRODUCTS
> 1990 2478
> 1995 3192
> 2000 NA
> 2005 1594
>
>
> I wanto to replace NA values, in the PRODUCTS column, with 0.
>
>
On Wed, 2007-03-28 at 16:21 +0200, Sergio Della Franca wrote:
> Dearl R-Helpers,
>
>
> I have the following data set:
>
> YEAR PRODUCTS
> 1990 2478
> 1995 3192
> 2000 NA
> 2005 1594
>
>
> I wanto to replace NA values, in the PRODUCTS column, with 0.
>
>
> How can i obtain t
After reading the help files (?Startup) and using
RSiteSearch("defaultPackages"), I have been trying to add several packages
to my default startup list using the following code:
local({
old<-getOption("defaultPackages")
options(defaultPackages=c(old,"lattice","RODBC")
})
Subsequently, when I qu
Mich beschäftig folgende Fragestellung. Ich kenne die Verteilung
(lognormal) zusätzlich weiss ich das 99%, das 90% und das 1% Quantil.
Gibt es in R eine Möglichkeit die Lognormalverteilung zu finden, das
heisst den korrespondierenden logmean und logsd?
Vielen Dank für ihre Hilfe
Gruss
Yvonne
_
Dearl R-Helpers,
I have the following data set:
YEAR PRODUCTS
1990 2478
1995 3192
2000 NA
2005 1594
I wanto to replace NA values, in the PRODUCTS column, with 0.
How can i obtain this?
Thak you in advance.
Sergio Della Franca
[[alternative HTML version deleted]
Dear all,
I have some C code using the zlib uncompress. It compiles fine under Win XP
(with SP2) with R up to 2.3.1 and generates four files (Makedeps, *.d,
*.dll, and *.o):
C:\Data\R-packages\cmg>rcmd SHLIB cmg.c -o cmg.dll
latex: not found
making cmg.d from cmg.c
gcc -IC:/PROGRA~1/R/R-2.3.1/
Alberto Monteiro wrote:
> Pete Cap wrote:
>
>> I'm writing a tcl/tk gui tool to wrap around RMySQL for some co-workers.
>>
>>
> Good luck; I find the documentation on the library(tcltk) very
> poor, lacking examples for most of the functions.
>
>
Well, they are interface functions, and w
On Wed, 2007-03-28 at 15:46 +0300, [EMAIL PROTECTED] wrote:
> Hallo,
>
> I'm trying to sample a matrix with simple random sampling without
> replacement but seem to have a problem with the matrix length. Both
> sample() and srswor() use length() which returns the number of columns in
> matrix. Thi
Dear All:
If y~Ga(a,b)
What's about the (n-1)y~Ga(a, (n-1)b), is this right?
Many Thanks
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the po
Thanks Brian.
Working for one of these large corporations, I receive a considerable amount of
mails requesting me to verify whether software I use (and which was made
avalailable via company wide installers on my request) is Vista compatible.
It's not that these big companies will roll out Vis
Gabor Grothendieck wrote:
>
> Try this:
>
> library(gsubfn)
> pat <- "([[:upper:]][[:lower:]]*) "
> s <- "Aaaa 3 x 0 Bbbb"
>
> strapply(s, pat, backref =-1)[[1]]
>
Ok, I think I got it.
For example, if I want to retrieve Aaaa and Bbbb (which
could be any strings) or the numbers 3 and 0, I would
Thank you Gabor. This works. And you are right, bxp is actually
on the boxplot help page. I am sorry I missed it. Thanks again.
A.
- Original Message -
From: "Gabor Grothendieck" <[EMAIL PROTECTED]>
To: "AA" <[EMAIL PROTECTED]>
Cc:
Sent: Tuesday, March 27, 2007 7:49 PM
Subject: Re: [R] l
On Wed, 28 Mar 2007, [EMAIL PROTECTED] wrote:
> The Vista issue is not innocent as it threatens the life of R within
> large corporations. So any posts on how R runs under Vista and what has
> to be done to make it work and what cannot be done etc will be very
> useful.
See the rw-FAQ in the R
Please read the last line on every post to r-help.
On 3/28/07, Alfonso Sammassimo <[EMAIL PROTECTED]> wrote:
> Is there a way of aggregating 'zoo' daily data according to day of week? eg
> all Thursdays
>
> I came across the 'nextfri' function in the documentation but am unsure how
> to change thi
On Wed, 2007-03-28 at 02:42 -0700, francogrex wrote:
> Does anyone know of an R package that is equivalent of S+SeqTrial for
> analysis of clinical trials using group sequential methods? Thanks.
I don't know that there are fully functional equivalents, but you might
want to look at the following:
On Wed, 2007-03-28 at 13:22 +0200, Schmitt, Corinna wrote:
> Hallo,
>
> Sorry for more details before but I was in a hurry.
>
> One workspace comes from Matlab and one should come from Bioconductor.
> At the moment I can just work with the Matlab-workspace. I just know
> that the data is stored
Try this:
library(gsubfn)
pat <- "([[:upper:]][[:lower:]]*) "
s <- "Aaaa 3 x 0 Bbbb"
strapply(s, pat, backref =-1)[[1]]
or (not quite as general but works in this case):
pat <- "([[:upper:]][[:lower:]]*) "
s <- "Aaaa 3 x 0 Bbbb"
s.idx <- gregexpr(pat, s)[[1]]
substring(s, s.idx, s.idx + attr(s
Hallo,
I'm trying to sample a matrix with simple random sampling without
replacement but seem to have a problem with the matrix length. Both
sample() and srswor() use length() which returns the number of columns in
matrix. This means that to function it seems that the sample size exceeds
the matri
Dear list,
I use a lot of data files (at the moment 35 files by 34 MB each) for analysis.
My own processing steps will (hopefully) improve with increasing number of
files I have.
So I am looking for an suitable workstation configuration, e.g. 64bit
architecture, 16 GB RAM etc.
Does anyone
Hi folks,
The Vista issue is not innocent as it threatens the life of R within large
corporations. So any posts on how R runs under Vista and what has to be done to
make it work and what cannot be done etc will be very useful.
Chris
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[
Is there a way of aggregating 'zoo' daily data according to day of week? eg
all Thursdays
I came across the 'nextfri' function in the documentation but am unsure how
to change this so any day of week can be aggregated.
I have used POSIX to arrange the data (not as 'zoo' series) according to day
This works:
grep("([A-Za-z]*) ", "Aaaa 3 x 0 Bbbb") # 1
This also works:
grep("([A-Za-z]*) ", "Aaaa 3 x 0 Bbbb", value=T) # Aaaa 3 x 0 Bbbb
However, I want a grep that returns the _matched_ pattern, which, in this
case, would be Aaaa. How can I do it?
Alberto Monteiro
___
I am not 100% sure of what you want but maybe this will help. I
have modified your example data to omit one point so that there
is a holiday among the rows:
library(zoo)
# read data
Lines <- "DateOpenHigh Low CloseNifty
2004-01-01 1880.35 1917.05 1880.351912.25
2004-01-02 1912.25
The solution in my post has the advantage of not using
eval or character conversions (except in setting up
L where you want such character conversions to
build up the names, as your more general sitution
shows). The following is the same as in that post except
the line setting L in that post is re
Hello All
I have time series data like given hereinbelow :-
The data is having dates. I tried using ts() but it requires frequency. As
i have holidays, saturdays sunday market closed, and also sometimes the
data may not be recorded. I want to create a time series object but i want
that where e
Hallo,
Sorry for more details before but I was in a hurry.
One workspace comes from Matlab and one should come from Bioconductor.
At the moment I can just work with the Matlab-workspace. I just know
that the data is stored in an Excel-File --> I have a table.
I know that there exists an RMatlab
> On 3/25/07, John Bullock
> <[EMAIL PROTECTED]> wrote:
>>
>> I am trying to eliminate panel borders from my
>> lattice plots. By default, they always print.
>> For
>> example:
>>
>>library(lattice)
>>x <- seq(-3,3,length=1000)
>>y1 <- dnorm(x)
>>y2 <- dnorm(x, sd=.5)
>>data <-
[EMAIL PROTECTED] wrote:
> Hi All
>
> Does R supports Vista Business and Vista Home Premium. I am planning to
> upgrade my system. Before that I just wanted to confirm whether R will
> work on Vista or not.
Although R Core does not officially support Vista, AFAIK, we have heard
of working in
Hi All
Does R supports Vista Business and Vista Home Premium. I am planning to
upgrade my system. Before that I just wanted to confirm whether R will
work on Vista or not.
Thanks for this help in advance
-gaurav
Sayonara With Smile & With Warm Regards :-)
G a u r a v Y a d a v
Dear R-Helpers,
I want to perform a kmeans clustering with ward method.
How can i obtain this?
Thank yuo in advance.
Serigo Della Franca
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mai
Does anyone know of an R package that is equivalent of S+SeqTrial for
analysis of clinical trials using group sequential methods? Thanks.
--
View this message in context:
http://www.nabble.com/R-equivalent-of-S%2BSeqTrial--tf3478858.html#a9708877
Sent from the R help mailing list archive at Nabb
you can still use scale() (as you have been told), look at the help
page for more info, especially at the Arguments section, e.g.,
mat <- matrix(rnorm(100*10), 100, 10)
rng <- apply(mat, 2, range)
scale(mat, scale = rng[2, ] - rng[1, ])
or you could even use apply() directly, e.g.,
apply(mat,
Now it works, I made no changes to the functions in question, must
have been some other mistake in my code and I obviously misinterpreted
the crash report.
defineVar(install(str),mkans(x),rho) works fine. The content of the
variable 'str' is assigned in the given environment in R as it should.
I
Hi sergio,
Sergio Della Franca wrote:
> Dear R-Helpers,
>
>
> I want to perform a standardization of a variable with range method.
>
> i.e.:
>
> Standardization (range) ==> (var-min(var))/(max(var)-min(var))
>
>
> Do you konw how can i develop this?
>
>
As you do ... but don't use var which is
Dear R-Helpers,
I want to perform a standardization of a variable with range method.
i.e.:
Standardization (range) ==> (var-min(var))/(max(var)-min(var))
Do you konw how can i develop this?
Thank you in advance.
Sergio Della Franca
[[alternative HTML version deleted]]
__
Hi everybody, I'm a newbye with lots of problems :). I'm trying to use
nlsystemfit, but I recieve two error messages whose origin that I don't
understand.
1) When I try to reproduce the example reported in the systemfit package
manual, that is
library( systemfit )
data( ppine )
hg.formula <-
Thank you very much for your good suggestions.
I have chosen to pursue the suggestion by Peter which worked like a dream :)
However, my problem is slightly more complicated still. I apologize
for not mentioning this in the initial question.
I have to do the evaluation inside a loop, not knowing ex
100 matches
Mail list logo