Re: [R] removed data is still there!

2010-09-23 Thread Peter Ehlers
On 2010-09-21 5:51, Nikhil Kaza wrote: example(factor) iris1$Species- factor(iris1$Species, drop=T) will get you what you need. Hmm, doesn't work for me. ?factor does not list a 'drop=' argument. -Peter Ehlers Nikhil Kaza Asst. Professor, City and Regional Planning University of North

Re: [R] Plotting multiple animal tracks against Date/Time

2010-09-23 Thread Jonathan Christensen
Include individual as a factor in your dataset, and use ggplot2: library(ggplot2) ggplot(aes(x=Date, y=Distance, color=Individual), data=data) + geom_line() ought to do it. Jonathan On Thu, Sep 23, 2010 at 9:31 AM, Struve, Juliane j.str...@imperial.ac.ukwrote: Sorry for posting this

[R] how to make point character thicker in xyplot?

2010-09-23 Thread array chip
Is there anyway to make plotting point character being thicker in xyplot? I mean not larger which can achieved by cex=2, but thicker. I tried lwd=2, but it didn't work. I know lwd works in regular plot() not only for lines, but also for points. For example plot(1:10, lwd=2) Thanks John

Re: [R] how to make point character thicker in xyplot?

2010-09-23 Thread Greg Snow
There is probably a simpler way, but if you want full customization, look at panel.my.symbols in the TeachingDemos package. -Original Message- From: array chip arrayprof...@yahoo.com Sent: Thursday, September 23, 2010 2:52 PM To: r-help@r-project.org r-help@r-project.org Subject: [R] how

Re: [R] plotting multiple animal tracks against Date/Time

2010-09-23 Thread Gabor Grothendieck
On Thu, Sep 23, 2010 at 5:53 PM, Struve, Juliane j.str...@imperial.ac.ukwrote: Hello, thank you very much for replying. I have tried this, but I get error message Error in .subset(x, j) : invalid subscript type 'list' after z1 - read.zoo(textConnection(Lines1), skip = 1, index = list(1,

Re: [R] how to make point character thicker in xyplot?

2010-09-23 Thread array chip
Thank you Greg. I also got it work by using panel.points (lwd=2) instead of using panel.xyplot() - Original Message From: Greg Snow greg.s...@imail.org To: array chip arrayprof...@yahoo.com; r-help@r-project.org r-help@r-project.org Sent: Thu, September 23, 2010 2:48:06 PM

[R] how to make point character thicker in the legend of xyplot?

2010-09-23 Thread array chip
Now I got point character thicker using panel.points(lwd=2), But I can't make it thicker in the legend of the plot. Here is an example: xyplot(1:10~1:10,groups=rep(1:5,2),col=1:2,pch=0:1,type='p',cex=2,panel=panel.points,lwd=3, key=list(x=0.7,y=0.2,corner=c(0,0),type='p',

Re: [R] how to make point character thicker in the legend of xyplot?

2010-09-23 Thread Peter Ehlers
On 2010-09-23 17:15, array chip wrote: Now I got point character thicker using panel.points(lwd=2), But I can't make it thicker in the legend of the plot. Here is an example: xyplot(1:10~1:10,groups=rep(1:5,2),col=1:2,pch=0:1,type='p',cex=2,panel=panel.points,lwd=3,

Re: [R] how to make point character thicker in the legend of xyplot?

2010-09-23 Thread array chip
Yes, it does what I want. Thank you Peter! Just wondering what else grid.pars controls? not just the symbol in legend, right? John - Original Message From: Peter Ehlers ehl...@ucalgary.ca To: array chip arrayprof...@yahoo.com Cc: r-help@r-project.org r-help@r-project.org Sent: Thu,

Re: [R] how to make point character thicker in the legend of xyplot?

2010-09-23 Thread Peter Ehlers
On 2010-09-23 17:57, array chip wrote: Yes, it does what I want. Thank you Peter! Just wondering what else grid.pars controls? not just the symbol in legend, right? John You can have a look at ?gpar (after loading package 'grid'). For example, your original request could be handled with

Re: [R] non-linear integer optimization?

2010-09-23 Thread Hans W Borchers
darckeen darckeen at hotmail.com writes: This is an example of the type of problem, and how i'm currently using optim() to solve it. I would classify this as an integer programming (IP) problem, it is not even non-linear, as there are no continuous variables. Your approach with optim() will

[R] Solving equations involving gamma functions

2010-09-23 Thread Shant Ch
Hi, I want to find a value of n1. I used the following code but I am getting the error - Error in as.vector(x, mode) : cannot coerce type 'closure' to vector of type 'any' n=10 a_g-(1/(n*(n-1)))*((pi/3)*(n+1)+(2*sqrt(3)*(n-2))-4*n+6) a_s-function(n1) { t1=(n1-1)/2;

Re: [R] empirical df, cdf and crossing points

2010-09-23 Thread Michael Bedward
Hello, Not sure about packages to suggest but some of this can be done with base stats... # Generate densities using common from, to and n args male - rnorm(100, 80, 10) female - rnorm(100, 60, 10) male.d - density(male, from=0, to=150, n=1024) female.d - density(female, from=0, to=150, n=1024)

<    1   2