[R] legend and lend (line end style)

2007-06-22 Thread Camarda, Carlo Giovanni
], pch=15, cex=2.5, col=1) points(lege$rect$left+0.03,lege$text$y[1], pch=15, cex=2.5, col=1) points(lege$rect$left+0.05,lege$text$y[1], pch=15, cex=2.5, col=1) points(lege$rect$left,lege$text$y[2], pch=15, cex=2.5, col=2) === Camarda Carlo Giovanni PhD

[R] Extracting from a matrix w/o for-loop

2006-07-28 Thread Camarda, Carlo Giovanni
=== Camarda Carlo Giovanni PhD-Student Max Planck Institute for Demographic Research Konrad-Zuse-Strasse 1 18057 Rostock, Germany Tel: +49 (0)381 2081 172 Fax: +49 (0)381 2081 472 [EMAIL PROTECTED] === -- This mail has been sent through the MPI

Re: [R] Help on plotting a 3 dimensional surface

2006-06-03 Thread Camarda, Carlo Giovanni
, a possible function would be plot.surface(fields). Ciao, Carlo Giovanni Camarda === Camarda Carlo Giovanni PhD-Student Max Planck Institute for Demographic Research Konrad-Zuse-Strasse 1 18057 Rostock, Germany Tel: +49 (0)381 2081 172 Fax: +49 (0)381 2081 472

[R] optim with specific constraints on possible values

2006-05-16 Thread Camarda, Carlo Giovanni
=== Camarda Carlo Giovanni PhD-Student Max Planck Institute for Demographic Research Konrad-Zuse-Strasse 1 18057 Rostock, Germany Tel: +49 (0)381 2081 172 Fax: +49 (0)381 2081 472 [EMAIL PROTECTED] === -- This mail has

Re: [R] draw a graphic using segments to represent data

2006-05-10 Thread Camarda, Carlo Giovanni
(Data), cex=2) axis(2) for(i in 1:nrow(Data)){ segments(x0=0, y0=Data[i,1], x1=Data[i,1] + Data[i,2], y1=Data[i,1], col=i) } Best, Carlo Giovani Camarda === Camarda Carlo Giovanni PhD-Student Max Planck Institute for Demographic Research Konrad-Zuse-Strasse

Re: [R] How to get a grid behind a boxplot

2006-04-29 Thread Camarda, Carlo Giovanni
=white) grid(ny=5) boxplot(a, col=red, add = TRUE) Ciao, Carlo Giovanni Camarda === Camarda Carlo Giovanni PhD-Student Max Planck Institute for Demographic Research Konrad-Zuse-Strasse 1 18057 Rostock, Germany Tel: +49 (0)381 2081 172 Fax: +49 (0)381 2081

[R] Different deviance residuals in a (similar?!?) glm example

2006-02-27 Thread Camarda, Carlo Giovanni
Dear R-users, I would like to show you a simple example that gives an overview of one of my current issue. Although my working setting implies a different parametric model (which cannot be framed in the glm), I guess that what I'll get from the following example it would help for the next steps.

[R] Changing labels on the scale of image.plot(fields)

2006-02-08 Thread Camarda, Carlo Giovanni
Dear R-Users, I am searching a way of changing the labels of the scale of a image.plot function in library(fields). Taking the following example I would like to add my own labels on the right side of the scale and it would be nice to decide also where to have the ticks. In this particular template

[R] Lexis maps in R

2006-02-08 Thread Camarda, Carlo Giovanni
Dear R-users, Sorry for bothering you twice in a day, but I was wondering whether there is any R-function which can easily plot the single elements of a Lexis diagram. I though that Lexis.diagram(Epi) could be the case, but it just plot life times in a frame. In particular I have been searching

[R] glm-logistic on discrete-time methods with individual and aggregated data

2006-02-01 Thread Camarda, Carlo Giovanni
Dear R-Users, without going into details I tried to prepare a simple example to show you where I would need help. In particular I prepare two examples-template for a study I'm conduction on discrete-time methods for survival analysis. Each of this example has two datasets which are basically

[R] Subsetting a matrix without for-loop

2006-01-30 Thread Camarda, Carlo Giovanni
Dear R-users, I'm struggling in R in order to squeeze a matrix without using a for-loop. Although my case is a bit more complex, the following example should help you to understand what I would like to do, but without the slow for-loop. Thanks in advance, Carlo Giovanni Camarda A - matrix(1:54,

[R] Solid and dotted lines saving an xyplot

2005-10-18 Thread Camarda, Carlo Giovanni
Dear R-users, I would a small problem saving a graph with postscript after using xyplot. Whereas in the R-screen the two lines are solid, once the image is saved as .eps file the second lines (blue in the following example) is dotted. Here there is a simple example. Thanks in advance for your

[R] Historical England and Wales Shape Files

2005-10-12 Thread Camarda, Carlo Giovanni
Dear R-user, I would like to apply spatial statistics on some historical data, in particular from 1850 to 1900 by registration district of England and Wales. I have searched in the R-archive, but unsuccessfully. Do you know whether R contains those shape files? Or where would be possible to

[R] Plotting rows (or columns) from a matrix in different graphs, not using par

2005-06-14 Thread Camarda, Carlo Giovanni
Dear R-users, I would like to ask whether it's possible (for sure it would be), to plot each rows (or columns) in different graphs and in the same figure region without using the function par and then struggling around with axes and labels etc. Luckily, I would always have rows + columns = even

[R] Save parameters from optim during iteration procedure

2005-05-15 Thread Camarda, Carlo Giovanni
[...] X = 0.0942275 -0.0942275 [...] X = -0.00110391 0.00110391 [...] X = 1.93949e-006 -1.93949e-006 [...] X = -4.98008e-013 4.98008e-013 [...] Thanks in advance, Carlo Giovanni Camarda Camarda Carlo Giovanni Max Planck Institute for Demographic Research Konrad-Zuse-Strasse 1 18057 Rostock

[R] persp, array and colors

2004-08-10 Thread Camarda, Carlo Giovanni
Dear R-users, I'd like to plot a three-dimensional surface and at the meantime I'm using an array. I would like to have the values of my first matrix in the heights of the plot and the colors of the single facet taking into account the second matrix. I hope that the next code will help all of you

[R] optimisation procedure with flat log-likelihood

2004-07-30 Thread Camarda, Carlo Giovanni
Dear R-friends, I use optim(par=c(mystartingpoints), fn=myloglikelihoodfunction, gr=NULL, method=c(L-BFGS-B), ## I would like to do not use any bounds control=list(trace=6, ## just to see what it's going on maxit=c(2)),

[R] histbackback for population pyramid

2004-04-22 Thread Camarda, Carlo Giovanni
Hi R-people, I hope someone could help me. I have to draw a population pyramid, but using histbackback (in Hmisc library) function I've got just histograms of density or frequency of my distribution. How can I get a population pyramid without playing on barplot function? I give a simple example