[R] problems with too many NA in the function ideal() from pscl package.

2008-09-19 Thread Antonio P. Ramos
Hi all, I'm trying to run some monte carlo simulation for my roll call data using the ideal() function, which resides in the pscl package. However, I'm receiving an error message that I don't understand. Error in ideal(a, maxiter = 1000, thin = 10, burnin = 50, store.item = TRUE, : NA/NaN/Inf

Re: [R] Mixed effects model with binomial errors - problem

2008-09-19 Thread RFTW
anyone? RFTW wrote: ok... the model now runs properly (say, without errors). Now about the result. These are the averages per treatments tapply(VecesArbolCo.VecesCo.C1,T2,mean) a b c d 0.49 0.56 0.45 0.58 I run this very simple model

Re: [R] help on sampling from the truncated normal/gamma distribution on the far end (probability is very low)

2008-09-19 Thread Peter Dalgaard
Daniel Davis wrote: Hi, guys, I am trying to sample from a truncated normal/gamma distribution. But only the far end of the distribution (where the probability is very low) is left. e.g. mu = - 4; sigma = 0.1; The distribution is Normal(mu,sigma^2) truncated on [0,+Inf]; How can I get a

Re: [R] how to keep up with R?

2008-09-19 Thread Robin Hankin
Hi Wensei. Why not do as I do? Find an interesting area of numerical computation (perhaps not statistical) that has not been implemented in open-source. Then write an R package for it, under GPL-2, then write an article about the new package in R-news or JSS. works for me. Best wishes

Re: [R] How to do knn regression?

2008-09-19 Thread Hans W. Borchers
Shengqiao Li shli at stat.wvu.edu writes: Hello, I want to do regression or missing value imputation by knn. I searched r-help mailing list. This question was asked in 2005. ksmooth and loess were recommended. But my case is different. I have many predictors (p20) and I really want try

Re: [R] how to keep up with R?

2008-09-19 Thread Barry Rowlingson
2008/9/19 Wensui Liu [EMAIL PROTECTED]: Dear Listers, I've been a big fan of R since graduate school. After working in the industry for years, I haven't had many opportunities to use R and am mainly using SAS. However, I am still forcing myself really hard to stay close to R by reading

[R] Type I SS and Type III SS problem

2008-09-19 Thread leeznar
Dear all: I m a newer on R.  I have some problem when I use anova function.  I use anova function to get Type I SS results, but I also need to get Type III SS results.  However, in my code, there is some different between the result of Type I SS and Type III SS.  I don’t know why the “seqe”

Re: [R] How to find a shift between two curves or data sets

2008-09-19 Thread Hans W. Borchers
Dear Sébastien, identifying similarity in curves or time series is one of the main tasks in the quite recent field of 'Functional Data analysis' (FDA). See the 2005 book by Silverman from Springer Verlag or the corresponding Web page at url{http://www.psych.mcgill.ca/misc/fda/}. The 'fda'

[R] frequency table across multiple variables

2008-09-19 Thread Ralikwen
Dear R users, I have a dataframe like this: x1-c(1,2,3,4,NA ,NA ,NA, 3, 1, 1, 1, 1, 2, 2, 3, 4, 4) x2-c(2,3,4,3,4,3,4,2,2,3,4,NA,NA,NA,NA,4,3) x3-c(1,1,1,1,2,2,2,2,3,3,3,3,4,4,4,1,2) m-data.frame(x1,x2,x3) I would like to create a frequency table like this: x1 x2 x3 NA 1 2 3 4 where

Re: [R] PDF fonts problem

2008-09-19 Thread Mihalicza Péter
Dear Dr. Murrel, Thank you for all the clarifications! Paul Murrell írta: Hi #CMS pdf(tryfont-cms.pdf, family=CMS) grid.text(gg\u151hh\uF6ii\uF3jj kk\u171ll\uFCmm\uFAnn) dev.off() #u151 and u171 doesn't show, though the other accented ones do embedFonts(tryfont-cms.pdf,

Re: [R] Mixed effects model with binomial errors - problem

2008-09-19 Thread ONKELINX, Thierry
First of all I'm forwarding this mail to the R-SIG-mixed-models, which is more appropriate to your question. Remember that family = bionomial uses by default the logit link. Hence all parameters will be on the logit scale. So you will need to backtransform them for comparison. Then you'll see

Re: [R] frequency table across multiple variables

2008-09-19 Thread Philipp Pagel
I have a dataframe like this: x1-c(1,2,3,4,NA ,NA ,NA, 3, 1, 1, 1, 1, 2, 2, 3, 4, 4) x2-c(2,3,4,3,4,3,4,2,2,3,4,NA,NA,NA,NA,4,3) x3-c(1,1,1,1,2,2,2,2,3,3,3,3,4,4,4,1,2) m-data.frame(x1,x2,x3) I would like to create a frequency table like this: x1 x2 x3 NA 1 2 3 4 where

Re: [R] Type I SS and Type III SS problem

2008-09-19 Thread Peter Dalgaard
leeznar wrote: Dear all: I m a newer on R. I have some problem when I use anova function. I use anova function to get Type I SS results, but I also need to get Type III SS results. However, in my code, there is some different between the result of Type I SS and Type III SS. I don’t

Re: [R] frequency table across multiple variables

2008-09-19 Thread Ivar Herfindal
Hi Here is an alternative solution which will include count of the NA. It may not be formally correct, but it seems to work: mapply(function(x) table(factor(ifelse(is.na(x), NA, x), levels=c(NA,1,2,3,4))), m) x1 x2 x3 NA 3 4 0 1 5 0 5 2 3 3 5 3 3 5 4 4 3 5 3 Hope

Re: [R] Use of colour in plots

2008-09-19 Thread ONKELINX, Thierry
Steve, - Use tranparancy to prevent overplotting: more details on p. 16 of the ggplot2 book: http://had.co.nz/ggplot2/book/ - You can choose your own colour with scale_manual(): http://had.co.nz/ggplot2/scale_manual.html - The backgroundcolor can be set with ggopt(background.color = white):

[R] panel data analysis possible with mle2 (bbmle)?

2008-09-19 Thread Erich STRIESSNIG
Dear R community, I want to estimate coefficients in a (non-linear) system of equations using 'mle2' from the bbmle package. Right now the whole data is read in as just one long time series, when it's actually 9 cross sections with 30 observations each. I would like to be able to test and correct

[R] Error message in lmer

2008-09-19 Thread Thomas Sattler
Dear list I try to run a bootstrap with lmer. I often, but not always, get the error message: Error in objective(.par, ...) : Leading minor of order 6 in downdated X'X is not positive definite (with the number (here 6) varying) In R-archives I came across some threads that treated this

Re: [R] Still no R-2.7.2 rpms for Redhat Enterprise Linux

2008-09-19 Thread Martyn Plummer
We'll see what we can do. I have forwarded your email to Bob Kinney, who builds the RHEL RPMS. Martyn On Thu, 2008-09-18 at 14:57 -0700, Waichler, Scott R wrote: There's already an announcement about R 2.8.0, and yet there are still no R-2.7.2 binaries (rpms) for Redhat Enterprise Linux 4 and

Re: [R] frequency table across multiple variables

2008-09-19 Thread Ralikwen
Hi, I went for a slight alteration of your solution x1-c(1,2,3,4,NA ,NA ,NA, 3, 1, 1, 1, 1, 2, 2, 3, 4, 4) x2-c(2,3,4,3,4,3,4,2,2,3,4,NA,NA,NA,NA,4,3) x3-c(1,1,1,1,aaa,2,2,2,3,3,3,3,4,4,4,1,2) m-data.frame(x1,x2,x3) m-replace(m,is.na(m),NA) levels=unique(as.vector(as.matrix(m)))

[R] migrating data from s-plus to R

2008-09-19 Thread Patrizio Frederic
Dear all, is there any way to transform a .Data directory created in S-plus 6.1 for windows in a .RData file? Thanks in advance, Patrizio Frederic __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

Re: [R] frequency table across multiple variables

2008-09-19 Thread Kenn Konstabel
or ... sapply(m,function(x) table(factor(x, levels=c(NA, 1:4), exclude=NULL))) On Fri, Sep 19, 2008 at 12:59 PM, Ralikwen [EMAIL PROTECTED] wrote: Hi, I went for a slight alteration of your solution x1-c(1,2,3,4,NA ,NA ,NA, 3, 1, 1, 1, 1, 2, 2, 3, 4, 4)

Re: [R] R on a eeepc 901 ?

2008-09-19 Thread Agustin Lobo
The one with Linux (Xandros) As an alternative in case I don't get through with the R installation(s) in Xandros, I've installed ubuntu-eee on a usb memory stick for booting from it and plan to go on with a regular ubuntu R installation, but have not decided yet whether I'll be erasing the

[R] R: Adding 1 month to a dataframe column

2008-09-19 Thread ANGELO.LINARDI
It works great! Thank you so much Gabor Angelo Linardi -Messaggio originale- Da: Gabor Grothendieck [mailto:[EMAIL PROTECTED] Inviato: giovedì 18 settembre 2008 21.16 A: LINARDI ANGELO Cc: r-help@r-project.org Oggetto: Re: [R] Adding 1 month to a dataframe column Try this: # example

[R] spatstat - tabulate point pattern

2008-09-19 Thread Paulo Cardoso
Hi, Is it possible to tabulate (count) points generated by spatstat ppp() using any regular Grid created with owin()? The idea is to count No of points in each grid cell. Any other method available? Paulo __ R-help@r-project.org mailing list

Re: [R] migrating data from s-plus to R

2008-09-19 Thread John Kane
You might find something about this in the R manual R Data Import/Export --- On Fri, 9/19/08, Patrizio Frederic [EMAIL PROTECTED] wrote: From: Patrizio Frederic [EMAIL PROTECTED] Subject: [R] migrating data from s-plus to R To: r-help@r-project.org r-help@r-project.org Received: Friday,

Re: [R] How to show complete time values in a plot x axis

2008-09-19 Thread Sébastien Durand
Dear Yihui, The resulting plot produced by dotchart is not what I wish for! I need to produce a time axis that is consequent to the time values. When dotchart is used, time values are simply considered as labels. What I wish for is a way to produce a numerical time axis that really takes into

Re: [R] How to show complete time values in a plot x axis

2008-09-19 Thread Duncan Murdoch
On 9/19/2008 9:14 AM, Sébastien Durand wrote: Dear Yihui, The resulting plot produced by dotchart is not what I wish for! I need to produce a time axis that is consequent to the time values. When dotchart is used, time values are simply considered as labels. What I wish for is a way to produce

[R] intToUtf8

2008-09-19 Thread Christian Hennig
Hi there, any explanation for this? intToUtf8(66) Error in intToUtf8(66) : argument 'x' must be an integer vector intToUtf8(c(66,55)) Error in intToUtf8(c(66, 55)) : argument 'x' must be an integer vector intToUtf8(c(66,55),multiple=TRUE) Error in intToUtf8(c(66, 55)) : argument 'x' must

Re: [R] Use of colour in plots

2008-09-19 Thread Steve Murray
Jim, Thanks for this - I've looked into cluster.overplot in particular which, judging by the help file, sounds quite useful (count.overplot seems less relevant). I'm finding however, that when I execute cluster.overplot, it simply returns many values (which total the number in the dataset

[R] getting line breaks with xtable

2008-09-19 Thread Erich Studerus
Sorry, for asking the same question again, but I got no reactions the last time. Maybe it was just overseen by the experts. I'm using the xtable function with Sweave and Lyx and I would like to know how to get automatic line breaks for long strings in a column of the table. I've learned from the

Re: [R] Use of colour in plots

2008-09-19 Thread Steve Murray
Thierry, Thanks - I've had a look into using the transparency option, but can't seem to work out where to place it within the command I'm using: ggplot(Jan, aes(x = PopDensity, y = Average.Burnt.Area.Fraction, colour = factor(Urban.Rural 1.25))) + geom_point() I'm assuming that it has to

[R] Plot availability

2008-09-19 Thread john crepezzi
Is it possible to see if a plot is already open before I call lines()? -- View this message in context: http://www.nabble.com/Plot-availability-tp19572437p19572437.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org

[R] misalignement of secondary axis in multhist

2008-09-19 Thread Angel Marley
Hi everybody, I'm trying a multiple histogram, and I'd like to add on it a secondary axis with a line (cumulative frequencies, or wathever lines), and the 0 intercept doesnt coincide in y1 and y2, which dificult interpretation and makes presentation not nice. I had seen some related answer for

Re: [R] detecting null values in a CSV file

2008-09-19 Thread Jason Thibodeau
On a related note, I am trying to do some matching using conditional statements. These NULL values are being brought in to my data frame as NA, as expected, but in a conditional if() statement, I cannot compare then to a integer value, it fails the program. Here is a small snippet of where the

Re: [R] getting line breaks with xtable

2008-09-19 Thread Gabor Grothendieck
Read the last line to every message to r-help to find out one reason you may be getting no responses. On Fri, Sep 19, 2008 at 9:39 AM, Erich Studerus [EMAIL PROTECTED] wrote: Sorry, for asking the same question again, but I got no reactions the last time. Maybe it was just overseen by the

Re: [R] intToUtf8

2008-09-19 Thread Duncan Murdoch
On 9/19/2008 9:32 AM, Christian Hennig wrote: Hi there, any explanation for this? intToUtf8(66) Error in intToUtf8(66) : argument 'x' must be an integer vector intToUtf8(c(66,55)) Error in intToUtf8(c(66, 55)) : argument 'x' must be an integer vector intToUtf8(c(66,55),multiple=TRUE)

Re: [R] intToUtf8

2008-09-19 Thread Christian Hennig
Dear Duncan, No, they're numeric. integer here is a type, not the mathematical thing. This works: intToUtf8(as.integer(c(66,55))) The docs don't mention this requirement, and it does seem somewhat unnecessary; I'll look into it. Thanks! You're lucky I didn't read to the end of your

Re: [R] Use of colour in plots

2008-09-19 Thread Steve Murray
Sorry - I should've maybe also pointed out that the command I've been trying to use is: alpha(col=green, 1/10) On its own this results in the following error: [1] #00FF001A and I haven't been able to successfully incorporate it into the main formula just yet (please see my previous message).

Re: [R] Plot availability

2008-09-19 Thread Duncan Murdoch
On 9/19/2008 9:32 AM, john crepezzi wrote: Is it possible to see if a plot is already open before I call lines()? I don't think so. You can see if a graphics device is open by looking at dev.cur(), but I don't think there's a test for plot.new(). I'd just wrap the call in try() if you're

Re: [R] getting line breaks with xtable

2008-09-19 Thread Erich Studerus
Ok. Sorry, here's a reproducible example: library(xtable) x-as.table(cbind(1:3,rep(this is an example for a long character string that I want break into several lines))) xtable(x) Regards Erich -Ursprüngliche Nachricht- Von: Gabor Grothendieck [mailto:[EMAIL PROTECTED] Gesendet:

Re: [R] Use of colour in plots

2008-09-19 Thread ONKELINX, Thierry
Steve, You want something like this: library(ggplot2) n - 1000 dataset - data.frame(x = round(rnorm(n), 2), y = round(rnorm(n), 1), z = rnorm(n)) ggplot(dataset, aes(x = x, y = y, colour = factor(z 1))) + geom_point() + scale_colour_manual(values = c(alpha(red, 1/4), green)) HTH, Thierry

Re: [R] intToUtf8

2008-09-19 Thread Duncan Murdoch
On 9/19/2008 9:58 AM, Christian Hennig wrote: Dear Duncan, No, they're numeric. integer here is a type, not the mathematical thing. This works: intToUtf8(as.integer(c(66,55))) The docs don't mention this requirement, and it does seem somewhat unnecessary; I'll look into it. Thanks!

Re: [R] How to do knn regression?

2008-09-19 Thread Shengqiao Li
Date: Fri, 19 Sep 2008 07:00:33 + (UTC) From: Hans W. Borchers [EMAIL PROTECTED] Subject: Re: [R] How to do knn regression? To: [EMAIL PROTECTED] Message-ID: [EMAIL PROTECTED] Content-Type: text/plain; charset=us-ascii Shengqiao Li shli at stat.wvu.edu writes: Hello, I want to do

Re: [R] how to keep up with R?

2008-09-19 Thread Richard . Cotton
The good thing about training and evangelization ... I spent quite a lot of time evangelizing about R when I first started my current job. Eventually my boss told me I was an R-Soul, or something that sounded like that anyway. ; ) Regards, Richie. Mathematical Sciences Unit HSL

Re: [R] How to show complete time values in a plot x axis

2008-09-19 Thread Yihui Xie
Sorry, I thought the time values were equal-spaced... Regards, Yihui -- Yihui Xie [EMAIL PROTECTED] Phone: +86-(0)10-82509086 Fax: +86-(0)10-82509086 Mobile: +86-15810805877 Homepage: http://www.yihui.name School of Statistics, Room 1037, Mingde Main Building, Renmin University of China, Beijing,

[R] plot order in multi-panel figure

2008-09-19 Thread Stephen Tucker
Hi, Does anyone know if there is a way to 'reset the plot number' on a traditional graphics device? For instance, I want to have two plots on stacked top of each other (mfrow=c(2,1)) but with underlying grid lines spanning both figures vertically. I can put the grid lines on top if I add

Re: [R] getting line breaks with xtable

2008-09-19 Thread Gabor Grothendieck
Try placing your text in a parbox: x-as.table(cbind(1:3,rep(\\parbox{5cm}{this is an example for a long character string that I want break into several lines}))) print(xtable(x), sanitize.text = force) You may need to fix up the vertical spacing as well. On Fri, Sep 19, 2008 at 10:11 AM, Erich

Re: [R] getting line breaks with xtable

2008-09-19 Thread Marc Schwartz
I don't believe that \linebreak will work within a tabular environment as you expect here. Plus, be aware, that it is print.xtable() and not LaTeX that is adding the additional $\backslash$ characters, as part of the text sanitization process. I think that you basically have two options: 1.

Re: [R] plot order in multi-panel figure

2008-09-19 Thread Mark Lyman
Stephen Tucker brown_emu at yahoo.com writes: Hi, Does anyone know if there is a way to 'reset the plot number' on a traditional graphics device? For instance, I want to have two plots on stacked top of each other (mfrow=c (2,1)) but with underlying grid lines spanning both figures

Re: [R] Multiple logical operations in a subscript

2008-09-19 Thread Greg Snow
Peter showed you the %in% operator, you may also want to look at the subset, transform, with, and within functions for future use as ways to reduce the need to type the name of an object multiple times. Hope this helps, -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain

[R] Lines between panels in lattice

2008-09-19 Thread Saptarshi Guha
Hello, I have a multi-page display each consisting of two-panels above each other. I need to draw a line from the top panel to bottom panel. Using current.vpTree() i find that plot1.panel.1.2.vp and plot1.panel.1.1.vp are the top and bottom ones respectively. I am using the following

Re: [R] Use of colour in plots

2008-09-19 Thread Steve Murray
Hi Thierry and all, Thanks very much for your suggestion. I've given it a go and played around with the transparency values, but seem to be having a problem in that some of the red values are made transparent, even though there are no green values being overplotted! The code I used to

[R] ggplot2 - accessing legend data

2008-09-19 Thread Pedro Barros
Hi All, I am trying to build a set of plots, and build a common legend to all of them. Because I want to use a multi-line, multi-column legend layout, I cannot use the common ggplot2 legends (if I can, please tell me how). So, I would like to know how can I access the properties of a legend, so

Re: [R] Multiple logical operations in a subscript

2008-09-19 Thread Ted Harding
On 19-Sep-08 16:03:45, Greg Snow wrote: Peter showed you the %in% operator, you may also want to look at the subset, transform, with, and within functions for future use as ways to reduce the need to type the name of an object multiple times. Hope this helps, Let me add my perhaps very

[R] evaluate return values of commands executed with system()

2008-09-19 Thread Rainer M Krug
Hi I would like to check the return value of commands executed with system() I am using Linux (Ubuntu Hardy) and I am executing GRASS commands and would like to call stop() when they return an error (i.e. return code = 1). Thanks Rainer -- Rainer M. Krug, PhD (Conservation Ecology, SUN),

[R] Novice question about getting data into R

2008-09-19 Thread Ted Byers
I found it easy to use R when typing data manually into it. Now I need to read data from a file, and I get the following errors: refdata = read.table(K:\\MerchantData\\RiskModel\\refund_distribution.csv, header = TRUE) Error in scan(file, what, nmax, sep, dec, quote, skip, nlines,

[R] Cbind help

2008-09-19 Thread Rajasekaramya
Hi there I want to cbind a vector(contains 20 values) to a dataframe containg (18 rows).The missing value in the dataframe is however a repeated values but still i want to retain the repeats in the vector and the dataframe values can even get enterd twice. how should i go about it Example:

Re: [R] Novice question about getting data into R

2008-09-19 Thread Peter Dalgaard
Ted Byers wrote: I found it easy to use R when typing data manually into it. Now I need to read data from a file, and I get the following errors: refdata = read.table(K:\\MerchantData\\RiskModel\\refund_distribution.csv, header = TRUE) Error in scan(file, what, nmax, sep, dec,

Re: [R] Novice question about getting data into R

2008-09-19 Thread John Kane
Try read.csv(K:\\MerchantData\\RiskModel\\refund_distribution.csv,header = TRUE) --- On Fri, 9/19/08, Ted Byers [EMAIL PROTECTED] wrote: From: Ted Byers [EMAIL PROTECTED] Subject: [R] Novice question about getting data into R To: r-help@r-project.org Received: Friday, September 19, 2008,

Re: [R] Lines between panels in lattice

2008-09-19 Thread Saptarshi Guha
Got it. It should be plot1.panel.1.2.off.vp and likewise. Thanks Saptarshi On Sep 19, 2008, at 12:08 PM, Saptarshi Guha wrote: Hello, I have a multi-page display each consisting of two-panels above each other. I need to draw a line from the top panel to bottom panel. Using

[R] reproduce this graph in ggplot2 (code and data included)

2008-09-19 Thread Juliet Hannah
How can I reproduce this graph in ggplot2 (regression lines and data point superimposed). Thanks, Juliet filename=http://personality-project.org/r/datasets/heating.txt; heating=read.table(filename,header=TRUE) symb=c(19,25,3,23) colors=c(black,red,green,blue)

Re: [R] reproduce this graph in ggplot2 (code and data included)

2008-09-19 Thread hadley wickham
Here's one way: qplot(degreedays, therms, data=heating, colour = Location, shape = Location) + geom_smooth(method = lm, fullrange = T, se = F) Hadley On Fri, Sep 19, 2008 at 12:45 PM, Juliet Hannah [EMAIL PROTECTED] wrote: How can I reproduce this graph in ggplot2 (regression lines and data

Re: [R] Coloring spaces between lines in xyplot

2008-09-19 Thread hadley wickham
On Thu, Sep 18, 2008 at 5:14 PM, Seth W Bigelow [EMAIL PROTECTED] wrote: Greetings: I wish to create a stacked line graph in xyplot, adding color to the spaces between the lines. For example, the code below creates a plot with two lines extending across it, and I want to color the rhomboid

Re: [R] R on a eeepc 901 ?

2008-09-19 Thread Agustin Lobo
Sigbert, When I do: apt-get install ksmserver kicker I get a message stating that there is no header in the package (Xandros came in Spanish, so I have to translate the message, which is not easy) Aterwards, I have to Edit sources.list file and create a preferences file. As the repositories

Re: [R] Novice question about getting data into R

2008-09-19 Thread Silvia Lomascolo
refdata = read.table(K:\\MerchantData\\RiskModel\\refund_distribution.csv, header = TRUE) Error in scan(file, what, nmax, sep, dec, quote, skip, nlines, na.strings, : line 1 did not have 42 elements refdata = read.table(K:\\MerchantData\\RiskModel\\refund_distribution.csv) Error in

Re: [R] how to keep up with R?

2008-09-19 Thread Adaikalavan Ramasamy
I agree! The best way to learn (and remember for longer) is to teach someone else about it. And there is not reason not to repeat some of the anlysis done on SAS with R. That way you can verify your outputs or compare the presentations. If you consistently find differences in the outputs,

[R] Plot Accessibility

2008-09-19 Thread john crepezzi
I'm not so sure the original way I stated my issue was clear enough, so I'll attempt to elaborate a little bit. I'd like to make a function that is passed the name of a plot object, and a lines/point specification, and graphs them all on the same plot. I don't want to redraw the plot each time,

Re: [R] Plot availability

2008-09-19 Thread john crepezzi
I'm not so sure the original way I stated my issue was clear enough, so I'll attempt to elaborate a little bit. I'd like to make a function that is passed the name of a plot object, and a lines/point specification, and graphs them all on the same plot. I don't want to redraw the plot each time,

Re: [R] Coloring spaces between lines in xyplot

2008-09-19 Thread Seth W Bigelow
Ok, thanks to Carl Witthoft I now know that to color spaces between and below a pair of lines in xyplot, I will need to redefine the lines as polygons, and use the lpolygon panel function, as in the following library(lattice) p - c(1,10,10,1) # vector of x values q - c(4,29.2,16,2.5)

[R] Weighted polytomous logistic regression

2008-09-19 Thread Yinghui Wang
Hi, I am trying to fit a polytomous logistic regression with inverted sampling probabity as weight. I've tried multinom(nnet) function. While the results look not right. For example, the results were not the same for the following two model. Should I normalize the weights first? Or there were

[R] Position WINDOWS TKTOPLEVEL

2008-09-19 Thread ermimi
Hello friends!!! I´m spanish student, I am finishing the finish project in computer science. I am creating a R package about CircularStatistics, In october i will upload it to CRAN. I need to resolve one problem. I need to put a windows created with tktoplevel into screen center. How I could

Re: [R] Novice question about getting data into R

2008-09-19 Thread Tom Backer Johnsen
Silvia Lomascolo wrote: refdata = read.table(K:\\MerchantData\\RiskModel\\refund_distribution.csv, header = TRUE) Error in scan(file, what, nmax, sep, dec, quote, skip, nlines, na.strings, : line 1 did not have 42 elements refdata =

[R] Fwd: Lines between panels in lattice

2008-09-19 Thread Saptarshi Guha
Forgot to send to the list. Begin forwarded message: From: Saptarshi Guha [EMAIL PROTECTED] Date: September 19, 2008 1:43:50 PM GMT-04:00 To: Gabor Grothendieck [EMAIL PROTECTED] Subject: Re: [R] Lines between panels in lattice Certainly.

Re: [R] Plot availability

2008-09-19 Thread hadley wickham
On Fri, Sep 19, 2008 at 1:00 PM, john crepezzi [EMAIL PROTECTED] wrote: I'm not so sure the original way I stated my issue was clear enough, so I'll attempt to elaborate a little bit. I'd like to make a function that is passed the name of a plot object, and a lines/point specification, and

Re: [R] Novice question about getting data into R

2008-09-19 Thread Ted Byers
Thanks one and all. Actually, I used OpenOffice's spreadsheet to creat the csv file, but I have been using it long enough to know to specify how I wanted it, and sometimes, when that proves annoying, I'll use Perl to finess it the way I want it. It seems my principle error was to assume that it

Re: [R] Plot availability

2008-09-19 Thread Gabor Grothendieck
On Fri, Sep 19, 2008 at 2:00 PM, john crepezzi [EMAIL PROTECTED] wrote: I'm not so sure the original way I stated my issue was clear enough, so I'll attempt to elaborate a little bit. I'd like to make a function that is passed the name of a plot object, and a lines/point specification, and

Re: [R] plot order in multi-panel figure

2008-09-19 Thread Greg Snow
You can reset the figure number with the 'mfg' parameter to 'par', but to do what you want you may be able to do using other tools. Look at the examples for the cnvrt.coords function in the TeachingDemos package (you can do similar things with the grconvertX and grconvertY functions) or look

Re: [R] caret package: arguments passed to the classification or regression routine

2008-09-19 Thread Max Kuhn
Forgot to cc... On Fri, Sep 19, 2008 at 3:22 PM, Max Kuhn [EMAIL PROTECTED] wrote: A new version 3.41 is on https://r-forge.r-project.org/projects/caret/ Until later tonight, you will have to get it via svn checkout svn://svn.r-forge.r-project.org/svnroot/caret and build it yourself.

Re: [R] migrating data from s-plus to R

2008-09-19 Thread Mike Prager
Patrizio Frederic [EMAIL PROTECTED] wrote: Dear all, is there any way to transform a .Data directory created in S-plus 6.1 for windows in a .RData file? You might try the functions dump or dput, which create text representations readable by either S or R. At least, they do so for many common

Re: [R] Plot availability

2008-09-19 Thread john crepezzi
I believe I've found a way to wrap and use dev.cur()[[1]] to have the desired effect. I'm just leaving my terminal now, but I'll post these results on monday so everyone can see how it turned out. If you have any ideas or qualms with this method, please let me know. Thanks Duncan! --john On

[R] lme: NaNs in hierarchical model

2008-09-19 Thread melanie.r
Dear R-users, I'm somehow new to the topic of mixed models. I'd liked to simulate data with a specific correlation-structure and analyze it with lme(). I have a fixed effect (group) and a random effect (block): ygroup block -5.2158969-2 1 -2.4174197-2 1

[R] gbm package on 64 bit machine

2008-09-19 Thread jtf2076
I'm having trouble getting the gbm package to load on my new 64 bit machine running xp professional OS. I installed version 2.72 and loaded the package, although I got an error message: The downloaded packages are in C:\Documents and Settings\Administrator\Local

[R] ggplot2 and lattice

2008-09-19 Thread stephen sefick
I am in the process of learning lattice graphics and have looked at ggplot2 a little. I would like to know if there is a tutorial that shows how to convert lattice code into ggplot code and vise versa. I am finally discovering the power of these two packages and would like suggestions to lessen

[R] two questions on Sweave (help in the output + warnings in the output)

2008-09-19 Thread Domenico Vistocco
Dear All, I am working on some slides using LaTeX/Beamer and R/Sweave. I have the two questions below (sorry if they are stupid or already solved but I didn't find solutions on the web). 1) Using the following code: \begin{frame}[containsverbatim] echo=TRUE= help('dim') @ \end{frame} I have

Re: [R] ggplot2 and lattice

2008-09-19 Thread hadley wickham
On Fri, Sep 19, 2008 at 3:27 PM, stephen sefick [EMAIL PROTECTED] wrote: I am in the process of learning lattice graphics and have looked at ggplot2 a little. I would like to know if there is a tutorial that shows how to convert lattice code into ggplot code and vise versa. I am finally

Re: [R] [R-gui] Tinn-R/Latex output

2008-09-19 Thread Frank E Harrell Jr
Mike Prager wrote: David Carslaw [EMAIL PROTECTED] wrote: I find Tinn-R to be an excellent editor for R, but I have one question I have not been able to answer. I wish to include some R code in Latex. Using the Edit/Copy formatted (to export)/TeX does provide Latex-type output. However,

Re: [R] ggplot2 and lattice

2008-09-19 Thread stephen sefick
yes a parallel coordinates plot- I understand that it is for multivariate data, but I am having a hard time figuring out what it is telling me. Thanks for your help. On Fri, Sep 19, 2008 at 5:02 PM, hadley wickham [EMAIL PROTECTED] wrote: On Fri, Sep 19, 2008 at 3:27 PM, stephen sefick [EMAIL

Re: [R] [R-gui] Tinn-R/Latex output

2008-09-19 Thread hadley wickham
On Fri, Sep 19, 2008 at 4:05 PM, Frank E Harrell Jr [EMAIL PROTECTED] wrote: Mike Prager wrote: David Carslaw [EMAIL PROTECTED] wrote: I find Tinn-R to be an excellent editor for R, but I have one question I have not been able to answer. I wish to include some R code in Latex. Using the

Re: [R] [R-gui] Tinn-R/Latex output

2008-09-19 Thread Frank E Harrell Jr
hadley wickham wrote: On Fri, Sep 19, 2008 at 4:05 PM, Frank E Harrell Jr [EMAIL PROTECTED] wrote: Mike Prager wrote: David Carslaw [EMAIL PROTECTED] wrote: I find Tinn-R to be an excellent editor for R, but I have one question I have not been able to answer. I wish to include some R code

Re: [R] Manual to learn TCL - TK

2008-09-19 Thread Duncan Murdoch
On 19/09/2008 6:25 PM, ermimi wrote: Hello! Anybody knows any manual to learn TCL - TK? I need it to finish my finsish projet. There are a lot of books available, and the docs are included with the Windows distribution of R, in R_HOME/Tcl/doc/tcltk85.chm. Not sure about other platforms.

Re: [R] how to keep up with R?

2008-09-19 Thread Steven McKinney
One thing that R does very well that SAS does not is graphics - graphical portrayal of data is important, and you can keep up with R by supplementing your SAS analyses with R graphics. Steve McKinney -Original Message- From: [EMAIL PROTECTED] on behalf of Barry Rowlingson Sent: Fri

Re: [R] selecting dataframe values that are not nulls

2008-09-19 Thread Adaikalavan Ramasamy
Ramya, you sent four near identical emails with different subject lines. Since the list is run by unpaid volunteers, please avoid wasting people's time (and yours too) with such redundancies. Please read http://www.r-project.org/posting-guide.html and search the mailing lists and

[R] AUC / ROC for presence only.

2008-09-19 Thread milton ruser
Dear all, I have a probability of presence of distribution of a species of interest (varying from 0 to 1 in continuous form) and I have a set of points where I know that species really occurs. But I don´t have points of absence. So, for each true presence I know the estimated presence. I

Re: [R] Design lrm function

2008-09-19 Thread Frank E Harrell Jr
milicic.marko wrote: Hi, Is it possible to get ROC and accuracy ratio/gini straight out of the Design package? Thanks The print method for lrm prints the ROC area (labeled C). lrm does not print the other 2 measures you listed. It computes a generalized R^2 (much more powerful than