similar to: Making a case for using R in Academia

Greetings, I am trying to add a boxplot to the bottom of a histogram, right between the histogram bars and the x axis. Here is the code I am using at the moment (the par line is probably not relevant for our discussion): hs <- hist(x, breaks = 20, plot = F) par(mar = c(3,3,2,1)) hist(x, breaks = 20, main = NULL, ylim = c(-2, max(hs$counts))) boxplot(x, horizontal = T, axes = T, add =

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:

Dear R-Help, I'm new to R and struggling with weighting data when I run regression. I've tried to use search to solve my problem but haven't found anything helpful so far. I (successfully) import data from SPSS (15) and try to run a linear regression on a subset of my data file where WEIGHT is the name of my weighting variable (numeric), e.g.: library(foreign)

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 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

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 )

I ran into a weird, to me at least, problem, and hoping someone can shed some light into it. In a nutshell, there seems to be some problem when one calls plot with a formula, from within another function, using ... to pass arguments, and one of those arguments being xlim (and only xlim shows this problem). Here is an example: > plotw <- function(obj,...) { + plot(k~j,

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,

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

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 <-

Hi all, I'm trying to write a plugin for R Commander, following the model of the TeachingDemos plugin. I am struggling trying to even add items by editing the menus.txt file. I would welcome any help from anyone who has messed with it. Essentially the problem I am having is: I cannot create a submenu to a newly created menu. Here is what the menus.txt file looks like: # type

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

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 =

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

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

R users: I am trying to replicate the boxplot output I achieve with Minitab in R. I realize that R gives the user many more options on the algorithm used to calculate the IQR than Minitab, so I concentrated on type=6 when using the quantile() function in R. The problem I am having is setting the upper and lower limit of the whisker based on the nearest actual data that should be included. If

Hi, Below I have a function mlogl_k, later it's called with "nlm" . __BEGIN__ vsamples<- c(14.7, 18.8, 14, 15.9, 9.7, 12.8) mlogl_k <- function( k_func, x_func, theta_func, samp) { tot_mll <- 0 for (comp in 1:k_func) { curr_mll <- (- sum(dgamma(samp, shape = x_func, scale=theta_func, log = TRUE))) tot_mll <- tot_mll + curr_mll }

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

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 =