density plots using density.lf, data.frame and sort.int errors

An embedded and charset-unspecified text was scrubbed...
Name: not available
URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20120531/1337f8da/attachment.pl>
My guess (unconfirmed) is that read.table() gives you a data frame but
density.lf expects an atomic (= not a list = not a data frame) vector.
Perhaps try density.lf(x[,1]) to just send the column -- the "drop"
behavior should make sure this is an atomic vector.

If that doesn't help, please do provide us with the output of

dput(head(x, 30))

and the package from which the density.lf function comes from.

Best,
Michael

On Thu, May 31, 2012 at 11:35 AM, Andrea S Sequeira
Dear R help group: I am attempting to produce a density plot from a list of
20000 values using the density.lf function and would appreciate any help, I
hope I have done my homework reading the documentation but I still seem to
be missing something basic.

I have read the data as a table using read.table, with header=TRUE (I
excluded 2000 values), when calling the objects it appears to be there and
I can see the values

this is what I get when doing

density.lf (x=logN0, n=50, window="gaussian")

Error in `[.data.frame`(x, order(x, na.last = na.last, decreasing =
decreasing)) :

?undefined columns selected

then assigned a column name (using colname so it is called first), then
assigned it as a vector using

assign ("x", c (logN0))

the error I get is

density.lf (x, n=50, window="gaussian")

Error in sort.int(x, na.last = na.last, decreasing = decreasing, ...) :

?'x' must be atomic

the traceback produces:

traceback ()
5: stop("'x' must be atomic")

4: sort.int(x, na.last = na.last, decreasing = decreasing, ...)

3: sort.default(x)

2: sort(x)

1: density.lf(x, n = 50, window = "gaussian")

Thanks in advance, Andrea

--
Andrea Sequeira
Associate Professor
Department of Biological Sciences
Wellesley College, Wellesley MA 02481

? ? ? ?[[alternative HTML version deleted]]

______________________________________________
R-help at 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.