[R] Starting R within an VBA makro

2007-06-11 Thread Maja Schröter
Hello everybody, at work I want to start R within an VBA program. First I calculate something in Excel 2003. After that I want to start within my VBA makro R wich should start an R file, say i.g, superplot.R which plots me the data well. So is it possible? Maybe something like that: sub

Re: [R] Tools For Preparing Data For Analysis

2007-06-11 Thread Chris Evans
(Ted Harding) sent the following at 10/06/2007 09:28: ... much snipped ... (As is implicit in many comments in Robert's blog, and indeed also from many postings to this list over time and undoubtedly well known to many of us in practice, a lot of the problems with data files arise at the

[R] GMM estimation

2007-06-11 Thread Sebastian Kruk
Dear everyone: I have to finish my thesis to graduate as Bs. in Economics. I choose to estimate a New Keynesian Phillips Curve (NKPC) for Uruguay using Generalized Moment Method (GMM). I do not know programming or R but I would like to use it. Should I use gee, geepack or gam? Thanks in

[R] generalized moment method

2007-06-11 Thread Sebastian Kruk
Dear everyone: I have to finish my thesis to graduate as Bs. in Economics. I choose to estimate a New Keynesian Phillips Curve (NKPC) for Uruguay using Generalized Moment Method (GMM). I do not know programming or R but I would like to use it. Should I use gee, geepack or gam? Thanks in

[R] Textpad help

2007-06-11 Thread Mary Royerr
I have installed textpad and tried running R code. But it gives me the following error message. The filename, directory name, or volume label syntax is incorrect. Tool completed with exit code 1 Can you provide any help? I am not a technical person. So the help in detail will be appreciated.

Re: [R] icc from GLMM?

2007-06-11 Thread Gregor Gorjanc
Shinichi Nakagawa S.Nakagawa at sheffield.ac.uk writes: ... I am a little confused which one to trust and use. Or there are no easy form to do this? I am guessing formula would change depending on what distribution you use and what link function as well? I want to calculate icc from GLMM with

Re: [R] Starting R within an VBA makro

2007-06-11 Thread Roger Bivand
On Mon, 11 Jun 2007, Maja Schröter wrote: Hello everybody, at work I want to start R within an VBA program. First I calculate something in Excel 2003. After that I want to start within my VBA makro R wich should start an R file, say i.g, superplot.R which plots me the data well. So is it

[R] epitools and R 2.5

2007-06-11 Thread Pietro Bulian
At work after updating to R 2.5 I get an error using epitab from package epitools, when at home (R 2.4) I get no error. Could someone help me? Thanks Pietro Bulian Servizio di Onco-Ematologia Clinico-Sperimentale I.R.C.C.S. Centro di Riferimento Oncologico Via Franco Gallini 2 33081 Aviano

Re: [R] epitools and R 2.5

2007-06-11 Thread Prof Brian Ripley
On Mon, 11 Jun 2007, Pietro Bulian wrote: At work after updating to R 2.5 I get an error using epitab from package epitools, when at home (R 2.4) I get no error. Could someone help me? The maintainer: this is a long-standing bug in the package. But you have enough information from the error

[R] Looking for R-code for non-negative matrix factorization in the presence of Gaussian or Poisson noise

2007-06-11 Thread christian.ritter
Hi all, Has any of you implemented code for non-negative matrix factorization to solve Y=T P' +E; dim(Y)=n,p ; dim(T)=n,nc; dim (P)=(p,nc); dim(E)=n,p where T and P must be non-negative and E either Gaussian or Poisson noise. I'm looking for two variants: 1. Easy (I think), T is known (that

Re: [R] Lines in dotchart dotplot ?

2007-06-11 Thread John Kane
--- [EMAIL PROTECTED] wrote: On 6/9/07, John Kane [EMAIL PROTECTED] wrote: Is it possible to use dotchart or dotplot and set the lines in such a way that they only extend from the left y-axis to the data point? Yes (sort of) in dotplot at least. E.g., dotplot(VADeaths, groups =

[R] system() and R BATCH

2007-06-11 Thread Meinhard Ploner
If I start from within R a new R batch job by using something like system(R CMD BATCH --no-save --quiet Rin.txt Rout.txt, intern=FALSE, ignore.stderr=TRUE, wait=FALSE, input=NULL) the job runs fine and smooth. However, when, for any reason, I put twice ctrl+C in the calling R,

Re: [R] Textpad help

2007-06-11 Thread John Kane
Hi Mary, You really have not given us much to go on. An example of the code that you were trying to run would be a great help (as it says down below : PLEASE do read the posting guide. The point about code is important. However the first thing to check is your file path. This is my first guess

Re: [R] generalized moment method

2007-06-11 Thread Markus Jäntti
Sebastian Kruk wrote: Dear everyone: I have to finish my thesis to graduate as Bs. in Economics. I choose to estimate a New Keynesian Phillips Curve (NKPC) for Uruguay using Generalized Moment Method (GMM). I do not know programming or R but I would like to use it. Should I use gee,

[R] problem with xlsReadWrite package

2007-06-11 Thread hassen62
Hi friends, I have installed R 2.4.0 in my pc. I have a file xls entitled dali following this directory:c://programfiles//R 2.4.0. Recently I have installed xlsreadwrite 1.3.2. but , when I wrote the following lines: library(xlsReadWrite) read.xls( file, colNames = TRUE, sheet = 1, type =

[R] problem with xlsReadWrite package

2007-06-11 Thread hassen62
Hi friends, I have installed R 2.4.0 in my pc. I have a file xls entitled dali following this directory:c://programfiles//R 2.4.0. Recently I have installed xlsreadwrite 1.3.2. but , when I wrote the following lines: library(xlsReadWrite) read.xls( file, colNames = TRUE, sheet = 1, type =

[R] Error using mgcv package

2007-06-11 Thread spime
Hi all, I need some solution in the following problem. The following error appears when i use mgcv package for implementing GAM. But the same formula works fine in gam package. model.gam - gam(formula = RES ~ CAT01+s(NUM01,5)+CAT02+CAT03+s(NUM02,5)+CAT04+ +

Re: [R] {spam?} Nonlinear Regression

2007-06-11 Thread Martin Henry H. Stevens
Hi tronter, PLEASE do read the posting guide http://www.R-project.org/posting- guide.html and provide commented, minimal, self-contained, reproducible code. Hank On Jun 7, 2007, at 5:50 PM, tronter wrote: Hello I followed the example in page 59, chapter 11 of the 'Introduction to R'

[R] Rounding?

2007-06-11 Thread BXC (Bendix Carstensen)
I was a bit puzzed by: formatC(6.65,format=f,digits=1) [1] 6.6 So I experimented and found: formatC(6.6501,format=f,digits=1) [1] 6.6 formatC(6.651,format=f,digits=1) [1] 6.7 round(6.6501,1) [1] 6.7 round(6.651,1) [1] 6.7 version

Re: [R] Tools For Preparing Data For Analysis

2007-06-11 Thread Barry Rowlingson
Chris Evans wrote: Thanks Ted, great thread and I'm impressed with EpiData that I've discovered through this. I'd still like something that is even more integrated with R but maybe some day, if EpiData go fully open source as I think they are doing (A full conversion plan to secure this and

Re: [R] Rounding?

2007-06-11 Thread jim holtman
your number 6.6501 is to large to fit in a floating point number. It takes 56 bits and there are only 54 in a real number so the system see it as 6.65 and does the rounding to an even digit; 6.6 6.651 does fit into a real number (takes 54 bits) and this will now round to

Re: [R] Rounding?

2007-06-11 Thread jim holtman
A good reference is: *What Every Computer Scientist Should Know About Floating-Point* *...*http://docs.sun.com/source/806-3568/ncg_goldberg.html http://docs.sun.com/source/806-3568/ncg_goldberg.html On 6/11/07, BXC (Bendix Carstensen) [EMAIL PROTECTED] wrote: I was a bit puzzed by:

[R] simultaneous computing

2007-06-11 Thread Markus Schmidberger
Hello, which possibilities are available in R for simultaneous or parallel computing? I only could find biopara (http://cran.r-project.org/src/contrib/Descriptions/biopara.html) Are there other possibilities? Are there special groups working on simultaneous computing with R? Thanks Markus --

[R] Selecting all values smaller than X in a dataframe

2007-06-11 Thread Christoph Scherber
Dear R users, I have a correlation matrix for a dataframe called synth, for which I now want to select only those cells that have correlations larger than +/-0.6: synth=data.frame(x=rnorm(10,1),y=rnorm(10,2),z=rnorm(10,0.5)) w=cor(synth,use=pairwise.complete.obs) w=as.data.frame(w)

Re: [R] simultaneous computing

2007-06-11 Thread Mike Lawrence
There's RMPI tutorial: http://ace.acadiau.ca/math/ACMMaC/Rmpi/index.html On 11-Jun-07, at 9:11 AM, Markus Schmidberger wrote: Hello, which possibilities are available in R for simultaneous or parallel computing? I only could find biopara

[R] Recoding

2007-06-11 Thread David Scott
I want to do some recoding of variables: code Age into groups and recode a factor into a smaller number of levels. There are a couple of options for recode functions, in the car package and in memisc, and I think in gmisc. Does anyone have any opinions on the the easiest, most reliable

Re: [R] Selecting all values smaller than X in a dataframe

2007-06-11 Thread Vladimir Eremeev
Christoph Scherber-2 wrote: Dear R users, I have a correlation matrix for a dataframe called synth, for which I now want to select only those cells that have correlations larger than +/-0.6: synth=data.frame(x=rnorm(10,1),y=rnorm(10,2),z=rnorm(10,0.5))

[R] Odp: Selecting all values smaller than X in a dataframe

2007-06-11 Thread Petr PIKAL
Hi [EMAIL PROTECTED] napsal dne 11.06.2007 14:09:45: Dear R users, I have a correlation matrix for a dataframe called synth, for which I now want to select only those cells that have correlations larger than +/-0.6: synth=data.frame(x=rnorm(10,1),y=rnorm(10,2),z=rnorm(10,0.5))

Re: [R] simultaneous computing

2007-06-11 Thread Luke Tierney
The package snow available from CRAN is one possibility. Best, luke On Mon, 11 Jun 2007, Markus Schmidberger wrote: Hello, which possibilities are available in R for simultaneous or parallel computing? I only could find biopara

Re: [R] Looking for R-code for non-negative matrix factorization in the presence of Gaussian or Poisson noise

2007-06-11 Thread Thomas Lumley
On Mon, 11 Jun 2007, [EMAIL PROTECTED] wrote: Hi all, Has any of you implemented code for non-negative matrix factorization to solve Y=T P' +E; dim(Y)=n,p ; dim(T)=n,nc; dim (P)=(p,nc); dim(E)=n,p where T and P must be non-negative and E either Gaussian or Poisson noise. I'm looking

Re: [R] Selecting all values smaller than X in a dataframe

2007-06-11 Thread Ingo Michaelis
Christoph Scherber Christoph.Scherber at agr.uni-goettingen.de writes: Dear R users, I have a correlation matrix for a dataframe called synth, for which I now want to select only those cells that have correlations larger than +/-0.6:

Re: [R] Recoding

2007-06-11 Thread David Barron
I've used the one in car with no problems. For grouping a continuous variable you might also consider cut. On 11/06/07, David Scott [EMAIL PROTECTED] wrote: I want to do some recoding of variables: code Age into groups and recode a factor into a smaller number of levels. There are a couple

[R] biplot package II

2007-06-11 Thread terra
Dear all, I've been learning biplot (Gabriel, 1971) and some days ago I sent for this list a procedural function with invitation for a collaborative package. Jari Oksanen made some suggestions and I agree with all. So, I reworked the function under the object-oriented programming (OOP/S3). I

[R] Gini coefficient in R

2007-06-11 Thread MICHELE DE MEO
If I use the Ineq library and the Gini function in this way: Gini(c(100,0,0,0)) I obtain the result 0.75 instead of 1 (that is the perfect inequality). I think Gini's formula in Ineq is based on a formula as reported here: http://mathworld.wolfram.com/GiniCoefficient.html but in the case of

Re: [R] Gini coefficient in R

2007-06-11 Thread Dimitris Rizopoulos
try this version instead: gini - function(x, unbiased = TRUE, na.rm = FALSE){ if (!is.numeric(x)){ warning('x' is not numeric; returning NA) return(NA) } if (!na.rm any(na.ind - is.na(x))) stop('x' contain NAs) if (na.rm) x - x[!na.ind] n -

Re: [R] Rounding?

2007-06-11 Thread Peter Dalgaard
jim holtman wrote: your number 6.6501 is to large to fit in a floating point number. It takes 56 bits and there are only 54 in a real number so the system see it as 6.65 and does the rounding to an even digit; 6.6 6.651 does fit into a real number (takes 54 bits) and

[R] design package

2007-06-11 Thread Qin Qin
Hi all: I tried to install design package for R used under the window. But it dose not work. Here is the message from R. library(Design) Error in library(Design) : 'Design' is not a valid package -- installed 2.0.0? I did not find 2.x version of Design package for window 95. I installed

Re: [R] p-value from GEE why factor 2*pnorm?

2007-06-11 Thread Carmen Meier
I got an answer for the other question (thank you) But there is another question (I am afraid this is a basic question ...) In this tread there is a hint hwo to calculate the p-vlue of an GEE: _http://finzi.psych.upenn.edu/R/Rhelp02a/archive/74150.html_ Then, get the P values using a normal

Re: [R] p-value from GEE why factor 2*pnorm?

2007-06-11 Thread Benilton Carvalho
the recommendation was to use lower.tail=FALSE. b On Jun 11, 2007, at 11:21 AM, Carmen Meier wrote: I got an answer for the other question (thank you) But there is another question (I am afraid this is a basic question ...) In this tread there is a hint hwo to calculate the p-vlue of

Re: [R] p-value from GEE why factor 2*pnorm?

2007-06-11 Thread Carmen Meier
Benilton Carvalho schrieb: the recommendation was to use lower.tail=FALSE. b O but then the results are significant and this does not match the observation. The results are matching the observations if the formula is

[R] lmekin() function in kinship package

2007-06-11 Thread Qiong Yang
Hi, I had a problem with the lmekin() in kinship package: lmekin() can not be wrapped into another function library(kinship) #creat an example dataset xx-rnorm(100) yy-rnorm(100) id-1:100 test.dat-as.data.frame(cbind(xx,yy,id)) rm(xx,yy,id)

[R] autoregressive spectral density estimate by andrews' plug-in method?

2007-06-11 Thread Martin Ivanov
Hello! I would like to ask if there is in R a function that estimates the spectral density function of a stochastic series at frequency zero by the plug-in method, advocated by Andrews in his paper Heteroscedasticity and Autocorrelation Consistent Covariance Matrix Estimation, Econometrica,

Re: [R] open .r files with double-click

2007-06-11 Thread Bert Gunter
However, do note (on Windows) that you can use an external text/programming editors (see CRAN's listings)and can register .r / .R files to open automatically in the chosen editor when clicked on.At least some of these editors (eg TINN-R) can be configured to automatically and simultaneously open

Re: [R] p-value from GEE why factor 2*pnorm?

2007-06-11 Thread Robert A LaBudde
At 11:21 AM 6/11/2007, Carmen wrote: snip In this tread there is a hint hwo to calculate the p-vlue of an GEE: _http://finzi.psych.upenn.edu/R/Rhelp02a/archive/74150.html Then, get the P values using a normal approximation for the distribution of z: / 2 *

Re: [R] Fwd: Using odesolve to produce non-negative solutions

2007-06-11 Thread Jeremy Goldhaber-Fiebert
Hi Spencer, Thank you for your response. I also did not see anything on the lsoda help page which is the reason that I wrote to the list. From your response, I am not sure if I asked my question clearly. I am modeling a group of people (in a variety of health states) moving through time (and

Re: [R] design package

2007-06-11 Thread Uwe Ligges
Which version of R For R-2.5.0, Design (with D, not d) does not pass the checks under Windows. Hence there is no Windows version for R-2.5.0 currently available on CRAN. For older versions of R (down to R-1.7.x): Please either use the Windows binary version, r.g. by typing

Re: [R] p-value from GEE why factor 2*pnorm?

2007-06-11 Thread Benilton Carvalho
Well, AFAIK, the definition of a p-value is the probability of observing something at least as extreme as the observed data. If you observed z, and Z follows a std-normal p-value = P( Z -abs(z) ) + P( Z abs(z) ) = 2*P ( Z abs(z) ) = 2*pnorm(z, lower.tail=FALSE) try z=0 (you should

[R] how to ignore error messages?

2007-06-11 Thread Jianping Jin
Dear group: I wrote a code to iterate a non-linear fit with a set of data. The entire process didn't implemented to the end because an error message, singular gradient. I knew that some sub-sets (columns) do not fit my formula well and may result in parameters going to infinity. It is pretty

[R] history in R.app on Mac

2007-06-11 Thread Weiwei Shi
hi, I am really confused by history function in R.app for Mac: Here is some test from command-line: ls() character(0) ls() character(0) ls() character(0) history() ls() ls() ls() history() That is what I expected. But from R.app, it does not show anything. _ platform

Re: [R] Fwd: Using odesolve to produce non-negative solutions

2007-06-11 Thread Martin Henry H. Stevens
Hi Jeremy, First, setting hmax to a small number could prevent a large step, if you think that is a problem. Second, however, I don't see how you can get a negative population size when using the log trick. I would think that that would prevent completely any negative values of N (i.e.

Re: [R] how to ignore error messages?

2007-06-11 Thread Xiaohui
see ?try Jianping Jin wrote: Dear group: I wrote a code to iterate a non-linear fit with a set of data. The entire process didn't implemented to the end because an error message, singular gradient. I knew that some sub-sets (columns) do not fit my formula well and may result in

Re: [R] how to ignore error messages?

2007-06-11 Thread Dieter Menne
Jianping Jin jjin at email.unc.edu writes: ... I wrote a code to iterate a non-linear fit with a set of data. The entire process didn't implemented to the end because an error message, singular gradient. I knew that some sub-sets (columns) do not fit my formula well and may result in

Re: [R] history in R.app on Mac

2007-06-11 Thread Rob J Goedman
If you click on the history icon in the toolbar in R.app you will see them. Or up and down arrows. R.app implements some extra features, such as multiline command retrieval. Through preference settings you can control if you want to see just a single ls() or multiple, etc. Rob On Jun 11,

Re: [R] Fwd: Using odesolve to produce non-negative solutions

2007-06-11 Thread Spencer Graves
in line Martin Henry H. Stevens wrote: Hi Jeremy, First, setting hmax to a small number could prevent a large step, if you think that is a problem. Second, however, I don't see how you can get a negative population size when using the log trick. SG: Can lsoda estimate complex or imaginary

Re: [R] Fwd: Using odesolve to produce non-negative solutions

2007-06-11 Thread Martin Henry H. Stevens
Hi Spencer, I have copied Woody Setzer. I have no idea whether lsoda can estimate parameters that could take imaginary values. Hank On Jun 11, 2007, at 12:52 PM, Spencer Graves wrote: in line Martin Henry H. Stevens wrote: Hi Jeremy, First, setting hmax to a small number could prevent a

Re: [R] history in R.app on Mac

2007-06-11 Thread Weiwei Shi
got it. Thanks. I was always trying to find it at menu and did not expect that in my toolbar :) I just assumed that any tool icon should also be shown on menu, which is just a routine to me. -w On 6/11/07, Rob J Goedman [EMAIL PROTECTED] wrote: If you click on the history icon in the toolbar

[R] Problem with RSVGTipsDevice

2007-06-11 Thread mister_bluesman
Hi there. I am still trying to get the RSVGTipsDevice to work, yet I can not. I have copied the first example from RSVGTipsDevice documentation: library(RSVGTipsDevice) devSVGTips(C:\\svgplot1.svg, toolTipMode=1, title=SVG example plot 1: shapes and points, tooltips are title + 1 line)

Re: [R] Weighted least squares

2007-06-11 Thread Tim Hesterberg
As John noted, there are different kinds of weights, and different terminology: * inverse-variance weights (accuracy weights) * case weights (frequencies, counts) * sampling weights (selection probability weights) I'll add: * inverse-variance weights, where var(y for observation) = 1/weight

Re: [R] Fwd: Using odesolve to produce non-negative solutions

2007-06-11 Thread Ravi Varadhan
Spencer, Lsoda does not estimate any parameters (nlmeODE does parameter estimation). It just computes the solution trajectory, at discrete times, of a dynamical systems (i.e. set of differential equations). It only works with real numbers, as far as I know. Ravi.

Re: [R] Problem with RSVGTipsDevice

2007-06-11 Thread Mark Difford
Hi Mister_Bluesman, Sorry to sadden you further, but your example svg works perfectly on my machine (R 2.5 running on Vista). All the ToolTips display. Perhaps it's your browser [setup]: I'm using Opera 9.20 (Build 8771). HTH, Mark Difford. mister_bluesman wrote: Hi there. I am still

Re: [R] Problem with RSVGTipsDevice

2007-06-11 Thread Mark Difford
Mister_Bluesman, Perhaps I should have been more precise: your included svgplot1.svg displays fine...! mister_bluesman wrote: Hi there. I am still trying to get the RSVGTipsDevice to work, yet I can not. I have copied the first example from RSVGTipsDevice documentation:

[R] debugger library. Error?

2007-06-11 Thread Jose Quesada
Hi, running the following example code (taken from the docs) will prodice an error in my system (R 2.5.0, win): library(debug) mtrace(glm) # turns tracing on names( tracees) # glm check.for.tracees( package:base) # glm glm(stupid.args) # voila le debugger Error in all.levs[[j]] : subscript out of

Re: [R] Fwd: Using odesolve to produce non-negative solutions

2007-06-11 Thread Ravi Varadhan
Hi Jeremy, A smaller step size may or may not help. If the issue is simply truncation error, that is the error involved in discretizing the differential equations, then a smaller step size would help. If, however, the true solution to the differential equation is negative, for some t, then the

Re: [R] Fwd: Using odesolve to produce non-negative solutions

2007-06-11 Thread Jeremy Goldhaber-Fiebert
Hi Ravi, Thanks for your response. I tried this in Berkeley Madonna and in Matlab. In Berkeley Madonna I did not have the problem (RK4 solver). In Matlab (ode45 solver), I had the problem if I did not use their NonNegative option. My thought was that NonNegative uses something like an

Re: [R] Problem with RSVGTipsDevice

2007-06-11 Thread mister_bluesman
Ah. Now that's intersting. It works in Opera. But do you get an annoying 'Null' label by the cursor when you place it over the svg file? MANY thanks Mark Difford wrote: Mister_Bluesman, Perhaps I should have been more precise: your included svgplot1.svg displays fine...!

Re: [R] Fwd: Using odesolve to produce non-negative solutions

2007-06-11 Thread Setzer . Woodrow
By the way, if someone could forward the original question to me (I'm subscribed to but not currently receiving R-help, as I found I was spending too much time reading it!) I might think of something more useful. (alternatively, when was it posted; I can find it on gmane, too). Woody R. Woodrow

Re: [R] Fwd: Using odesolve to produce non-negative solutions

2007-06-11 Thread Setzer . Woodrow
Hi, all. lsoda can certainly not handle complex parameters. You can try (as Hank suggested) limiting hmax. You can also crank up relative and absolute precision by specifying smaller values of rtol and atol. I've seen similar problems in which the state variable becomes negative, with very

Re: [R] Problem with RSVGTipsDevice

2007-06-11 Thread Mark Difford
Hi Mister_Bluesman, Sadly it is; and it stays up (in what's probably a bold, sans serif), lurking behind the ToolTip, when that comes up. This is not my area of X, but any standard txt editor will open an svg file. If you search through your file you will find lots and lots of null arguments

[R] R vs. Splus in Pharma/Devices Industry

2007-06-11 Thread Cody_Hamilton
Following up to some extent on Friday's discussion regarding the 'validation' of R, could I ask the list group's opinion on possible advantages of R over Splus from a pharma/devices perspective? I wish to exclude the obvious price difference, which doesn’t seem to carry as much weight as I

[R] Reading old S-plus dmp files

2007-06-11 Thread francogrex
I'm sorry that this question has been asked before but I ask it again because in the archives I didn't see a solution. It's an old S-plus dmp file for a hierarchical bayes linear model program written by DuMouchel and available publicly and freely at:

Re: [R] R vs. Splus in Pharma/Devices Industry

2007-06-11 Thread Frank E Harrell Jr
[EMAIL PROTECTED] wrote: Following up to some extent on Friday's discussion regarding the 'validation' of R, could I ask the list group's opinion on possible advantages of R over Splus from a pharma/devices perspective? I wish to exclude the obvious price difference, which doesn’t seem to

Re: [R] Fwd: Using odesolve to produce non-negative solutions

2007-06-11 Thread Ravi Varadhan
Hi Jeremy, A smaller step size may or may not help. If the issue is simply truncation error, that is the error involved in discretizing the differential equations, then a smaller step size would help. If, however, the true solution to the differential equation is negative, for some t, then the

Re: [R] Problem with RSVGTipsDevice

2007-06-11 Thread mister_bluesman
Ah. Now that's intersting. It works in Opera. But do you get an annoying 'Null' label by the cursor when you place it over the svg file? MANY thanks -- View this message in context: http://www.nabble.com/Problem-with-RSVGTipsDevice-tf3902760.html#a11065042 Sent from the R help mailing list

[R] Can I access the filename of the active editor window?

2007-06-11 Thread Paul Johnson
When I run a script from an open editor window (using Ctrl-A, Ctrl-R), I would like the filename of the script to be automatically written into the program output, to keep up with frequent version changes. Is there a way to access the filename (+ path) of the open script (the active one, if there

[R] Rearranging Capture History Data in R

2007-06-11 Thread cox
What code can i use to convert a table like this: Tag#Date 1 1 2 1 3 1 4 1 2 2 4 2 1 3 2 3 4 4 Into one like this: Tag 1 2 3 4 #Date header 1 1 0 0 1 2 1 1 1 0 3 1 0 0

[R] A Question about R

2007-06-11 Thread jwang
Hi Sir/Madam, I'm a researcher in university of Guelph, Canada and now considering using R to do some data analysis. I'm wondering whether there is a library available in R that includes algorithms for archetypal analysis? This is a method quite similar to principal components analysis

Re: [R] Lines in dotchart dotplot ?

2007-06-11 Thread Deepayan Sarkar
On 6/11/07, John Kane [EMAIL PROTECTED] wrote: --- [EMAIL PROTECTED] wrote: On 6/9/07, John Kane [EMAIL PROTECTED] wrote: Is it possible to use dotchart or dotplot and set the lines in such a way that they only extend from the left y-axis to the data point? Yes (sort of) in

Re: [R] Fwd: Using odesolve to produce non-negative solutions

2007-06-11 Thread Ravi Varadhan
Jeremy, You should examine the steady-state solution to your system of equations, by setting the time-derivatives to zero and then solving/analyzing the resulting algebraic equations. This should give you some insights. Let us say you have 3 groups, A,B, and C, with initial conditions:

Re: [R] Rearranging Capture History Data in R

2007-06-11 Thread Benilton Carvalho
date = c(1, 1, 1, 1, 2, 2, 3, 3, 4) tag = c(1, 2, 3, 4, 2, 4, 1, 2, 4) table(factor(tag, levels=1:4), factor(date, levels=1:4)) (not sure how you got Tag 1/Date 4 = 1) On Jun 11, 2007, at 3:11 PM, [EMAIL PROTECTED] wrote: What code can i use to convert a table like this: Tag#Date 1

Re: [R] simultaneous computing

2007-06-11 Thread Ramon Diaz-Uriarte
Dear Markus, You might want to check Rmpi, papply, snow, rpvm, and nws. Best, R. On 6/11/07, Markus Schmidberger [EMAIL PROTECTED] wrote: Hello, which possibilities are available in R for simultaneous or parallel computing? I only could find biopara

[R] selecting characters from a line of text

2007-06-11 Thread Tim Holland
Is there a way in R to select certain characters from a line of text? I have some data that is presently in a large number of text files, and I would like to be able to select elements of each text file (elements are always on the same line, in the same position) and organize them into a table.

[R] ylim settings

2007-06-11 Thread Judith Flores
Hi, I need to know what is the plus/minus adjustment (proportion) that ylim applies if I have something like this: ylim=range(c(x,z)) meaning what's the x-(proportion) and the z+ (proportion)? Thank you, Judith

[R] Package update announcements

2007-06-11 Thread Martin Henry H. Stevens
Hi Folks, I was wondering what everyone thought about adding a sentence to each package update announcement that described what the package did. R extensions are so numerous that it is difficult to keep up with them. Would it be appropriate to ask package developers to add a brief sentence

[R] Overlaying lattice graphs

2007-06-11 Thread Seb
Hello I apologize in advance if this question has already be posted on the list, although I could not find a relevant thread in the archives. I would like to overlay xyplots using different datasets for each plot. I typically work on the following data.frame (mydata) structure mydata

Re: [R] ylim settings

2007-06-11 Thread Marc Schwartz
On Mon, 2007-06-11 at 13:57 -0700, Judith Flores wrote: Hi, I need to know what is the plus/minus adjustment (proportion) that ylim applies if I have something like this: ylim=range(c(x,z)) meaning what's the x-(proportion) and the z+ (proportion)? Thank you, Judith If

Re: [R] Overlaying lattice graphs

2007-06-11 Thread Deepayan Sarkar
On 6/11/07, Seb [EMAIL PROTECTED] wrote: Hello I apologize in advance if this question has already be posted on the list, although I could not find a relevant thread in the archives. I would like to overlay xyplots using different datasets for each plot. I typically work on the following

[R] Rgobbi and colours question

2007-06-11 Thread mister_bluesman
Hi I am using R to implement a multidimensional algorithm which maps places based on the distances between each other. This is presented really well in Rggobi. However, for each place plotted, I have another numerical statistic between 0 and 1 which I would like to represent by colouring the

Re: [R] selecting characters from a line of text

2007-06-11 Thread Stephen Tucker
Maybe substring() is what you're looking for? Some examples: substring(textstring,1,5) [1] texts substring(textstring,3) [1] xtstring substring(textstring,3,nchar(textstring)) [1] xtstring --- Tim Holland [EMAIL PROTECTED] wrote: Is there a way in R to select certain characters from a line

[R] barplot and map overlay

2007-06-11 Thread Héctor Villalobos
Hi, I wonder if it is possible with the graphics package to overlay one or several plots (barplots, for example) over a map. Data for the map is in a data frame with the latitude and longitude coordinates, and then: plot(map$long, map$lat, type =l) produces the map. I want to put each barplot

[R] Components in the control list of optim()

2007-06-11 Thread Daniel Young
Because the function optim() does not return the values of the components in the 'control' list, I am seeking help to uncover the blackbox for some of these components: parscale: par/parscale is used for ndeps, but how does optim() set parscale at the first place? abstol: no default information

Re: [R] Unable to load RMySQL

2007-06-11 Thread Joe W. Byers
All, I have been down for the past two weeks moving from Tulsa to Houston for a new job. I just got internet access yesterday at our new home. I will be working on the RMySQL binary later this week if possible. There is a hitch in compiling the binary. I have the binaries for an older

Re: [R] Textpad help

2007-06-11 Thread Mary Royerr
Thx for the response John and Jim. I had heard before that textpad is a good editor for R codes and one can run the R code from textpad itself and once we do that, the output is displayed on the textpad itself. The reason I wanted to use textpad is to make use of the syntax highlighting a la SAS.

[R] question about data availale in .RData file using the biobase package

2007-06-11 Thread ramakanth reddy
Hi all I am analyzing micro array data and I have R workspace images as my source of the data(.Rdata format).That was in the biobase package format,so I used some commands from the bio base package manual and could write the data into excel files. The data I am working on is the cancer data.

Re: [R] Textpad help

2007-06-11 Thread jim holtman
I use TextPad for editting but have not tried to use it to submit code. I use tinn-r which does the highlighting of the syntax and does submit to R without any problems. I will take a look at the link and see what happening when I try to use TextPad for submission. On 6/11/07, Mary Royerr

Re: [R] Textpad help

2007-06-11 Thread jim holtman
I looked at the link and it would never work for me since with TextPad it starts a new session with each submission. With tinn-R it just sends it to the current GUI session which makes debugging much easier since if an error occurs you have the GUI to see exactly what was sent. I would strongly

Re: [R] Rounding?

2007-06-11 Thread Prof Brian Ripley
On Mon, 11 Jun 2007, jim holtman wrote: your number 6.6501 is to large to fit in a floating point number. It takes 56 bits and there are only 54 in a real number so the system see it as 6.65 and does the rounding to an even digit; 6.6 I'll take it you mean a IEC60559 double,