search for: charilao

...ee them next to each other, to develop a feeling for what histograms might correspond to what boxplots. Is there perhaps some reason why I should avoid showing those graphs to them like that, that I am not aware of? Or just a reason why I haven't seen them combined like this much? TIA Charilaos Skiadas Department of Mathematics Hanover College P.O.Box 108 Hanover, IN 47243

Hello all, I ran into the following, to me unexpected, behavior. I have (for reasons that don't necessarily pertain to the question at hand, hence I won't go into them) the need/desire to use an empty string for the name of a vector entry. Perhaps I did not read ?"[" very carefully, but it seems to me that he following lines should return "1" at the end:

Greetings, I have a problem that I am sure is very straightforward, but I just can't wrap my head around it. I've read the help pages on text, plotmath, expression, substitute, but somehow I can't find the answer to this simple question. Basically consider the following example: plot( NULL, xlim = c(0,2), ylim = c(0,2) ) expressions <- expression( -infinity, infinity )

...ing used outside of academia, or in general any other resources that talk about R as a practical alternative to the other non-free statistical packages. Perhaps some statistics, or particular examples of use? Any links would be greatly appreciated. Thanks for any thoughts/input into this. Charilaos Skiadas Department of Mathematics Hanover College P.O.Box 108 Hanover, IN 47243

hello all sorry for the following "none" R related question. does anyone know of a reference to calculate the following identity: |I + ABC| where I is an identity matrix and A, B,C may not have to be square matrices? you help will be greatly appreciated. H. V. Henderson; S. R. Searle SIAM Review, Vol. 23, No. 1. (Jan., 1981), pp. 53-60. provides a result to

I have a data frame named "database" with panel data, a little piece of which looks like this: Symbol Name Trial Factor1 Factor2 External 1 548140 A 1 -3.87 -0.32 0.01 2 547400 B 1 12.11 -0.68 0.40 3 547173 C 1

Hi, I'm looking for some way to pick up the numbers which are contained and buried in a long character. For example, outtree.new="(((B:1204.25,E:1204.25):7581.11,F:8785.36):8353.85,C:17139.21);" num.char =

...compelled to advertise it. Didn't see any previous postings about > it. Hope no offence taken! > > Best wishes, and well done with book. Looking forward to receiving > it... > > Mark > > 2008/5/2 Deepayan Sarkar <deepayan.sarkar at gmail.com>: >> On 4/30/08, Charilaos Skiadas <cskiadas at gmail.com> wrote: >> >>> Actually it's been out for a couple of weeks now at least. >> >> Yes, it's been out since March 12, actually. >> >> >>> I just finished >>> my first reading of it, and I must say i...

I just acquired a copy of "Statistical Models in S", I guess most commonly known as the "white book", and realized to my dismay that most of the code is not directly executable in R, and I was wondering if there was a source discussing the things that are different and what the new ways of calling things are. For instance, the first obstacle was the solder.balance data

I am getting some unexpected results from some functions of igraph and it is possible that I am misinterpreting the vertex numbers. Eg., the max betweenness measure seems to be from a vertex that is not connected to a single other vertex. Below if my code snippet: require(igraph) my.graph <- graph.adjacency(adjmatrix = my.adj.matrix, mode=c("undirected")) most.between.vert <-

After many hours of debugging code, I came to the conclusion that I have a fundamental misunderstanding regarding eval, and hope that someone here can explain to me, why the following code acts as it does: foo <- function(expr) { eval(substitute(expr), envir=list(a=5), enclos=parent.frame()) } bar <- function(er) { foo(er) } > foo(a) [1] 5 > bar(a) Error in eval(expr,

Several people have given great advice on how to successfully use X11 on Mac to so we can use the dvi previewer to view latex() output from Hmisc. Now after a version upgrade of X11 and X11sdk we are getting lots of crashes. I noticed a Mac dvi previewer with an executable stored as /Applications/TeX/TeXShop.app/Contents/MacOS/TeXShop Can we bypass X11 problems and point Hmisc to this

Hoping this is not off topic... I am in the process of writing some tutorials for my students for learning R, and naturally I'm using Sweave for this. So suddenly a question occurred to me: LaTeX has a recommended way of typesetting the TeX and LaTeX symbols, via the \TeX and \LaTeX commands. Is there a similar command for the R symbol, or in general are there any

I would like to get rid of the row numbers using xtable and latex. The commands d = cbind(1:10,rep(1:2,5)) ans = xtable(d) latex(ans) gives output containing \begin{tabular}{rrr} \hline & 1 & 2 \\ \hline 1 & 1.00 & 1.00 \\ 2 & 2.00 & 2.00 \\ 3 & 3.00 & 1.00 \\ 4 & 4.00 & 2.00 \\ 5 & 5.00 & 1.00 \\ 6 & 6.00 &

Dear all, I haven't seen this mentioned and presumably Dr. Deepayan Sarkar is too polite to advertise! I just received a flyer from Springer: A new book on Lattice released today!! http://www.springer.com/statistics/computational/book/978-0-387-75968-5?cm_mmc=NBA-_-Apr-08_UK_1753460-_-product-_-978-0-387-75968-5

I'm trying to do a linear regression between the columns of matrices. In example below I want to regress column 1 of matrix xdat with column1 of ydat and do a separate regression between the column 2s of each matrix. But the output I get seems to give correct slopes but incorrect intercepts and another set of slopes with value NA. How do I do this correctly? I'm after the slope and

I have the following relatively simple problem. Say we have three factors, and we want to create a cross-tabulation against each of the other two: x <- factor(rbinom(5, 1, 1/2)) y <- factor(rbinom(5, 1, 1/2)) z <- factor(rbinom(5, 1, 1/2)) table(x,y) table(x,z) This looks like: y x 0 1 0 2 0 1 1 2 z x 0 1 0 1 1 1 2 1 I would like to get (surely this will

Hello list, the following approach did not work: clustersA <- pam(distances, nkA, diss=TRUE); gc(); filenameclu = paste("filenameclu", ".txt"); write.table(clustersA , file=filenameclu,sep=","); although it worked with clustersA <- hclust(distances, method="ward"); and a consecutive kclassA <- cutree(clustersA, k=nkA); filename =

Hi everybody, I have to recognize that my post certainly shows my lack of skills in really navigating the R web page .... i am surprised that only one or two persons wrote me about the "Task Views" - it is what i was after - although maybe too general for my lazy taste - but hei - it is there and certainly makes my life easier. So - Thank you John and Gabor and Roger for pointing this

Hi, I can get from a string to a function with this name: >f1 <- function(x){ mean(x) } >do.call("f1",list{1:4}) >get("f1") etc... But how do I get from a function to its name? >funcVec <- c(f1,median) >funcVec [[1]] function(x){ mean(x) } > str(funcVec) List of 2 $ :function (x) ..- attr(*, "source")= chr "function(x){ mean(x)