[R] [R-pkgs] stringr: version 0.4

2010-08-25 Thread Hadley Wickham
Strings are not glamorous, high-profile components of R, but they do play a big role in many data cleaning and preparations tasks. R provides a solid set of string operations, but because they have grown organically over time, they can be inconsistent and a little hard to learn. Additionally, they

Re: [R] SEM : Warning : Could not compute QR decomposition of Hessian

2010-08-25 Thread John Fox
Dear Anne, I started to diagram your model but stopped when I noticed some problems: (1) Some variables, such as pays_alti, are clearly endogenous, since they have arrows pointing to them, yet are declared as fixed exogenous variables; that clearly doesn't make sense. (2) You've placed

Re: [R] output values from within a function

2010-08-25 Thread Greg Snow
Put the line: cat(z1,'\n') in your function. You may also want to put the flush.console() command right after that. Hope this helps, -- Gregory (Greg) L. Snow Ph.D. Statistical Data Center Intermountain Healthcare greg.s...@imail.org 801.408.8111 -Original Message- From:

Re: [R] Removing inter-bar spaces in barchart

2010-08-25 Thread Greg Snow
Using the barplot function in base graphics you just set space=0, but that function does not have a box.ratio argument which would imply that you are using something else. If you let us know which function (and which package it is in) then it is easier (possible) for us to help you, even

Re: [R] What does this warning message (from optim function) mean?

2010-08-25 Thread Prof Brian Ripley
You mean 'TRUE': 'T' is a variable in R, with initial value TRUE. On Wed, 25 Aug 2010, Sally Luo wrote: Hi R users, I am trying to use the optim function to maximize a likelihood funciton, and I got the following warning messages. Could anyone explain to me what messege 31 means exactly? Is

Re: [R] approxfun-problems (yleft and yright ignored)

2010-08-25 Thread Greg Snow
The plots did not come through, see the posting guide for which attachments are allowed. It will be easier for us to help if you can send reproducible code (we can copy and paste to run, then examine, edit, etc.). Try finding a subset of your data for which the problem still occurs, then send

Re: [R] lattice help required

2010-08-25 Thread Kay Cichini
.. thanks again, richard. and you swiftly saw the next problem comming up - when using par.settings = list(layout.widths = list(axis.panel = c(1, 0))) getting rid of the double tick labeling would be natural - but i'll leave it at that for today. many thanks, kay Richard M. Heiberger wrote:

Re: [R] GLM outputs in condensed versus expanded table

2010-08-25 Thread Peter Dalgaard
On 08/25/2010 05:17 PM, francogrex wrote: Hi I'm having different outputs from GLM when using a condensed table V1V2 V3 Present Absent 0 0 0 3 12 0 0 1 0 0 0 1 0 0 0 0 1 1 1 0 1 0

Re: [R] several odfWeave questions

2010-08-25 Thread Max Kuhn
Ben,  1a. am I right in believing that odfWeave does not respect the 'keep.source' option?  Am I missing something obvious? I believe it does, since this gets passed directly to Sweave.  1b. is there a way to set global options analogous to \SweaveOpts{} directives in Sweave? (I looked at

Re: [R] package MuMI

2010-08-25 Thread Ben Bolker
Marino Taussig De Bodonia, Agnese agnese.marino09 at imperial.ac.uk writes: I am using the package MuMI to run all the possible combinations of variables in my full model, and select my best models. When I enter my variables in the original model I write them like this lm(y~ a +b +c

Re: [R] Change value of a slot of an S4 object within a method.

2010-08-25 Thread Steve Lianoglou
Howdy, On Wed, Aug 25, 2010 at 1:21 PM, Joris Meys jorism...@gmail.com wrote: Hi Steve, thanks for the tip.  I'll definitely take a closer look at your solution for implementation for future use.  But right now I don't have the time to start rewriting my class definitions. Luckily, I found

Re: [R] Draw a perpendicular line?

2010-08-25 Thread Dennis Murphy
Hi: The function below plots the line segment between two points A and B as well as the normal from C to AB, with a dot as the intersection point (D), which is returned with the function call. The aspect ratio is kept at one so that the orthogonality between the two lines is not distorted by the

Re: [R] What does this warning message (from optim function) mean?

2010-08-25 Thread Ravi Varadhan
Hi, You did not give us any information about your likelihood function, f, nor did you provide a reproducible example. So, I cannot tell for sure whether the parameter estimates are reliable. Ravi. -Original Message- From: r-help-boun...@r-project.org

[R] Checking a point behind two polygons

2010-08-25 Thread CZ
Hi, I am working on a problem of checking whether a point is behind two convex hulls. By 'behind', I mean a point which is outside of two convex hulls can't reach either of these two hulls without reaching the other. My idea is to find the segments between the point and the vertices of a

Re: [R] What does this warning message (from optim function) mean?

2010-08-25 Thread Cheng Peng
the deteminant is a nonpositive value. log(det(...)) produce NaNs... -- View this message in context: http://r.789695.n4.nabble.com/What-does-this-warning-message-from-optim-function-mean-tp2338689p2338719.html Sent from the R help mailing list archive at Nabble.com.

Re: [R] how do R calculates the number of intervals between tick-marks

2010-08-25 Thread Antonio Olinto
Hi, With axTicks I did what I wanted. Thank you all for the attention. Below is the code I did to have a plot with three y axes. Suggestions to improve the routine are welcome. All the best, Antonio Olinto Data in a spreadsheet YEARL NFB NFT 200526,352158 1500

Re: [R] Removing inter-bar spaces in barchart

2010-08-25 Thread Dennis Murphy
Hi: On Wed, Aug 25, 2010 at 11:28 AM, Greg Snow greg.s...@imail.org wrote: Using the barplot function in base graphics you just set space=0, but that function does not have a box.ratio argument which would imply that you are using something else. If you let us know which function (and which

Re: [R] Draw a perpendicular line?

2010-08-25 Thread CZ
Thank you all for your suggestions and especially the codes. Now I am able to solve it by finding the intercepts and slopes for both lines first, and then I can find the coordinates of x and y at the intersection. At first I forgot how to use C to get the intercept of line CD... Thanks

[R] accessing the attr(*,label.table) after importing from spss

2010-08-25 Thread moleps
Dear all, I just received a file from a colleague in spss. The read.spss could not finish the file due to an error (Unrecognized record type 7, subtype 18 encountered in system file) so instead I converted the file using stat-transfer. Looking at my data I see that most labels are in the

Re: [R] Removing inter-bar spaces in barchart

2010-08-25 Thread David Winsemius
On Aug 24, 2010, at 10:20 PM, Jonathan Greenberg wrote: Rhelpers: I'm trying to make a barchart of a 2-group dataset (barchart(x~y,data=data,groups=z,horizontal=FALSE)). My problem is that I can't, for the life of me, seem to get rid of the inter-bar space -- box.ratio set to 1 doesn't

Re: [R] package MuMIn

2010-08-25 Thread Ben Bolker
[cc'ing back to r-help: this is good etiquette so that the responses will be seen by others/ archived for future reference.] On 10-08-25 04:35 PM, Marino Taussig De Bodonia, Agnese wrote: Yes, I meant MuMIn the global formula I introduced was: rc4.mod-lm(central$hunting~ central$year +

Re: [R] Secant Method Convergence (Method to replicate Excel XIRR/IRR)

2010-08-25 Thread Ravi Varadhan
Hi Adrian, Your code for the secant method is correct. I just tweaked it a bit w/o the calendar time feature (assuming that the cash flows will be available on an annual basis) as follows: ANXIRR - function (cashFlow, guess, tol=1.e-04){ npv - function (cashFlow, irr) { n -

Re: [R] Secant Method Convergence (Method to replicate Excel XIRR/IRR)

2010-08-25 Thread Ravi Varadhan
Another approach is to use `uniroot' to find the zero of the NPV function: npv - function (cashFlow, irr) { n - length(cashFlow) sum(cashFlow / (1 + irr)^{0: (n-1)}) } uniroot(f=npv, interval=c(0,1), cashFlow=cashFlow) However, there may be situations where there are no

[R] OT: R for iPhone/iPad OS?

2010-08-25 Thread Carl Witthoft
No, seriously: I've had more than one person at work wonder what math toolset could be loaded onto iOS. So, before Matlab, FreeMat, Mathematia, SciLab, Octave, or numpy (:-) ) produces a version for iPad, any chance someone is working on R for iPad?

Re: [R] OT: R for iPhone/iPad OS?

2010-08-25 Thread Greg Snow
I believe that the iPad conditions for producing apps are incompatible with the GPL (gnuGo was removed from the apps store recently for this reason), so don't hold your breath. There may be some possibility for an app on the iPad that would access R on a server somewhere else, but that is

Re: [R] OT: R for iPhone/iPad OS?

2010-08-25 Thread Erik Iverson
Instructions for installing R on jailbroken devices: http://rwiki.sciviews.org/doku.php?id=getting-started:installation:iphone Carl Witthoft wrote: No, seriously: I've had more than one person at work wonder what math toolset could be loaded onto iOS. So, before Matlab, FreeMat,

Re: [R] OT: R for iPhone/iPad OS?

2010-08-25 Thread Steve Lianoglou
Hi, On Wed, Aug 25, 2010 at 5:37 PM, Greg Snow greg.s...@imail.org wrote: I believe that the iPad conditions for producing apps are incompatible with the GPL (gnuGo was removed from the apps store recently for this reason), so don't hold your breath. There may be some possibility for an

Re: [R] Secant Method Convergence (Method to replicate Excel XIRR/IRR)

2010-08-25 Thread Ravi Varadhan
Yes, the secant method (like Newton Raphson) is not guaranteed to converge, unlike the bisection method, but it has a superlinear convergence (not that this matters much!). Brent's method, which is used in `uniroot', is a reliable and fast method, which is why I suggested it in my previous email.

Re: [R] OT: R for iPhone/iPad OS?

2010-08-25 Thread Marc Schwartz
On Aug 25, 2010, at 4:23 PM, Carl Witthoft wrote: No, seriously: I've had more than one person at work wonder what math toolset could be loaded onto iOS. So, before Matlab, FreeMat, Mathematia, SciLab, Octave, or numpy (:-) ) produces a version for iPad, any chance someone is working on

[R] Looking for an image (R 64-bit on Linux 64-bit) on Amazon EC2

2010-08-25 Thread noclue_
I have found an existing image on Amazon EC2 including R. But unfortunately, it is 32-bit R on 32-bit Linux. Does anybody know if there exists an mage (R 64-bit on Linux 64-bit) on Amazon EC2? Or how can I install 64-bit R on my own Linux instance there? Thanks. -- View this message in

Re: [R] OT: R for iPhone/iPad OS?

2010-08-25 Thread Brian Diggs
On 8/25/2010 2:23 PM, Carl Witthoft wrote: No, seriously: I've had more than one person at work wonder what math toolset could be loaded onto iOS. So, before Matlab, FreeMat, Mathematia, SciLab, Octave, or numpy (:-) ) produces a version for iPad, any chance someone is working on R for iPad?

Re: [R] Looking for an image (R 64-bit on Linux 64-bit) on Amazon EC2

2010-08-25 Thread stephen sefick
You have a 64 bit Linux? If so... Dowload the sources make sure you have all of the dependencies unpack tarball cd to de-compressed directory issue ./configure make sudo make install or maybe use your distros packages managment tool. Am I missing something? On Wed, Aug 25, 2010 at 4:39

Re: [R] OT: R for iPhone/iPad OS?

2010-08-25 Thread Carl Witthoft
Hey, Many thanks to all who responded! I'll pass along the links to my pals. Personally I can't imagine doing R on a virtual keyboard in the first place, but to each his own. Carl __ R-help@r-project.org mailing list

Re: [R] how to plot y-axis on the right of x-axis

2010-08-25 Thread elaine kuo
Thank you. The answers provides the right direction and a code was written accordingly. One more request, if the label of axis X wants to be drawn from 5 to 1 (left to right) rather than 1 to 5, is it fine to change axis (4, at = NULL) ? If so, which value should be input ? Thanks again. Elaine

Re: [R] how to plot y-axis on the right of x-axis

2010-08-25 Thread elaine kuo
Yes, I appreciated your answers which well hit my questions. (esp the perfect model parts). About the plot part, one more question. Is it possible to make the two plots (northern and southern richness) sharing the same Y axis (latitude = 0, the equator) ? In other words, southern richness would

Re: [R] Secant Method Convergence (Method to replicate Excel XIRR/IRR)

2010-08-25 Thread Ravi Varadhan
The secant method converges just fine. Your problem might have occurred due to improper conversion of dates to elapsed time. You want to calculate IRR using year as the time unit, not days. Here is the secant function (modified to account for irregular times) and the results for your example:

Re: [R] lattice help required

2010-08-25 Thread Felix Andrews
Deepayan Sarkar has a function combineLimits() in the development version of the latticeExtra package (i.e. the version on r-forge.r-project.org) which will set common scales in each row or column of your layout. It can also remove the internal axes. # Felix On 26 August 2010 04:43, Kay Cichini

Re: [R] Adding column to dataframe

2010-08-25 Thread nzcoops
To update on this. I ran the same command on a grid of computers with 32gb ram, and it completed in 15 seconds, compared to the ~20 minutes on my desktop. Simply a ram issue as suspected by a few on the list here. Thanks -- View this message in context:

Re: [R] Rotate x-axis label on log scale

2010-08-25 Thread Tim Elwell-Sutton
Hi Jim Thanks so much for this. The updated staxlab function works perfectly! Tim -Original Message- From: Jim Lemon [mailto:j...@bitwrit.com.au] Sent: Wednesday, August 25, 2010 6:36 PM To: tesutton Cc: r-help@r-project.org Subject: Re: [R] Rotate x-axis label on log scale On

[R] How to obtain the graph of fitted values against one variable after estimation?

2010-08-25 Thread Le Wang
Hi there, I have a question regarding the plot command after estimation. Specifically, I estimate a model, say regressing y on x and z. And after estimation, I would like to plot the fitted values against x, but I don't need that for z. The following command always gives two graphs, for both

Re: [R] Secant Method Convergence (Method to replicate Excel XIRR/IRR)

2010-08-25 Thread Ravi Varadhan
You should use cfDate[1] as the time origin. You cannot use 08-24-2010 as the time origin, since that will yield negative times. Here is the correct solution. ANXIRR - function (cashFlow, dates, guess, tol=1.e-04){ npv - function (cashFlow, times, irr) { n - length(cashFlow) sum(cashFlow

Re: [R] How to obtain the graph of fitted values against one variable after estimation?

2010-08-25 Thread David Winsemius
On Aug 25, 2010, at 10:46 PM, Le Wang wrote: Hi there, I have a question regarding the plot command after estimation. Specifically, I estimate a model, say regressing y on x and z. And after estimation, I would like to plot the fitted values against x, but I don't need that for z. The

Re: [R] Secant Method Convergence (Method to replicate Excel XIRR/IRR)

2010-08-25 Thread Ravi Varadhan
I have written a general root-finder using the secant method. secant - function(par, fn, tol=1.e-07, itmax = 100, trace=TRUE, ...) { # par = a starting vector with 2 starting values # fn = a function whose first argument is the variable of interest # if (length(par) != 2) stop(You must specify

[R] list of closures

2010-08-25 Thread Stephen T.
Hi, I wanted to create a list of closures. When I use Map(), mapply(), lapply(), etc., to create this list, it appears that the wrong arguments are being passed to the main function. For example: Main function: adder - function(x) function(y) x + y Creating list of closures with Map(): plus -

Re: [R] list of closures

2010-08-25 Thread Gabor Grothendieck
On Thu, Aug 26, 2010 at 12:04 AM, Stephen T. obsessiv...@hotmail.com wrote: Hi, I wanted to create a list of closures. When I use Map(), mapply(), lapply(), etc., to create this list, it appears that the wrong arguments are being passed to the main function. For example: Main function:

Re: [R] list of closures

2010-08-25 Thread Philippe Grosjean
Unless for learning R, you should really consider R.oo or proto packages that may be more convenient for you (but you don't provide enough context to tell). Best, Philippe On 26/08/10 06:28, Gabor Grothendieck wrote: On Thu, Aug 26, 2010 at 12:04 AM, Stephen T.obsessiv...@hotmail.com wrote:

Re: [R] Secant Method Convergence (Method to replicate Excel XIRR/IRR)

2010-08-25 Thread Adrian Ng
Hi Ravi, Using days and dividing it by 365 effectively converts the number to years anyway and allows for the irregular times to be specific to the days. Also, when I replace dates[1] in your line: times - as.numeric(difftime(dates, dates[1], units=days) / 365.24) with 2010-08-24 I think I am

Re: [R] creation package

2010-08-25 Thread anderson nuel
Dear r-help, I took your advice into consideration and i tried to solve my errors. But , when I use the command R CMD check namepackage, I find an error in this file C:/Rpack/namepackage.Rcheck/00check.txt : * using log directory 'C:/Rpack/namepackge.Rcheck' * using R version 2.10.1

Re: [R] Looking for an image (R 64-bit on Linux 64-bit) on Amazon EC2

2010-08-25 Thread noclue_
You have a 64 bit Linux? If so... Dowload the sources Do you mean download Linux kernel source code and then compile it on Amazon EC2? -- View this message in context: http://r.789695.n4.nabble.com/Looking-for-an-image-R-64-bit-on-Linux-64-bit-on-Amazon-EC2-tp2338938p2339072.html Sent

[R] Specify non-default libiconv, readline, libpng, tiff etc during R compilation

2010-08-25 Thread Derrick LIN
Hi guys, I am trying to compile 64-bit R 2.11.1 on Solaris 10. I mainly follow the guide in here https://www.initworks.com/wiki/pages/viewpage.action?pageId=6521038 The guide suggests that install the customised libiconv, readline under the designated R installation folder and become the

<    1   2