[R] How to store the parameter estimated by nls( ) to a variable?

2007-08-12 Thread Yuchen Luo
Dear Colleagues. I believe this should be a problem encountered by many: nls( ) is a very useful and efficient function to use if we are just to display the estimated value on screen. What if we need R to store the estimated parameter in a variable? For example: x=rnorm(10, mean=1000, sd=10)

[R] Convert factor to numeric vector of labels

2007-08-12 Thread Falk Lieder
Hi, I have imported a data file to R. Unfortunately R has interpreted some numeric variables as factors. Therefore I want to reconvert these to numeric vectors whose values are the factor levels' labels. I tried as.numeric(factor), but it returns a vector of factor levels (i.e. 1,2,3,...) instead

Re: [R] Convert factor to numeric vector of labels

2007-08-12 Thread Prof Brian Ripley
See the FAQ Q7.10 (and please study the posting guide) On Sun, 12 Aug 2007, Falk Lieder wrote: Hi, I have imported a data file to R. Unfortunately R has interpreted some numeric variables as factors. Therefore I want to reconvert these to numeric vectors whose values are the factor levels'

[R] question on glmmML compared to NLMIXED

2007-08-12 Thread Fluss
Hello! Can anyone help me. I am using the posterior.mode from the result of glmmML. It apears to be different from the BLUe estimate of the RANDOM statement in PROC NLMIXED in SAS. Why is that? Thank you Ronen [[alternative HTML version deleted]]

Re: [R] R Weka and cobweb

2007-08-12 Thread Hans W. Borchers
adschai at optonline.net writes: Hi, I never use cobweb before and I'm quite new to this. I have a couple of questions around the cobweb implementation in R Weka. If you could supply answer or insight, I would really appreciate. 1. From Fisher's paper in 1987, it seems that Cobweb only

Re: [R] Artifacts in pdf() of image() (w/o comments)

2007-08-12 Thread Duncan Murdoch
Michael Kubovy wrote: Dear r-helpers, In my previous message there were comments in the code that may have made cutting and pasting awkward. Here it is w/o them. I have two questions: (1) The following produces a pdf with artifacts. How do I prevent them? What artifacts do you see?

Re: [R] How to store the parameter estimated by nls( ) to a variable?

2007-08-12 Thread Duncan Murdoch
Yuchen Luo wrote: Dear Colleagues. I believe this should be a problem encountered by many: nls( ) is a very useful and efficient function to use if we are just to display the estimated value on screen. What if we need R to store the estimated parameter in a variable? For example:

[R] Write values on y axe

2007-08-12 Thread akki
Hi, I have values on y axe from 0.0001 to 3.086. When I do plot I have writen values: 0.001, 0.050,1.000 ..., but how I can write on graph the minimum value and maximum value, with all decimals (I don't want to use the format 1e-0x)? I am using log scale. For example, if I have the values:

Re: [R] Write values on y axe

2007-08-12 Thread jim holtman
Does this do what you want: x - runif(10) plot(x) # put min/max in red axis(2, at=round(range(x), 4), col.axis='red', las=2) On 8/12/07, akki [EMAIL PROTECTED] wrote: Hi, I have values on y axe from 0.0001 to 3.086. When I do plot I have writen values: 0.001, 0.050,1.000 ..., but how I can

[R] Legend on graph

2007-08-12 Thread akki
Hi, I have a problem when I want to put a legend on the graph. I do: legend(topright, names(o), cex=0.9, col=plot_colors,lty=1:5, bty=n) but the legend is writen into the graph (graphs' top but into the graph), because I have values on this position. How can I write the legend on top the graph

Re: [R] need help to manipulate function and time interval

2007-08-12 Thread Uwe Ligges
Matthew Walker wrote: Uwe Ligges wrote: Henrique Dallazuanna wrote: Hi, Try whit: if(time[j] = 18:00:00 23:59:59) This code is obviously wrong and does not help for the next few lines in the questioner's message, please do not post unsensible stuff. Uwe Ligges

[R] competing risk p value in a figure

2007-08-12 Thread phguardiol
Dear R users Is there a way to include in a figure the 2 sided p value obtained with the crr function in the package cmprsk Such as the one for cuminc: text(3.5,0.05, adj=0, paste(Gray test: p-value = , round(fit1$Tests[1,2],3))) thanks for your help ? Phil [[alternative HTML

Re: [R] Help using gPath

2007-08-12 Thread hadley wickham
Here's a partial extract from a sample session after running your code (NOTE this is using the development version of R; grid.ls() does not exist in R 2.5.1 or earlier): Inspect the grob tree with grid.ls() (similar to Hadley's current.grobTree(), but with different formatting) ... (I'll

Re: [R] Artifacts in pdf() of image() (w/o comments)

2007-08-12 Thread Michael Kubovy
On Aug 12, 2007, at 6:24 AM, Duncan Murdoch wrote: Michael Kubovy wrote: Dear r-helpers, In my previous message there were comments in the code that may have made cutting and pasting awkward. Here it is w/o them. I have two questions: (1) The following produces a pdf with artifacts.

Re: [R] How to store the parameter estimated by nls( ) to a variable?

2007-08-12 Thread Yuchen Luo
Dear Professor Murdoch. Thank you so much Best Wishes Yuchen Luo On 8/12/07, Duncan Murdoch [EMAIL PROTECTED] wrote: Yuchen Luo wrote: Dear Colleagues. I believe this should be a problem encountered by many: nls( ) is a very useful and efficient function to use if we are just to

Re: [R] Legend on graph

2007-08-12 Thread jim holtman
If you are asking to have the values plotted on top of the legend, then you can do the following: plot(x, y, type='n', ...) # create plot, but don't plot legend('topright', ...) lines(x,y) # now plot the data If you want it outside the plot, check the archives for several examples. On 8/12/07,

Re: [R] deldir package - voronoi

2007-08-12 Thread Rolf Turner
On 12/08/2007, at 1:22 PM, zubin wrote: Hello! I am using the deldir package to visualize my data, pretty neat. However, i need to fill the colors using polycol =, fill with colors like a heatmap - more of a gradient fill. The only colors i get are very blocky - how do i assign the

[R] How to control the number format on plot axes ?

2007-08-12 Thread Sébastien
Dear R-users, Basically, everything is in the title of my e-mail. I know that some threads from the archives have already addressed this question but they did not really give a clear solution. Here is a series of short codes that will illustrate the problem: # First a-1:10

Re: [R] Colours in R, was: deldir package - voronoi

2007-08-12 Thread Wolfgang Huber
Dear Rolf and Binabina perhaps this is of use to some: Colour for Presentation Graphics. Ross Ihaka. www.stat.auckland.ac.nz/~ihaka/colour/color.pdf Choosing Color Palettes for Statistical Graphics Achim Zeileis and Kurt Hornik. eeyore.ucdavis.edu/stat250/epub-wu-01_abd.pdf Best wishes

Re: [R] Colours in R, was: deldir package - voronoi

2007-08-12 Thread Rolf Turner
On 13/08/2007, at 11:27 AM, Wolfgang Huber wrote: Dear Rolf and Binabina perhaps this is of use to some: Colour for Presentation Graphics. Ross Ihaka. www.stat.auckland.ac.nz/~ihaka/colour/color.pdf Choosing Color Palettes for Statistical Graphics Achim Zeileis and Kurt Hornik.

Re: [R] How to control the number format on plot axes ?

2007-08-12 Thread jim holtman
Here is a way that you can put the formatting that you want; you were not clear on exactly what you were after. You can setup the 'labels' argument for whatever you want. a-1:10 myTicks-c(0.1,1,2,5,10) # set ylim to range of myTicks that you want plot(x=a,y=a,log=y,type=p,yaxt=n,

[R] Data structure in R

2007-08-12 Thread adschai
Hi, I have a question around how to build data structure in R. I have to implement a tree structure of data. I'm wondering if R already has something like a tree where I can extend from or I need to start from scratch. If so, what would be the most effective way to represent parent and child

[R] How to write to a table column by column?

2007-08-12 Thread Yuchen Luo
Dear friends. Every loop of my program will result in a list that is very long, with a structure similar to the one below: Lst - list(name=Fred, wife=Mary, daily.incomes=c(1:850)) Please notice the large size of daily.incomes. I need to store all such lists in a csv file so that I can easily

Re: [R] How to control the number format on plot axes ?

2007-08-12 Thread Sébastien
Thanks Jim, This sprintf function was exactly what I was looking for. jim holtman a écrit : Here is a way that you can put the formatting that you want; you were not clear on exactly what you were after. You can setup the 'labels' argument for whatever you want. a-1:10

Re: [R] Data structure in R

2007-08-12 Thread J. Burgos
Hi Adschai, R has a class for tree-like structures. See the 'tree' package. Julian On Mon, 13 Aug 2007, [EMAIL PROTECTED] wrote: Hi, I have a question around how to build data structure in R. I have to implement a tree structure of data. I'm wondering if R already has something like a

Re: [R] Artifacts in pdf() of image() (w/o comments)

2007-08-12 Thread Duncan Murdoch
Michael Kubovy wrote: On Aug 12, 2007, at 6:24 AM, Duncan Murdoch wrote: Michael Kubovy wrote: Dear r-helpers, In my previous message there were comments in the code that may have made cutting and pasting awkward. Here it is w/o them. I have two questions: (1) The following

[R] To succeed in business!2129

2007-08-12 Thread Calvin Mckenzie
2050 Mynes C Fin is ONE of the leading providers a consulting at the world. Our success depends both on high quality of services and on quality managed reliable business processes. This is the reason why quality of our base concern. However, the only way to reach supreme quality at our

[R] A clean way to initialize class slot of type numeric vector

2007-08-12 Thread adschai
Hi, I have a class definition like this: setClass(foo, representation(members=numeric), prototype(members=c())) I intend my class to have members, a slot whose value should be a vector of integer. When I initialize this class, I don't have any member yet. So my member is blank. But if I

Re: [R] A clean way to initialize class slot of type numeric vector

2007-08-12 Thread Seth Falcon
[EMAIL PROTECTED] writes: Hi, I have a class definition like this: setClass(foo, representation(members=numeric), prototype(members=c())) I intend my class to have members, a slot whose value should be a vector of integer. When I initialize this class, I don't have any member yet. So my