search for: kapat

Full_Name: P Kapat Version: 2.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.e...

...m the "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,...

OK, here is the issue: I have a graphcis device (say X11) which I can set as active by dev.set(n). Is there ANY way to identify/infer whether this X11 device was obtained from a lattice function (like xyplot, levelplot...) or a standard function (like plot, scatterplot,...). A related question is whether plot.new (or frame) was called to obtain this graphics device... Thanks. PK

...tive to 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 wh...

...# som cylinder om h?jden ?r <3 m if(datamat$LangdHojd<30){ h?gstvol<- ( (datamat$Diameter/2) /100) ^2 * pi * (datamatris$LangdHojd/10) } # och som kapat trad om h?jden ?r h?gre else # f?rst r?knas ursprunglig h?jd ut urspr.tradh?jd<- exp(S?derberg(datamat, pars2))/10 # omr?knat till m #h?jdkvot i<- datama...

Hi, I have two matrices, A (axd) and B (bxd). I want to get another matrix C (axb) such that, C[i,j] is the Euclidean distance between the ith row of A and jth row of B. In general, I can say that C[i,j] = some.function (A[i,], B[j,]). What is the best method for doing so? (assume a < b) I have been doing some exploration myself: Consider the following function: get.f, in which,

Hi, Excuse me for using this list as a proxy for "r-sig-linux"! Is there anything analogous to `basename $0` (more importantly, "$0") of bash in R? Of course, this only makes sense in either BATCH mode or source() or shebang scripts (e.g using littler or Rscript)? Thanks for being considerate ;) -- Prasenjit

Hi, A. In a nutshell: The training error, obtained as "error (ret)", from the return value of a ksvm () call for a eps-svr model is (likely) being computed wrongly. "nu-svr" and "eps-bsvr" suffer from this as well. I am attaching three files: (1) ksvm.R from the the kernlab package, un-edited, (2) ksvm_eps-svr.txt: (for easier reading) containing only eps-svr