[R] Gamma

2009-05-16 Thread Kon Knafelman
Hey, I am having trouble graphing the following function √2Γ(n/2)/[√n − 1Γ((n − 1)/2 for the values of n between 2 and 50. i know that Γ(n) = (n-1)!, which in R is factorial(n-1) When i type that into R, using y - function(n). and then plot(y,2,50), it doesnt give me anything

[R] How to turn the xyplot graphics upside-down?

2009-05-16 Thread lau kue
Hi everyone, I'm learning to use R. The below script works well, but i want to turn the y-axis and the graphics upside-down, how should i do? Appreciate if any tips! library(lattice)

[R] Gamma function

2009-05-16 Thread Kon Knafelman
Hi Guy, I am having trouble graphing the following function √2Γ(n/2)/[√n − 1Γ((n − 1)/2 for the values of n between 2 and 50. i know that Γ(n) = (n-1)!, which in R is factorial(n-1) When i type that into R, using y - function(n). and then plot(y,2,50), it doesnt give me

[R] Fwd: Cannot allocate a new database connection error

2009-05-16 Thread Moumita Das
-- Forwarded message -- From: Moumita Das das.moumita.onl...@gmail.com Date: Sat, May 16, 2009 at 2:26 PM Subject: Cannot allocate a new database connection error To: r-help-requ...@r-project.org Hi friends, why do i keep getting this error?The program runs, twice and every third

Re: [R] Gamma

2009-05-16 Thread Johannes Huesing
Kon Knafelman konk2...@hotmail.com [Sat, May 16, 2009 at 08:49:44AM CEST]: Hey, I am having trouble graphing the following function ???2??(n/2)/[???n ??? 1??((n ??? 1)/2 for the values of n between 2 and 50. i know that ??(n) = (n-1)!, which in R is factorial(n-1) When i type that

[R] Installing R on Ubuntu 8.10

2009-05-16 Thread Maria I. Tchalakova
Hi, I am trying to install R on Ubuntu 8.10. I have added this line: deb http://cran.at.r-project.org/bin/linux/ubuntu intrepid/ in my /etc/apt/sources.list file. However when I Reload, I get the following error message: W: GPG error: http://cran.at.r-project.org intrepid/ Release: The

Re: [R] Installing R on Ubuntu 8.10

2009-05-16 Thread Nordlund, Dan (DSHS/RDA)
-Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Maria I. Tchalakova Sent: Saturday, May 16, 2009 2:37 AM To: R-help Forum Subject: [R] Installing R on Ubuntu 8.10 Hi, I am trying to install R on Ubuntu 8.10. I have added

Re: [R] drawing arrows

2009-05-16 Thread Jim Lemon
christophe dutang wrote: Hi, I would like to draw arrows in a classic 2D plot. Which package should I use? is there R base functions that do job? On google, I could not find any useful discussion about this topic, except a link to the function 'grid.arrows' of the grid package. My problem is

Re: [R] Installing R on Ubuntu 8.10

2009-05-16 Thread Maria I. Tchalakova
Hi Dan, thanks for the reply. Yes, I have read everything there, but it didn't help. I've written to the Ubuntu forums as well, from where I was told to type first: gpg --keyserver keyserver.ubuntu.com --recv D67FC6EAE2A11821 gpg --export --armor D67FC6EAE2A11821 | sudo apt-key add - before

[R] newbie: closing unused connection + readline

2009-05-16 Thread Aval Sarri
Hello; I am new to R and trying to read a line from socket connection at a time but at the end of script I am getting closing unused connection warning. I am not able to understand how to solve this. I want to read a line from socket and then use read.table/scan on that line but it looks like I am

Re: [R] newbie: closing unused connection + readline

2009-05-16 Thread Gabor Grothendieck
Try explicitly closing it: close(reader.socket) or closeAllConnections() On Sat, May 16, 2009 at 8:34 AM, Aval Sarri aval.sa...@gmail.com wrote: Hello; I am new to R and trying to read a line from socket connection at a time but at the end of script I am getting closing unused connection

[R] barplot in a loop not printing out to directory

2009-05-16 Thread Dimitri Liakhovitski
Hello! I am experiencing a problem with section 3 of the code below. I want to generate barplots (based on data generated in Sections 1 and 2) in a loop - for each variable in data - and save them as .emf files in my current directory. But it's not working - it's printing values to be plotted to

Re: [R] fitdistr for t distribution

2009-05-16 Thread Martin Maechler
l == lagreene lagreene...@gmail.com on Fri, 15 May 2009 04:22:59 -0700 (PDT) writes: l Thanks Jorge, l but I still don't understand where they come from. when I use: l fitdistr(mydata, t, df = 9) and get values for m and s, and the variance l of my data should be the

Re: [R] ggplot2: annotating plot with mathematical formulae

2009-05-16 Thread Paul Emberson
Hi Stephen, The problem is that the label on the graph doesn't get rendered with a superscript. I want the label on the graph to be rendered the same way as the label you have put on the axis. I am plotting a piecewise function and I wanted to label each section of it. Paul stephen sefick

Re: [R] barplot in a loop not printing out to directory

2009-05-16 Thread Jim Lemon
Dimitri Liakhovitski wrote: Hello! I am experiencing a problem with section 3 of the code below. I want to generate barplots (based on data generated in Sections 1 and 2) in a loop - for each variable in data - and save them as .emf files in my current directory. But it's not working - it's

Re: [R] newbie: closing unused connection + readline

2009-05-16 Thread Aval Sarri
On Sat, May 16, 2009 at 6:12 PM, Gabor Grothendieck ggrothendi...@gmail.com wrote: Try explicitly closing it: close(reader.socket) or closeAllConnections() Thank you Sir but is there a way that it does not get open each time? I mean following line creates a new text connection each time so

Re: [R] anova(cph(..) output

2009-05-16 Thread Frank E Harrell Jr
pompon wrote: Hello, I am a beginner in R and statistics, so my question may be trivial. Sorry in advance. I performed a Cox proportion hazard regression with 2 categorical variables with cph{design}. Then an anova on the results. the output is anova(cph(surv(survival, censor) ~ plant +

Re: [R] Fwd: Cannot allocate a new database connection error

2009-05-16 Thread Uwe Ligges
Open a connection as in con - myDbConnect(dbName,hostName) before the loop starts and afterwards just work on con. Uwe Ligges Moumita Das wrote: -- Forwarded message -- From: Moumita Das das.moumita.onl...@gmail.com Date: Sat, May 16, 2009 at 2:26 PM Subject: Cannot

Re: [R] Matrix package,solve() errors and crashes Please help

2009-05-16 Thread Martin Maechler
SS == Surendar Swaminathan surendar.swaminat...@gmail.com on Fri, 15 May 2009 15:55:23 -0700 writes: Hello All, SS Please help me with this problem.I have been having this problem for over a SS month now and I could not find any information.I later realised that error

Re: [R] Memory size problems in R

2009-05-16 Thread Uwe Ligges
Fernando José Nascimento Sebastião wrote: Dear all, I have some problems how to increase memory in R. I tried to start R by calling e.g. :(in the Properties, Shortcut tab, Target field) Rgui.exe --max-mem-size=2Gb (Is this correct?) But the error message still continue:

Re: [R] Plotting question re. cuminc

2009-05-16 Thread Uwe Ligges
K F Pearce wrote: Hello everyone, (This is my second question posted today on the R list). and you have still not read the posting guide? It asks you to provide commented, minimal, self-contained, reproducible code which also means that you need to tell which packages you are using. I

Re: [R] data summary and some automated t.tests.

2009-05-16 Thread Uwe Ligges
stephen sefick wrote: I would like to preform a t.test to each of the measured variables (sand.silt etc.) I am a big fan of applying t.test()s, but in this case: Are you really sure? The integers and particularly boxplot(x) do not indicate very well that the variables are somehow close to

Re: [R] (no subject)

2009-05-16 Thread Linlin Yan
On Sat, May 16, 2009 at 12:05 PM, Debbie Zhang debbie0...@hotmail.com wrote: Dear R users, Does anyone know how to write a function involving derivative? i.e. I want to implementing Newton's method in R, so my function is something like x- x-y/y' I am not sure how to write y' in my

Re: [R] Fw: Help with loops(corrected question)

2009-05-16 Thread Uwe Ligges
Amit Patel wrote: --- On Fri, 15/5/09, Amit Patel amitrh...@yahoo.co.uk wrote: From: Amit Patel amitrh...@yahoo.co.uk Subject: Help with loops To: r-help@r-project.org Date: Friday, 15 May, 2009, 12:17 PM Hi I am trying to create a loop which averages replicates in my data. The original

[R] maxLik pakage

2009-05-16 Thread amene kheradmandi
Hi all; I recently have been used 'maxLik' function for maximizing G2StNV178 function with gradient function gradlik; for receiving this goal, I write the following program; but I have been seen an error  in calling gradient  function; The maxLik function can't enter gradlik function

[R] MANOVA

2009-05-16 Thread Davide Floriello
Dear Sir, I am Davide Floriello and I am a student in mathematics. I have got a problem with the MANOVA commands. I write here what I have done: PV - read.table('PV.txt') PV - data.frame(PV) g - 12 p - 2 REL - factor(PV$REL) HPV - factor(PV$HPV) fit - manova(as.matrix(PV) ~ REL + HPV +

[R] vector fragment

2009-05-16 Thread jonas garcia
Dear R users: I have got a simple question that has been bothering me for a while. Given a certain character vector, I would like to get in a separate vector a fragment of text, in this case the 3rd and 4th letters of each element. So, if: v- c(“stratosphere”, “mesosphere”,

Re: [R] vector fragment

2009-05-16 Thread Jorge Ivan Velez
Dear Jonas, Try this: v- c('stratosphere', 'mesosphere', 'troposphere') substr(v,3,4) [1] ra so op See ?substr for more information. HTH, Jorge On Sat, May 16, 2009 at 9:52 AM, jonas garcia garcia.jona...@googlemail.com wrote: Dear R users: I have got a simple question that has been

Re: [R] maxLik pakage

2009-05-16 Thread spencerg
1. Have you worked through the examples in the maxLik help page? Your example is sufficiently complicated that I hesitate to try it myself, especially since I see characters in your email that are not simple ASCII. If you can get the examples in the maxLik help page to work, identify the

[R] bagged importance estimates in earth problem

2009-05-16 Thread Joseph Retzer
I was trying to produced bagged importance estimates of attributes in earth using the caret package with the following commands:   fit2 - bagEarth(loyalty ~ ., data=model1, B = 10)   bagImpGCV - varImp(fit2,value=gcv) My bootstrap estimates are produced however the second command varImp

Re: [R] barplot in a loop not printing out to directory

2009-05-16 Thread Dimitri Liakhovitski
Thanks a lot. The problem was that I was naming my files illegally. I was starting them with a number. As soon as I changed that, everything worked. Dimitri On Sat, May 16, 2009 at 9:12 AM, Jim Lemon j...@bitwrit.com.au wrote: Dimitri Liakhovitski wrote: Hello! I am experiencing a problem

Re: [R] ggplot2: annotating plot with mathematical formulae

2009-05-16 Thread hadley wickham
Hi Paul, Unfortunately that's not something that's currently possible with ggplot2, but I am thinking about how to make it possible. Hadley On Sat, May 16, 2009 at 7:48 AM, Paul Emberson em...@calidasoft.co.uk wrote: Hi Stephen, The problem is that the label on the graph doesn't get rendered

Re: [R] Checking a (new) package - examples require other package functions

2009-05-16 Thread Rebecca Sela
After adding the right functions to the NAMESPACE file, I am now getting a different error: Error: package/namespace load failed for 'REEMtree' Call sequence: 2: stop(gettextf(package/namespace load failed for '%s', libraryPkgName(package)), call. = FALSE, domain = NA) 1:

[R] MEXVAL

2009-05-16 Thread Mihai Nica
Greetings: I would like to kindly ask help with obtaining mexval statistics (marginal explanatory value - percentage increase in SEE if the variable were left out of the regression model) for a logit (glm) model with several continuous independent variables. I believe I can do it manually for

Re: [R] ggplot2: annotating plot with mathematical formulae

2009-05-16 Thread baptiste auguie
If you're desperate for a workaround, you might want to try this example using pgfSweave, http://ggplot2.wik.is/Mathematical_annotations On a similar vein, you could try psfrag replacements with a postscript device (there is some code for this on the list archives). Feel free to comment /

Re: [R] fitdistr for t distribution

2009-05-16 Thread Paul Johnson
On Fri, May 15, 2009 at 6:22 AM, lagreene lagreene...@gmail.com wrote: Thanks Jorge, but I still don't understand where they come from.  when I use: fitdistr(mydata, t, df = 9) and get values for m and s, and the variance of my data should be the df/s? I jsut want to be able to confirm how

Re: [R] MANOVA

2009-05-16 Thread John Fox
Dear Davide, -Original Message- From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On Behalf Of Davide Floriello Sent: May-16-09 8:51 AM To: r-help@r-project.org Subject: [R] MANOVA Dear Sir, I am Davide Floriello and I am a student in mathematics. I have

[R] converting numeric to integer

2009-05-16 Thread Thomas Mang
Hello, Suppose I have x, which is a variable of class numeric. The calculations performed to yield x imply that mathematically it should be an integer , but due to round-off errors, it might not be (and so in either direction). The error is however small, so round(x) will yield the

[R] Question about barplot: gridlines value labels

2009-05-16 Thread Dimitri Liakhovitski
Hello! I promise I looked into help files before asking. Still cannot figure it out. I think it's because I am totally confused what packages use lettice, which use trellis, etc. Sections 1 and 2 below produce the data and the data to plot. My question is about barplot in Section 3. I am trying

[R] How to draw a R plot in Java frame ?

2009-05-16 Thread Obi Wan
Hi guys, I am new to R and would like to display a R plot inside my Java app. Is there a possibility to do that? If so, what library do I need? May the force be with you [[alternative HTML version deleted]] __ R-help@r-project.org

Re: [R] data summary and some automated t.tests.

2009-05-16 Thread stephen sefick
Up and down are the treatments. These are replicates within date for percent cover of habiat. This is habitat data for a stream restoration - up is the unrestored and dn is the restored. I have looked at the density plots and they do not look gaussian - you are absolutely right. Even log(n+1)

Re: [R] newbie: closing unused connection + readline

2009-05-16 Thread Stavros Macrakis
On Sat, May 16, 2009 at 8:34 AM, Aval Sarri aval.sa...@gmail.com wrote: # Create a socket from which to read lines - one at a time (record) reader.socket -   socketConnection( host = 'localhost', 5000,                                     server = TRUE, blocking = TRUE,                          

Re: [R] newbie: closing unused connection + readline

2009-05-16 Thread Stavros Macrakis
On Sat, May 16, 2009 at 9:11 AM, Aval Sarri aval.sa...@gmail.com wrote: ...I tried something line this also: mydataframe - read.table (socket, sep=,); but does not work says no input lines. this also. mydataframe - read.table (readLine(socket), sep=,); Sorry, I didn't see this before my

Re: [R] Gamma function

2009-05-16 Thread Stavros Macrakis
What exactly is the R code you wrote for your function f? Without that, it will be hard to help you. -s On Sat, May 16, 2009 at 2:48 AM, Kon Knafelman konk2...@hotmail.com wrote: Hi Guy, I am having trouble graphing the following function √2Γ(n/2)/[√n - 1Γ((n - 1)/2 for the

[R] Lattice scales question: using at when log = TRUE

2009-05-16 Thread Afshartous, David
All, I have a simple lattice plot where I have set log = TRUE for the y scale. When I attempt to change the tick locations via the at argument within scales, the supplied numeric vector is not followed. Any suggestions much appreciated for the example below: y = c(10^1.5, 10^2, 10^3, 10^2) t =

[R] I can't update

2009-05-16 Thread Gilson Sanchez
To whom it concenrt I am trying to update the R program 2.9 This is the message that the program give is: update.packages(ask='graphics') Warning in install.packages(update[instlib == l, Package], l, contriburl = contriburl, : 'lib = C:/PROGRA~1/R/R-29~1.0/library' is not writable Erro em

[R] How to save R clean sessions in BATCH mode?

2009-05-16 Thread mcnda839
Thanks a lot for all of you that have reply me about opening and ending R workspaces in BATCH mode. However replies were a king general and I’m afraid I could not take the entire message from them. Therefore I chose to expose here a representative fraction of my work. I have 50 Rdata files

Re: [R] I can't update

2009-05-16 Thread Gilson Sanchez
To whom it concenrt I am trying to update the R program 2.9 This is the message that the program give is: update.packages(ask='graphics') Warning in install.packages(update[instlib == l, Package], l, contriburl = contriburl, : 'lib = C:/PROGRA~1/R/R-29~1.0/library' is not writable

Re: [R] How to save R clean sessions in BATCH mode?

2009-05-16 Thread Sarah Goslee
Would rm(list=ls()) not be an easier solution, if I understand your question correctly? Sarah On Sat, May 16, 2009 at 11:21 AM, mcnda...@mncn.csic.es wrote: Thanks a lot for all of you that have reply me about opening and ending R workspaces in BATCH mode. However replies were a king general

Re: [R] How to draw a R plot in Java frame ?

2009-05-16 Thread Liviu Andronic
On Sat, May 16, 2009 at 9:18 PM, Obi Wan wan_...@ymail.com wrote: I am new to R and would like to display a R plot inside my Java app. Is there a possibility to do that? If so, what library do I need? Perhaps library(JavaGD) is an option. It is currently used in JGR. Liviu

Re: [R] Lattice scales question: using at when log = TRUE

2009-05-16 Thread Gabor Grothendieck
Your parentheses are wrong. It should be y = list(log = TRUE, at = ...) On Sat, May 16, 2009 at 5:11 PM, Afshartous, David dafshart...@med.miami.edu wrote: All, I have a simple lattice plot where I have set log = TRUE for the y scale. When I attempt to change the tick locations via the

Re: [R] bagged importance estimates in earth problem

2009-05-16 Thread Gavin Simpson
On Sat, 2009-05-16 at 08:57 -0700, Joseph Retzer wrote: I was trying to produced bagged importance estimates of attributes in earth using the caret package with the following commands: fit2 - bagEarth(loyalty ~ ., data=model1, B = 10) bagImpGCV - varImp(fit2,value=gcv) My bootstrap

[R] Updating to 2.9 - ' no slot call in the result'

2009-05-16 Thread Raghu Naik
I updated R (to R-patched) and Zelig and the problem went away. Raghu [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

Re: [R] How to save R clean sessions in BATCH mode?

2009-05-16 Thread Douglas Bates
On Sat, May 16, 2009 at 10:21 AM, mcnda...@mncn.csic.es wrote: Thanks a lot for all of you that have reply me about opening and ending R workspaces in BATCH mode. However replies were a king general and I’m afraid I could not take the entire message from them. Therefore I chose to expose here

Re: [R] fitdistr for t distribution

2009-05-16 Thread spencerg
In addition to seeing the code by typing the name of the function (and copying it from there into a file), you can also enter debug(fitdistr), for example. Then the next time you use fitdistr, either directly or indirectly, it puts you in the environment of that function, and you can

Re: [R] Lattice scales question: using at when log = TRUE

2009-05-16 Thread Afshartous, David
Thanks, but even with typo corrected as below the supplied marks are not followed: xyplot(y ~ t, type = b, scales = list(y = list(log = TRUE, at = c(10^1,10^2, 10^2.5, 10^3 From: Gabor Grothendieck [ggrothendi...@gmail.com] Sent: Saturday, May

Re: [R] Lattice scales question: using at when log = TRUE

2009-05-16 Thread Deepayan Sarkar
On 5/16/09, Afshartous, David dafshart...@med.miami.edu wrote: Thanks, but even with typo corrected as below the supplied marks are not followed: xyplot(y ~ t, type = b, scales = list(y = list(log = TRUE, at = c(10^1,10^2, 10^2.5, 10^3 How so? Looks OK to me. -Deepayan

Re: [R] Lattice scales question: using at when log = TRUE

2009-05-16 Thread Gabor Grothendieck
It works for me on both of these (Windows Vista): R.version.string [1] R version 2.8.1 Patched (2008-12-26 r47350) packageDescription(lattice)$Version [1] 0.17-22 R.version.string [1] R version 2.9.0 Patched (2009-05-03 r48460) packageDescription(lattice)$Version [1] 0.17-22 On Sat, May

[R] Newton's method for finding roots

2009-05-16 Thread Kon Knafelman
Hey guys, i have a relatively simple problem. I need to use netwon's method to find the root of a polynomial, lets say x^3-2x-1 i start off with p - function(x) x^3-2*x-1 My method, which im sure is very amateur, is to type another function, which is the derivative of p, and after picking

Re: [R] Using column length in plot gives error

2009-05-16 Thread Stephanie Kovalchik
Is spectra.wavelength a factor? If so, plot will treat it as categorical and not draw a line. Try the following modification. plot(as.numeric(spectra.wavelength), cormat, type = l, ylim=c(-1,1), xlab=Wavelength (nm), ylab=Correlation) Quoting MikSmith m...@hsm.org.uk: Hi I'm trying to

Re: [R] converting numeric to integer

2009-05-16 Thread Linlin Yan
How about ceiling(x), which return the smallest integer not less than x? On Sun, May 17, 2009 at 2:49 AM, Thomas Mang thomas.m...@fiwi.at wrote: Hello, Suppose I have x, which is a variable of class numeric. The calculations performed to yield x imply that mathematically it should be an

Re: [R] Question about barplot: gridlines value labels

2009-05-16 Thread Duncan Mackay
Dimitri You mentioned lattice so I presumed that you tried a barchart in lattice if you found the class of table.a.percents will it comply with the lattice requirements of the data type ? class(table.a.percents) [1] matrix Using the lattice package and converting to a data.frame x -

Re: [R] ggplot2: annotating plot with mathematical formulae

2009-05-16 Thread Paul Emberson
Hi Baptiste, I think on this occasion I'll edit manually in inkscape with the textext plugin which is a somewhat clumsy but simple option. Thanks for the link to the wiki which I wasn't aware of until now. Paul baptiste auguie wrote: If you're desperate for a workaround, you might want to try

[R] projecting onto 2d *integer* grid?

2009-05-16 Thread haettulegur
I have a bunch of data points in R^2 (Euclidean plane). But I want to project these onto a 2d integer grid -- that is, there's at most one data point for each integer (x, y) coordinate, and points near each other in the integer grid should also be near each other in the original R^2 space. Is

Re: [R] How to turn the xyplot graphics upside-down?

2009-05-16 Thread Ben Bolker
lau kue wrote: Hi everyone, I'm learning to use R. The below script works well, but i want to turn the y-axis and the graphics upside-down, how should i do? Appreciate if any tips! Your example isn't reproducible (i.e. we don't have the data), so can't be sure, but try giving an

Re: [R] Row number of minimum value?

2009-05-16 Thread Ben Bolker
MikSmith wrote: This *must* be an insanely easy thing to work out, but I'm not too familiar with R syntax. So how do I work out the row number (if I pass a column) of the minimum value?? I can get the value itself from min(), but where can I get the row?? ?which ?which.min Ben Bolker

[R] wireframe 3-D problems

2009-05-16 Thread Xu Jun
Dear All R Users, I am trying to graph a 3-D graph of simulated data for logit models using the powerful wireframe command, but I got stuck. Here are the codes: x - seq(-4, 4, by=0.01) y - seq(-4, 4, by=0.01) p - 1/(1+exp(-0.12*x + 0.35*y)) mydata - cbind(x, y, p) require(lattice)