Displaying 1 result from an estimated 1 matches for "logpk".
Did you mean:
logk
2012 Apr 09
1
Panel.abline would not show beyond a certain slope value
...convert to match the names(freq), the character-based hash key of freq, which is the value
if (!(ichar %in% names(kfreq)))
kfreq[ichar] <- 0
}
sortedkeys <- as.character(k)
kfreq <- kfreq[sortedkeys]
pk <- kfreq / length(y)
logk <- log(k)
logpk <- log(pk)
logpk[logpk == -Inf] = "" # remove the -Inf or log(p(k)) = 0 values for lm function, NULL is 0-length so use "" instead that has length of one null character
logpk <- as.numeric(logpk) # "" converts all values to character, lm needs numeric...