You will find here a reprodusible example, and you will also see it "exploding"
myData<-seq(1:100000)
windowList<-list(1,2,3,4)
totalLength<- length(myData)
windowSize<-30
for (i in 2:length(windowList)) {
myMean<-array()
indexes<-seq(1,totalLength,by=windowSize*windowList[[i]])
for (j in c(1,length(indexes))){
myMean[j]<-mean(myData[indexes[j]:indexes[j+1]]) # Select windowSize
number of elements
# Then I will call a cdf and cast lines()
}
}
as you can see on a given vector I want to calculating the mean value by taking
each time different number of elements. The reason is I want to see how much
increasing the number of samples affects the cdf distribution
I would like to thank you in advance for your help
B.R
Alex
________________________________
From: R. Michael Weylandt <michael.weyla...@gmail.com>
Cc: R help <R-help@r-project.org>
Sent: Monday, March 12, 2012 4:35 PM
Subject: Re: [R] Window on a vector
It's not always true that going out of bounds in subscripting gives an
error (as you seem to assert in your original post)
x <- 1:3
x[4] # No error
and note that mean() has a na.rm argument.
Perhaps you should construct a *reproducible* example of what you
think will go wrong.
Michael
> Yep I did some mistake while I was typing the example
> What I want is this one
> cbind(2^(0:(n-1))) and the problem is that at the last "window" this is
> going to explode are there will be less elements than what the window asks
> for.
>
> How should I deal with that?
> Regards
> Alex
>
>
>
> ________________________________
> From: David Winsemius <dwinsem...@comcast.net>
>
> Cc: R help <R-help@r-project.org>
> Sent: Saturday, March 10, 2012 2:49 PM
> Subject: Re: [R] Window on a vector
>
>
> On Mar 10, 2012, at 7:44 AM, Alaios wrote:
>
>> Dear all,
>> I have a large vector (lets call it myVector) and I want to plot its value
>> with the logic below
>>
>> yaxis<-myVector[1]
>> yaxis<-c(xaxis,mean(myvector[2:3])
>> yaxis<-c(xaxis,mean(myvector[4:8])
>> yaxis<c(xaxis,mean(myvector[9:16])
>> yaxis<c(xaxis,mean(myvector[17:32])
>>
>> this has to stop when the new ..... yaxis<c(xaxis,mean(myvector[1024:2048])
>> will not find the correspondent number of elements, either wise it will stop
>> with an error.
This doesn't seem true:
x <- 1:3
x[4] # Not an error
>>
>>
>> How I can do something like that in R?
>
> This will generate two series that are somewhat like your index
> specification. I say "somewhat" because you appear to have changed the
> indexing strategy in the middle. You start with 2^0. 2^1 and 2^2 as you
> "begin" but then switch to 2^3+1, and 2^4+1.
>
> n=20
> cbind(2^(0:(n-1)), 2^(1:n)-1)
>
> You can decide what to use for n with logic like:
>
> which.max(20 >= 2^(1:10) )
>
> Then you can use sapply or mapply.
>
>
>> Alex
>> [[alternative HTML version deleted]]
>
> Please learn to post in plain text.
>
> --
> David Winsemius, MD
> West Hartford, CT
> [[alternative HTML version deleted]]
>
>
> ______________________________________________
> 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.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>
[[alternative HTML version deleted]]
______________________________________________
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.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
