similar to: save a plot from R graphics window (in Windows)

I recently posted a question concerning an inconsistency of expand.grid in defining the reference level of the factors. > > expand.grid(x = c("b", "a"), y = c(1, 2))$x > [1] b a b a > Levels: b a # reference level is b > > expand.grid(x = c("b", "a"))$x > [1] b a > Levels: a b # reference level is a Thank you

Dear All, I'm trying to maximize a likelihood with respect one parameter using "optimize" on simulated data (without error component). I've iterated the maximization procedure 1000 times and I should always obtain the same estimate of the parameter (equal to the simulated one) but, instead, i obtain different results (the likelihood function shouldn't be flat). Does

I'm a PhD student and I'm working with covariance function. I'm interested to know if exist some packages in R to calculate and plot the bidimensional Autocovariance Function. the input matrix is a matrix that describe a spatial location over a 2-D space and I want to use it in the same way I can use a time serie in the 1-D acf. Thanks, Nicola.

I'm a PhD student and I'm working with covariance function. I'm interested to know if exist some packages in R to calculate and plot the bidimensional Autocovariance Function. the input matrix is a matrix that describe a spatial location over a 2-D space and I want to use it in the same way I can use a time serie in the 1-D acf. Thanks, Nicola.

Please can you help me I have the following problem: I have selected an lm model through the step procedure which visualize for each step the AIC value; then I have calculated for the initial model and the selected one the AIC using the funnction AIC. The results are different.What's happened? Emilia Rocco Dipartimento di Statistica "G. Parenti" Università di Firenze e-mail:

Hi I've always updated and committed my cvsroot directory using rsync also under windows with rsync distributed with cygwin. After installed the new release of cywing with rsync 2.5.6 protocol version 26, updating my cvsroot does not work anymore (under windows): it locks with receiving file list ... I use this command to update: rsync -e ssh -avzu --delete -c -C --include=CVS

Dear R-users, I found a strange problem working with products of two matrices, say: a <- A[i, ] ; crossprod(a) where i is a set of integers selecting rows. When i is empty the result is in a sense random. After some trials the right answer (a matrix of zeros) appears. --------------- Illustration -------------------- R : Copyright 2003, The R Development Core Team Version 1.8.0

Hi I would need to synschronize some Netscape stuff to a server ("roaming profiles" are not included anymore in the recent versions of Netscape). In particular I'd like to synchronize bookmarks, address book and mail filters. So I need to synchronize ONLY the following three files (that are in the mozilla directory work/mozilla) to the remote server (in the directory

I have a matrix 2900X2900 and I have to solve it. But R says that it can't allocate 67899kb. How can I do? Appreciate your kindly help! Salvati Nicola Dipartimento di Statistica "G.Parenti" salvati@ds.unifi.it tel. 055.4237.224 fax 055.4223560 [[alternative HTML version deleted]]

Hi All, I have the following adjacency matrix for a directed graph: [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 0 0 0 0 0 0 0 0 [2,] 0 0 0 0 0 0 0 0 [3,] 1 0 0 0 0 0 0 0 [4,] 0 0 1 0 0 0 0 0 [5,] 0 0 1 0 0 0 0 0 [6,] 1 1 0 0 0 0 0 0 [7,] 0 0

Dear R-users, after seeking for help in R-search I did not find any hint on my particular problem. Countless help on "true" overlay, but nothing on this. Please consider the following: par(mfrow = c(2, 1)) T <- seq(0, 20, by = 0.01) ## PLOT 1 plot(T, 30*exp(-0.65*T), type = "l" ) points(T, 30*exp(-0.26 * T), type = "l", lty = "F8") points(T,

Dear R-users, consider the two following outputs, ## 1 and ## 2 \begin{Scode}{Setup, echo = FALSE, print = FALSE, eval = TRUE} with(expand.grid(Fact1 = 1:3, Fact2 = 1:40), table(Fact1, Fact2)) ## 1 xtable(with(expand.grid(Fact1 = 1:3, Fact2 = 1:40), table(Fact1, Fact2))) ## 2 \end{Scode} The first line with(expand.grid(Fact1 = 1:3, Fact2 = 1:40), table(Fact1, Fact2))

Rakknar ha scritto: > Hello everybody. I've been learning R for about a month to do a > econometric study and now i'm stuck with some problems to make R do the > things I want. Here I give the list of things I wanna do from the most > simple to the more complex (for me of course): > > 1. Make a log. I've been using Stata and there i have a great tool to >

Dear R-users, I am trying to plot in the same map both NUTS 2 (i.e. regional) and NUTS 3 (i.e. provincial) boundaries for Italy but I couldn't find how to plot more than one level boundary in the same map.// // I have shapefile both for regions and provinces in Italy and I tried to use the "maptool" and "maps" packages without success. Thank you in advance for any help.

Hi there, this is my first post to the list. I'm a newcomer to the linux world. I started using it on a regular basis on May and I must say I'm quite comfortable with it, even if I have to re-learn a lot of things. But this is not a problem, I will improve my knowledge with time. My main problem now is that I bought myself a new computer (described at the end of this e-mail) and I was

Hello, I write to know how can I modify the x axis : when I plot a survival object, R plots a graph with x values = 0, 10, 20, 30 while I want a graph with values 0, 6, 12, 18, 24 in the x axis. How can I do this? In R 2.1.1 version there was "time.inc" in survplot, but in version R 2.2.0 there isn't it! I am sorry for my english and I hope that you understand my problem. Thank you

Dear R users, I have to call fortran program from within R (R 2.8.1 on ubuntu 8.10 machine). Suppose I have a fortran code like this (this is only a toy model, my working model is far more complex, but input/output is similar) DOUBLE PRECISION FUNCTION model(times, alfa, beta) DOUBLE PRECISION alfa, beta, times model=beta*sin(times)+alfa*cos(times) END FUNCTION which

Dear R users, I would like to exactly reproduce a figure like the 1.5 or 1.9 or 4.13 from the book Mixed effects models in S and S-Plus. Not for the sake of it, but because I have my own data I would like to plot in that fashion (no colors) If I write plot(ergoStool) I can get a good informative plot with colors, but I would like to have a B&W one instead. I've played a little with

Hi R users, I'm here trying to understand correlated residuals in nonlinear estimation. I'm reading/studying the book Bates, D. M. and D. G. Watts, (1988), /Nonlinear regression analysis and its applications/, Wiley, NY. pages 92-94, trying to reproduce the figures and to find out the code in R to perform the necessary calculations. I also consulted Pinheiro and Bates, but without

Hi all, first of all thanks for the answer to my previous question on lattice some time ago. In particular to Patrick Connolly for advices on netiquette (I hope this time I'm doing right....) and to Deepayan Sarkar fro the help on lattice. Now, my nowaday problem. Please consider the following mydf<-cbind.data.frame( RESPONSE = c(rnorm(9,rep(2:4,each=3),10),