similar to: bar plot patterns

Hi, I am plotting a policy function (result from a dynamic stochastic optimization problem, discretized approximation). The policy function maps from an 2 x 2 x 2 x 3 x B x F state space to a B x F state space (B and F are usually between 4-6, and represent domestic and foreign savings. The other variables are income (Y), inflation (Pi), domestic and foreign interest rates (R and Z)). I

Dear R-helpers, I will like to know if there is a way to generate a stacked column graph using both patterns and colors to fill the bars. I have many categories for the number of color available in R, so I will like R to start with solid colors and then use patterns an colors. I have been using barplot and some code that I found in the archives, but I am still having bars with the

Dear list, I'd like to use the function cut() on matrices, ie that when I apply it to a matrix, it would return a matrix of the same dimensions instead of a vector. I wonder if there is a better (more elegant) solution than matrix(cut(a, ...), ncol=ncol(a), nrow=nrow(a)) because I would like to use cut on both vectors and matrices and avoid testing whether a is a matrix. Thanks, Tamas

Dear list, I have tried to setup my Xemacs for use with Sweave, which I indend to learn. I have followed the instructions in the Sweave FAQ, that is to say, I put (defun Rnw-mode () (require 'ess-noweb) (noweb-mode) (if (fboundp 'R-mode) (setq noweb-default-code-mode 'R-mode))) (add-to-list 'auto-mode-alist '("\\.Rnw\\'" . Rnw-mode)) (add-to-list

Hi, I am about to but a laptop, and have narrowed the choices down to a Dell Latitude 600 and an Apple Powerbook G4 Aluminium (Princeton provides these models at a discount for grad students). I am biased towards the Powerbook, and would like to run Debian on it. I have only used debian on i386 platforms so far. I use R quite frequently, so I would be interested in your experience of running R

Janusz Kawczak wrote: > You simply need to remove the stuff related to MS Win from zzz.R; > in partricular the lines after if( .... ) to clear your message. > As you can see, the info relates to the WinMenu under MS Win. I think people have been more than a little disingenuous in claiming that getting the Rmetrics package to go under Linux is transparent. If you have to dig into the

I need to find a computationally simple process for the movement of interest rates. In this simplified model, an interest rate can have 3--5 possible values, and its movement is characterized by a matrix of transition probabilities (ie, it is a Markov process). I would like to estimate this process from a given set of data. For example, let the interest rate time series be: 7 3 8 2 5 9 6

Dear List, I have a couple of (~200) 3x3 transition probability matrices (ie each defines a Markov chain). They are all estimated from the same underlying process, so it ie meaningful to take their elemetwise mean and standard deviation. [1] First question: supposing that they are given in a list l, how do I get their elementwise mean and standard deviation? Fortunately, the mean of trans. prob.

Hi, The department of economics at our university (Budapest) is planning a course on numerical methods in economics. They are trying to decide which software to use for that, and I would like to advocate R. The other alternative is Matlab. I have found comparisons in terms of computational time for matrix algebra, but I don't think that is relevant: the bottleneck for economists is usually

Hi, I am typesetting R code in TeX using the listings package. I have experimented with various settings for the text, but all look a bit ugly. This might be because I have no typographic experience ;-) I would really appreciate if people sent me the settings they use for the listings package (eg in the \usepackage line, or \lstset, I am thinking about choices for basicstyle, ... etc). I

Suppose I have A, an n x m matrix, each element is an integer (an index). I also have B, an n x l x m array. I need C, where C[n,m] = B[n, A[n, m], m] I am currently using loops, what would be the "R way" to do this? Another question: let A[n, m] <- argmax_l B[n, l, m] what would be the nicest way of doing this? Currently I am using max.col and a single loop, going though the

Hi, Sometimes I want to "reset" the R interpreter (as if I quit and started it again, with everything starting from scratch). The main reason for this is when I expreriment with calculations, and later form them into general functions, sometimes I forget to include variables or reset some settings which I don't notice because they are already set. This causes problems during batch

Is there a way in R to test if a given number is an integer, ie a whole number? I am not referring to the data type of a number, but to its value. That is to say, is.whole(pi-pi+2) would be TRUE, whereas is.whole(4/3) would be false. At the moment I am using is.whole <- function(a) { floor(a)==a } which is OK for real numbers, but not for complex ones (a+bi would be a whole number if both a

Hi, I have a dynamic programming problem with many state variables, let's call them l, n, m, etc. The transition probabilities are originally given in an array form, eg transtition[l,m,n,ll,mm,nn] give the probability of going from l,m,n to ll,mm,nn. However, the numerical solution is best handled when I "flatten" the L x M x N state space into a single one (call it S), ie a

Hello, In S-Plus Windows you can transform graphics to Powerpoint very easily, in R Windows you can use enhanced metafiles (.emf) and Powerpoint almost as easy. Is there a simular way with R in Linux to transform to the presentation program in StarOffice or OpenOffice or are you stuck with the pdf device? Fredrik Lundgren

Hi, I need to parse a data file (output of a measuring device) of the following format: BEGIN RECORD [first record data] RECORD [second record data] RECORD [third record data] END Line breaks can (and do ;-() occur anywhere. White space behaves very much like TeX, eg it is not important whether there are one or more spaces or linebreaks as long as there is one of them. It is a text file, not

Hi, I have the following problem. It has two parts. 1. I need to calculate the stationary probabilities of a Markov chain, eg if the transition matrix is P, I need x such that xP = x in other words, the left eigenvectors of P which have an eigenvalue of one. Currently I am using eigen(t(P)) and then pick out the vectors I need. However, this seems to be an overkill (I only need a single

R: I have a program that runs a For loop for days and I need to (if possible) pause the program. Any ideas on how to do that? If I use stop, how certain can I be that all the statements executed in the previous iteration of the for loop? Thanks, Jason Higbee Research Associate Federal Reserve Bank of St. Louis E: jason.l.higbee@stls.frb.org [[alternative HTML version deleted]]

R 1.8.1, Windows 2000. I am trying to find the "legend" for color coding and point types in a plot, which probably are standard for everybody but myself. Thanks for the tip! Mihai Nica Jackson State University "No good deed will ever remain unpunished" [[alternative HTML version deleted]]

Rational functions (ratios of polynomials) often provide good approximations to many functions. Does anyone know of any literature on nonlinear modeling with rational functions, sequential estimation, diagnostics, etc.? I know I can do it with "nls" and other nonlinear regression functions, but I'm wondering what literature might exist discussing how a search for an