Displaying 3 results from an estimated 3 matches for "miscl".
Did you mean:
misc
2010 Jul 27
2
lattice: option to sort x when type = l
..."right" way (of course, a personal opinion). Hence a feature request
to add a sort option when type = "l".
Before giving a reproducible example, is there anything even simpler
that I have missed?
Example:
Everything included:
> source ('http://www.stat.osu.edu/~pkapat/miscl/Code4xyplot_sort_type_l.R')
OR, step by step:
> library (lattice)
>
> # toy data
> D <- read.csv ('http://www.stat.osu.edu/~pkapat/miscl/Data4xyplot_sort_type_l.csv')
>
> # default behavior: obviously wrong
> xyp1 <- xyplot (y ~ x/1000 | C, groups = G, data...
2005 Jul 31
1
eigen gives NaN in $vectors (PR#8041)
....1.1 (2005-06-20)
OS: GNU/Linux 2.6.8-2-386, Debian testing
Submission from: (NULL) (65.24.56.41)
Relevant Bugs Ids : 7987, 7989
H is a 100x100 singular but symmetric matrix (a matrix defining the
neighbourhood structure for a spatial data) available from..
http://www.stat.ohio-state.edu/~pkapat/miscl/bugMatrix.RData
> load(file="bugMatrix.RData")
> ls()
"H"
> I <- diag(rep(1,100))
> g <- 0.1318467
> M <- I - g*H
> y <- eigen(M)
> y$vec[,91]
## all NaNs
Interestingly, as suggested in bug id 7989, this is very sensitive to the
scaling factor...
2010 Apr 30
1
Trouble using Ecdf () from the Hmisc library
...o Ecdf ()?
The details:
I've used Ecdf () _a lot_ over the past few years and I have learned
to live with its warnings. But I am running short on time and patience
now [*] Here is a reproducible example:
> library (Hmisc)
> x <- read.csv ( file = 'http://www.stat.osu.edu/~pkapat/miscl/Data4Ecdf.csv')
> table (x[, c('X1','X2')])
X2
X1 A B
6 0 1
13 7 20
15 13 26
16 13 13
17 10 19
18 12 10
As you can see (X1=6,X2=A) is a degenerate cell and (X1=6,X2=B) is
"near" degenerate! So, errors/warnings are certainly expected while
de...