Displaying 1 result from an estimated 1 matches for "logpsd".
Did you mean:
logpd
2012 Aug 05
1
trouble with looping for effect of sampling interval increase
...e(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 line
Betas <- rbind (Betas, -coef(Line)[2])
#Get values on the curve shape
BSkew <- skew (Betas)
BMean <- mean (Betas)
BMedian <- median (Betas)
Q <- quantile (Betas)
#store...