Instead of a for loop, why not use the vectorization inherent in R?
sigmasqaured <- 1
i <- complex(real = 0, imaginary =1)
f <- seq(0,0.5,0.1)
spectrum
<-
(sigmasqaured)/(abs(1-2.7607*exp(2*pi*i*f)+3.8106*exp(4*pi*i*f)-2.6535*exp(6*pi*i*f)+0.9258*exp(8*pi*i*f))^2)
> spectrum
[1] 9.632720e+00 1.411130e+03 2.947753e+00 6.479994e-02 1.295175e-02
8.042731e-03
On Tue, Feb 7, 2012 at 1:08 PM, Jaymin Shah <jayminsh...@live.com> wrote:
> I have mad a for loop to try and output values which i have named
> spectrum. However, I cannot seem to get the answers to come out as a
> vector which is what i need. They come out as separate values which I am
> then unable to join together. Thank you
>
> for(f in seq(0,0.5,0.1)) {
> sigmasqaured <- 1
> i = complex(real = 0, imaginary = 1)
> spectrum <-
> (sigmasqaured)/(abs(1-2.7607*exp(2*pi*i*f)+3.8106*exp(4*pi*i*f)-2.6535*exp(6*pi*i*f)+0.9258*exp(8*pi*i*f))^2)
> print(spectrum)
> }
> [[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.
