Displaying 1 result from an estimated 1 matches for "psdplot".
Did you mean:
cdplot
2012 Aug 05
1
trouble with looping for effect of sampling interval increase
...lUpperQ <- NULL
IntervalMean <- NULL
IntervalMedian <- NULL
#Here is the loop for betas from sampling interval increase
c <- WHOLESIZE[2]-1
for (i in 1:c)
{
x <- length(unique(index[,i]))
for (j in 1:x)
{
data <- WHOLE [WHOLE[,x]==j,1]
#get power spectral density
PSDPLOT <- spectrum (data, detrend = TRUE, plot = FALSE)
frequency <- PSDPLOT$freq
PSD <- PSDPLOT$spec
#log transform the power spectral density
Logfrequency <- log(frequency)
LogPSD<- log(PSD)
#fit my line to the data
Line <- lm (LogPSD ~ Logfrequency)
#store the slope of the l...