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