Re: [R] Is this correct?

Joe a écrit : Dear userR, With the following results, are they correct or acceptable? What can I do to correct them if they are not correct? http://cran.r-project.org/doc/FAQ/R-FAQ.html#Why-doesn_0027t-R-think-these-numbers-are-equal_003f see also ?round hih

Re: [R] greek symbols using pch

This is now well off the topic of the subject line, but I am afraid some misinformation has been propagated (and that is the `bug'). There _are_ bugs in the code shown: the postscript fonts support 32:255, not 1:256, and pch:0:31 are not taken from the font. It seems an uninformed

Re: [R] Writing to a file with fixed precision

On Mon, 10 Oct 2005, Marc Schwartz wrote: On Mon, 2005-10-10 at 19:50 -0400, Richard Hedger wrote: Hi, I'm trying to ouput to a filled with a fixed precision: eg. if I have data x=c(1.0,1.4,2.0), I want to be able to ouput the following to a file: 1.00 1.40

[R] aligning column of xyplots and removing space between them

The code below displays three graphs in three rows and one column but: 1. I want to remove the space between the graphs (I tried playing with position= arg to print.trellis but it seems quite difficult to get the right values and all my attempts had space between them or had overlapping graphs.

[R] Q: Suggestions for long-term data/program storage policy?

Dear list, we are a statistical/epidemiological departement that - after a few years of rapid growth - finally is getting around to formulate a general data storage and retention policy - mainly to ensure that we can reproduce results from published papers/theses easier in the future, but

[R] knn

Hello, Why do I get doubt (NA) in the factor of test classification even if I fix l (minimum vote)? By setting l, no doubt should be occurred. Look forward to your reply Haleh __ R-help@stat.math.ethz.ch mailing list

Re: [R] problem with lapply(x, subset, ...) and variable select argument

Gabor Grothendieck wrote: The problem is that subset looks into its parent frame but in this case the parent frame is not the environment in tt but the environment in lapply since tt does not call subset directly but rather lapply does. Try this which is similar except we have added the line

Re: [R] plot - no main title and missing abscissa value

Thanks Uwe. The R version I'm using is 2.1.1 on Mac OS 10.3.9 I was going to try and replicate this on a linux system at home but have not had the time so far. I'll let you know how it goes. Iain Uwe Ligges wrote: Iain Gallagher wrote: Hi. Sorry (esp to Uwe for the repeated messages!)

Re: [R] Multiple expressions, when using substitute()

Yes, I did get a very helpful reply from Marc Schwartz. I have had substitute() working in legend(), when the legend argument has length one. The challenge was to find some way to do the equivalent of substitute() when several expressions appear in parallel, as may be required for legend(). The

Re: [R] font=5 (Was: greek symbols using pch)

Thanks for that. Very instructive, and much appreciated. And sorry, yes, I strayed well off the original topic. The Greek symbols come out fine with font=5 in my locale, Locale: LC_CTYPE=en_GB.UTF-8; LC_NUMERIC=C; LC_TIME=en_GB.UTF-8; I was interested in some of the other nice characters, for

Re: [R] problem with lapply(x, subset, ...) and variable select argument

As Gabor said, the issue here is that subset.data.frame() evaluates the value of the `select' argument in the parent.frame(); Thus, if you create a local function within lapply() (or sapply()) it works: tt - function (n) { x - list(data.frame(a = 1, b = 2), data.frame(a = 3, b = 4))

Re: [R] Q: Suggestions for long-term data/program storage policy?

Alexander Ploner wrote: Dear list, we are a statistical/epidemiological departement that - after a few years of rapid growth - finally is getting around to formulate a general data storage and retention policy - mainly to ensure that we can reproduce results from published

Re: [R] problem with lapply(x, subset, ...) and variable select argument

Dimitris Rizopoulos [EMAIL PROTECTED] writes: As Gabor said, the issue here is that subset.data.frame() evaluates the value of the `select' argument in the parent.frame(); Thus, if you create a local function within lapply() (or sapply()) it works: It's more complicated than that: It

Re: [R] knn

Please do read the posting guide and supply a reproducible example as it asks. Here it matters critically what you mean by `fix l' (to what value?). Is this knn in package class? Have you read the help page? Have you read the references (especially the first)? On Tue, 11 Oct 2005, Haleh

Re: [R] Is this correct?

Joe wrote: Dear userR, With the following results, are they correct or acceptable? Yes. See the FAQ: 7.31 Why doesn't R think these numbers are equal? The only numbers that can be represented exactly in R's numeric type are integers and fractions whose denominator is a power of 2. Other

Re: [R] Q: Suggestions for long-term data/program storage policy?

Alexander Ploner wrote: Dear list, we are a statistical/epidemiological departement that - after a few years of rapid growth - finally is getting around to formulate a general data storage and retention policy - mainly to ensure that we can reproduce results from published

Re: [R] R on a supercomputer

On 10/10/05 3:54 PM, Kimpel, Mark William [EMAIL PROTECTED] wrote: I am using R with Bioconductor to perform analyses on large datasets using bootstrap methods. In an attempt to speed up my work, I have inquired about using our local supercomputer and asked the administrator if he thought R

Re: [R] Q: Suggestions for long-term data/program storage policy?

On Tue, 11 Oct 2005, Alexander Ploner wrote: we are a statistical/epidemiological departement that - after a few years of rapid growth - finally is getting around to formulate a general data storage and retention policy - mainly to ensure that we can reproduce results from published

Re: [R] Q: Suggestions for long-term data/program storage policy?

On 10/11/05 6:54 AM, Duncan Murdoch [EMAIL PROTECTED] wrote: Alexander Ploner wrote: Dear list, we are a statistical/epidemiological departement that - after a few years of rapid growth - finally is getting around to formulate a general data storage and retention policy - mainly to ensure

[R] random effects are mixture of normals

Dear All, I wonder if there is an R package to estimate the generalized linear mixed models but with a random effects having a mixture of normals as a prior distributinon .. Thank you, Abderrahim [[alternative HTML version deleted]] __

Re: [R] Writing to a file with fixed precision

On Tue, 2005-10-11 at 08:42 +0100, Prof Brian Ripley wrote: On Mon, 10 Oct 2005, Marc Schwartz wrote: On Mon, 2005-10-10 at 19:50 -0400, Richard Hedger wrote: Hi, I'm trying to ouput to a filled with a fixed precision: eg. if I have data x=c(1.0,1.4,2.0), I want to be able to ouput the

Re: [R] Vectorizing loop

I don't have anything specific to say. Only the following suggestion: Try and find out where in the code most of the time is spent. R has nice tools for that. See ?Rprof. If it's runifpoint() and Kest() that are taking most of the computing time, you may not be able to do much better. Andy

Re: [R] glm contrasts (was: no subject)

estimable in bundle gregmisc (package gmodels) should do this. (Kiebitzers: Hope I got the bundle/package/library definition correct) FWIW, we tried gregmisc as a bundle, but it proved to be too much of a pain, so all of the component packages are now separately provided. The 'gremgisc'

[R] Hmisc latex function

I'm using R 2.2.0 on an up-to-date version of Fedora Core 4 with the latest version of Hmisc. When I run an example from the latex function I get the following: x - matrix(1:6, nrow=2, dimnames=list(c('a','b'),c('c','d','enLine 2'))) x c d enLine 2 a 1 35 b 2 46 latex(x) #

Re: [R] problem with lapply(x, subset, ...) and variable select argument

On Tue, 11 Oct 2005, joerg van den hoff wrote: many thanks to thomas and gabor for their help. both solutions solve my problem perfectly. but just as an attempt to improve my understanding of the inner workings of R (similar problems are sure to come up ...) two more question: 1. why

Re: [R] problem with lapply(x, subset, ...) and variable select argument

Just one simple shortening of DR's solution: tt - function (n) { x - list(data.frame(a=1,b=2), data.frame(a=3,b=4)) print(sapply(x, function(...) subset(...), select = n)) } n - b tt(a) On 10/11/05, Dimitris Rizopoulos [EMAIL PROTECTED] wrote: As Gabor said, the issue here is that

[R] probs in installing packages with R 2.2.0

Dear R users. I was wondering if you could help me with this error message I get when I try to install packages. I just installed R (the latest version, R 2.2.0) on my laptop (windows), and I was trying to install packages intwo ways: 1) Packages Install Package(s)... 2) Packages Install

Re: [R] probs in installing packages with R 2.2.0

I expect you need to set a proxy. This is covered in the rw-FAQ that we do ask you to read before posting. See the item The internet download functions fail. On Tue, 11 Oct 2005, Giannitrapani, Marco GSUK-GSSC wrote: Dear R users. I was wondering if you could help me with this

[R] Sweave and Rnews

Hello list, I am writing a paper for Rnews. I use Sweave to do it. I did not find information about writing a paper for Rnews using Sweave and have some questions: - Is there a problem to use the environment 'Sinput' in the place of 'example' or 'smallexample'. It works fine but perhaps there

[R] Set Covering Problem (SCP) (in mathgraphs)

Dear All, Can I solve a Set Covering Problem (in a mathgraph) using R? Which package should I use for that? Thank you in advance. Sara Maltez Mouro [[alternative HTML version deleted]] __ R-help@stat.math.ethz.ch mailing list

Re: [R] font=5 (Was: greek symbols using pch)

ecatchpole [EMAIL PROTECTED] wrote in message news:[EMAIL PROTECTED] Thanks for that. Very instructive, and much appreciated. And sorry, yes, I strayed well off the original topic. The Greek symbols come out fine with font=5 in my locale, Locale: LC_CTYPE=en_GB.UTF-8; LC_NUMERIC=C;

Re: [R] Sweave and Rnews

Le 11 Octobre 2005 11:18, Stéphane Dray a écrit : - I have some long lines of code in schunk. I did not find any way to cut them in the Rnw file and they appear out of the column in the dvi file. The only solution I found is to cut these lines in the tex file generated by Sweave. Is there a

Re: [R] Q: Suggestions for long-term data/program storage policy?

A general comment. As usual, Brian is right on target. Indeed, this has been written, conferenced, agonized, kvetched, etc. about extensively in the computer science community (and no doubt, among many others ... like accountants). I seem to remember reading a Scientific American Magazine

[R] Two factor (or more) non-parametric comparison of means

Can anyone suggest a good non-parametric test, and an R implementation of the test, that allows for two or more factors? Wilcoxon signed rank allows for only one. Thanks, John John Sorkin M.D., Ph.D. Chief, Biostatistics and Informatics Baltimore VA Medical Center GRECC and University of

[R] Logistic Regression using glm

Hello everyone, I am currently teaching an intermediate stats. course at UCSD Extension using R. We are using Venables and Ripley as the primary text for the course, with Freund Wilson's Statistical Methods as a secondary reference. I recently gave a homework assignment on logistic

Re: [R] Two factor (or more) non-parametric comparison of means

?kruskal.test 2005/10/11, John Sorkin [EMAIL PROTECTED]: Can anyone suggest a good non-parametric test, and an R implementation of the test, that allows for two or more factors? Wilcoxon signed rank allows for only one. Thanks, John John Sorkin M.D., Ph.D. Chief, Biostatistics and

Re: [R] Two factor (or more) non-parametric comparison of means

I suspect that John may be looking for ?friedman.test, which I believe will allow for a two-way non-parametric test. kruskal.test() will perform a non-parametric test on two or more samples (or factor levels) as a generalization of wilcox.test() for one or two samples (or factor levels). HTH,

Re: [R] Hmisc latex function

On Tue, 2005-10-11 at 10:01 -0400, Rick Bilonick wrote: I'm using R 2.2.0 on an up-to-date version of Fedora Core 4 with the latest version of Hmisc. When I run an example from the latex function I get the following: x - matrix(1:6, nrow=2, dimnames=list(c('a','b'),c('c','d','enLine 2')))

[R] Problems with plot function

Hello all R users, My simulation function works correctly, but I have problems with plot function. You will find the following code using it. Thank you for your help ## simulation - function(k, n){ conc - seq(0,10,by=0.5) #choixg - seq(1,

Re: [R] Logistic Regression using glm

You're fitting two different models. The latter is saying: logit(p)=x+e, where e is a normal error, so that logit(p) is normal. lm fits a Linear Model, which uses normal error. The former says that p is Bernoulli; and p~Bernoulli does not imply logit(p) is normal. A Generalized Linear Model has

Re: [R] Logistic Regression using glm

One of these is modelling the mean of the logit of p, the other is modelling the logit of the mean of p. They aren't the same. -thomas On Tue, 11 Oct 2005, Daniel Pick wrote: Hello everyone, I am currently teaching an intermediate stats. course at UCSD Extension using R. We are

[R] Any way to add to data frame saved as .rData file?

Hi all, I've got a script that generates a few moderate-size data frames, and then puts them together into one big data frame at the end in order to write that data frame to disk, so that it may be re-opened later on... I'm trying to trim down memory requirements in this script, so I was

[R] non-zero sequence of numbers

Can anyone think of a way to create a pretty() sequence that excludes zero? Or a way to remove the zero from a sequence after using pretty()? Thanks, - Jason Jason Horn Boston University Department of Biology 5 Cumington Street Boston, MA 02215 [EMAIL PROTECTED] office: 617 353 6987 cell:

Re: [R] Two factor (or more) non-parametric comparison of means

Marc Schwartz (via MN) [EMAIL PROTECTED] writes: I suspect that John may be looking for ?friedman.test, which I believe will allow for a two-way non-parametric test. You could be right, but you could also be wrong... It depends quite a bit on what is meant by a two-factor comparison of means.

Re: [R] Two factor (or more) non-parametric comparison of means

Peter, Good points all around. Hopefully John will provide further clarification. Best regards, Marc On Tue, 2005-10-11 at 19:46 +0200, Peter Dalgaard wrote: Marc Schwartz (via MN) [EMAIL PROTECTED] writes: I suspect that John may be looking for ?friedman.test, which I believe will

Re: [R] non-zero sequence of numbers

Jason Horn [EMAIL PROTECTED] writes: Can anyone think of a way to create a pretty() sequence that excludes zero? Or a way to remove the zero from a sequence after using pretty()? The former is rather hard because zero is generally considered just about the prettiest number around... As for

[R] Need help write a function

Dear all, I am still learning R with write a small function for my self. I was wondering if someone can help me to write a R function formula below: Z_k (x) = \sum_{i=0}^{i=k} \binom{n}{i} (m-1)^i Thanks a million in advance, Sincerely, Jan Sabee [[alternative HTML version deleted]]

[R] Reading # in file with read.csv

I'm using read.csv to read in a csv file containing '#' characters. For example, say I'm reading the following file (test.csv): var1,var2,var3 a,b,c d,e#,f g,h,i It outputs: read.csv(Raw Data\\test.csv) var1 var2 var3 1abc 2de 3ghi Warning message:

[R] How to get aspect ratio as output from from plot()

Is there a way to get the aspect ratio as output from a plot() call or something similar in the base graphics system? I would like to note vertical exaggeration on an elevation profile. Thanks, Scott Waichler Pacific Northwest National Laboratory [EMAIL PROTECTED]

[R] a problem in random forest

Hi, there: I spent some time on this but I think I really cannot figure it out, maybe I missed something here: my data looks like this: dim(trn3) [1] 7361 209 dim(val3) [1] 7427 209 mg.rf2-randomForest(x=trn3[,1:208], y=trn3[,209], data=trn3, xtest=val3[, 1:208], ytest=val3[,209],

Re: [R] Reading # in file with read.csv

Michel Friesenhahn wrote: I'm using read.csv to read in a csv file containing '#' characters. For example, say I'm reading the following file (test.csv): var1,var2,var3 a,b,c d,e#,f g,h,i It outputs: read.csv(Raw Data\\test.csv) var1 var2 var3 1abc 2de

Re: [R] a problem in random forest

I am sorry for bother. I think I figured that out. the result of votes for test data is not rf$votes, but rf$test$votes On 10/11/05, Weiwei Shi [EMAIL PROTECTED] wrote: Hi, there: I spent some time on this but I think I really cannot figure it out, maybe I missed something here: my data

Re: [R] Need help write a function

Jan Sabee wrote: Dear all, I am still learning R with write a small function for my self. I was wondering if someone can help me to write a R function formula below: Z_k (x) = \sum_{i=0}^{i=k} \binom{n}{i} (m-1)^i Thanks a million in advance, Sincerely, Jan Sabee (This smells like a

[R] Reading # in file with read.csv

Never mind on my previous question below. read.csv(Raw Data\\test.csv,comment.char=) does it. Mike I'm using read.csv to read in a csv file containing '#' characters. For example, say I'm reading the following file (test.csv): var1,var2,var3 a,b,c d,e#,f g,h,i It outputs: read.csv(Raw

Re: [R] How to get aspect ratio as output from from plot()

Waichler, Scott R wrote: Is there a way to get the aspect ratio as output from a plot() call or something similar in the base graphics system? I would like to note vertical exaggeration on an elevation profile. Thanks, Scott Waichler Pacific Northwest National Laboratory [EMAIL

Re: [R] permutational Kolmogorov-Smirnov p-value for paired data

Here is one way to do a paired permutation test: perm1 - function(x,y){ rb - rbinom(length(x),1,0.5) xp - ifelse(rb==1, x, y) yp - ifelse(rb==1, y, x) ks.test(xp,yp)$statistic } my.x - rnorm(100) my.y - my.x + rnorm(100, 0.2, 0.1) mystat -

Re: [R] Reading # in file with read.csv

Michel Friesenhahn [EMAIL PROTECTED] writes: I'm using read.csv to read in a csv file containing '#' characters. For example, say I'm reading the following file (test.csv): var1,var2,var3 a,b,c d,e#,f g,h,i It outputs: read.csv(Raw Data\\test.csv) var1 var2 var3 1ab

Re: [R] How to get aspect ratio as output from from plot()

Sundar, Perhaps this will work for you? plot(1:10) w - par(pin)[1]/diff(par(usr)[1:2]) h - par(pin)[2]/diff(par(usr)[3:4]) asp - w/h Thank you for your help. For vertical exaggeration I will make a slight change to make it more intuitive (for me): w - diff(par(usr)[1:2]) / par(pin)[1]

[R] problems with levelplot and contourplot

Hello, Using the following code i want to make a level or contourplot of some data that I produced library(grid);library(lattice); mydata - read.table(avgee.dat); mymat - as.matrix(mydata); mymat -t(mymat) vals-as.vector(mymat); conc-c(0.0, 0.5, 1, 1.5, 2, 2.5, 3, 3.5, 4, 4.5, 5); a- c(0.01,

Re: [R] problems with levelplot and contourplot

Tom Lenaerts wrote: Hello, Using the following code i want to make a level or contourplot of some data that I produced library(grid);library(lattice); No need to explicitly load grid. This is done when attaching lattice. Also, you do not need any ; at the end of any lines. mydata

[R] an error in my using of nnet

Hi, there: I am trying nnet as followed: mg.nnet-nnet(x=trn3[,r.v[1:100]], y=trn3[,209], size=5, decay = 5e-4, maxit = 200) # weights: 511 initial value 13822.108453 iter 10 value 7408.169201 iter 20 value 7362.201934 iter 30 value 7361.669408 iter 40 value 7361.294379 iter 50 value 7361.045190

Re: [R] problems with levelplot and contourplot

Dear Sundar Thanks a lot, this resolves the problem of the warning messages. Yet it does not produce the plot: I only get a blank Quartz panel (i'm using R on Mac OS X 1.3). I attached the file and the corrected code. Maybe it is a problem of the data? library(lattice); mydata -

Re: [R] Any way to add to data frame saved as .rData file?

Ken Termiso wrote: Hi all, I've got a script that generates a few moderate-size data frames, and then puts them together into one big data frame at the end in order to write that data frame to disk, so that it may be re-opened later on... I'm trying to trim down memory requirements in

Re: [R] Q: Suggestions for long-term data/program storage policy?

Berton Gunter wrote: A general comment. As usual, Brian is right on target. Indeed, this has been written, conferenced, agonized, kvetched, etc. about extensively in the computer science community (and no doubt, among many others ... like accountants). I seem to remember reading a

Re: [R] aligning column of xyplots and removing space between them

On 10/11/05, Gabor Grothendieck [EMAIL PROTECTED] wrote: The code below displays three graphs in three rows and one column but: 1. I want to remove the space between the graphs (I tried playing with position= arg to print.trellis but it seems quite difficult to get the right values and all

Re: [R] problems with levelplot and contourplot

I have to correct my previous reply. The changes you propose work The strange thing is that if I upload the R-script to produce the plot from file using source(/Users/tomlenaerts/Programming/ObjC/Tools/AllEE/build/EEcont.R) The script does not work. When I copy-past the script myself into

Re: [R] problems with levelplot and contourplot

Tom Lenaerts wrote: I have to correct my previous reply. The changes you propose work The strange thing is that if I upload the R-script to produce the plot from file using source(/Users/tomlenaerts/Programming/ObjC/Tools/AllEE/build/EEcont.R) The script does not work. When I

Re: [R] non-zero sequence of numbers

Jason Horn [EMAIL PROTECTED] wrote: Can anyone think of a way to create a pretty() sequence that excludes zero? You could use except(pretty(x), 0), if you first defined the (quite useful) set-operation function: R except - function(a,b) unique(a[!match(a, b, 0)]) (Consider this a plug to

Re: [R] Problems with plot function

KOITA Lassana - STAC/ACE wrote: Hello all R users, My simulation function works correctly, but I have problems with plot function. You will find the following code using it. Thank you for your help ## simulation - function(k, n){

Re: [R] an error in my using of nnet

Weiwei Shi wrote: Hi, there: I am trying nnet as followed: mg.nnet-nnet(x=trn3[,r.v[1:100]], y=trn3[,209], size=5, decay = 5e-4, maxit = 200) # weights: 511 initial value 13822.108453 iter 10 value 7408.169201 iter 20 value 7362.201934 iter 30 value 7361.669408 iter 40 value

Re: [R] Under-dispersion - a stats question?

Martin Henry H. Stevens wrote: Hello all: I frequently have glm models in which the residual variance is much lower than the residual degrees of freedom (e.g. Res.Dev=30.5, Res.DF = 82). Is it appropriate for me to use a quasipoisson error distribution and test it with an F

Re: [R] Under-dispersion - a stats question?

On Mon, 10 Oct 2005, Martin Henry H. Stevens wrote: Hello all: I frequently have glm models in which the residual variance is much lower than the residual degrees of freedom (e.g. Res.Dev=30.5, Res.DF = 82). Is it appropriate for me to use a quasipoisson error distribution and test it with

Re: [R] an error in my using of nnet

Hi, I repost my question following the guide: trn3 is my dataset: str(trn3) `data.frame': 7361 obs. of 209 variables: $ V1 : num -0.931 -0.955 -0.843 -0.963 -0.967 ... $ V2 : num -0.4343 -0.0101 0.0707 -0.0909 -1. ... $ V3 : num -0.979 -0.979 -0.979 -0.980 -0.979 ... $ V4 : num -0.987 -0.987

[R] Sometimes having problems finding a minimum using optim(), optimize(), and nlm() (while searching for noncentral F parameters)

Ken Kelley [EMAIL PROTECTED] wrote I have a problem that I have been unable to determine either the best way to proceed and why the methods I'm trying to use sometimes fail. I'm using the pf() function in an optimization function to find a noncentrality parameter that leads to a specific

[R] adding 1 month to a date

Within an R dataset, I have a date field called “date_”. (The dates are in the format “-MM-DD”, e.g. “1995-12-01”.) How can I add or subtract “1 month” from this date, to get “1996-01-01” or “ “1995-11-01”. - [[alternative

Re: [R] Sweave and Rnews

Hi Stéphane Dray wrote: Hello list, I am writing a paper for Rnews. I use Sweave to do it. I did not find information about writing a paper for Rnews using Sweave and have some questions: - Is there a problem to use the environment 'Sinput' in the place of 'example' or 'smallexample'.

[R] R data import

Hello, I hope all is well. I've gone through and searched just about all the manuals, faqs, contribs available on the data import/export process and have not found any answer to a specific question. Hopefully, I will be able to fall back upon the valuable expertise in mailing list. Here goes:

[R] R and SAS pointer/informat functionality

Hello, I hope all is well. I've gone through and searched just about all the manuals, faqs, contribs available on the data import/export process and have not found any answer to a specific question. Hopefully, I will be able to fall back upon the valuable expertise in mailing list. Here goes:

Re: [R] aligning column of xyplots and removing space between them

On 10/11/05, Gabor Grothendieck [EMAIL PROTECTED] wrote: Thanks. That works although the alignment is still not perfect. I am attaching the saved image in both .png and .emf formats so you can see what I mean. Its not far off but its noticeable. In .emf format a portion of the bounding box

Re: [R] need suggestion about building formual

Have you considered writing a function to do the complex math and then calling nls referring to that function? Consider the following (not tried): expg - function(a., x, G0, R, T){ exp((a.[1]+(a.[2]*x)^(1/2)-G0)/(R*T)) } mynls-nls(formula=y~expg(a.=c(a, b), x=x, G0=G0, R=R,

Re: [R] symbol print-name too long

Have you received a reply or otherwise resolved this issue? If no, have you tried to construct the simplest possible example that comes to your mind. Sometimes in the course of trying that, I often figure out what the error message means. If you try this without getting the answer

Re: [R] R and SAS pointer/informat functionality

On Tue, 2005-10-11 at 18:55 -0500, Ratnendra Sharma wrote: Hello, I hope all is well. I've gone through and searched just about all the manuals, faqs, contribs available on the data import/export process and have not found any answer to a specific question. Hopefully, I will be able to

[R] bug checking

I have observed the following behavior, wondering if it is a bug before I submit a report. I am using the plot function with call: plot(X, Y, col=red, . . . ) where X is an object that inherits from classes 'dates' and 'times' (created with the 'dates' function from package 'chron') and

Re: [R] adding 1 month to a date

On Tue, 2005-10-11 at 16:26 -0700, t c wrote: Within an R dataset, I have a date field called date_. (The dates are in the format -MM-DD, e.g. 1995-12-01.) How can I add or subtract 1 month from this date, to get 1996-01-01 or 1995-11-01. There might be an easier way to do this, but

[R] monte carlo simulation

Dear R user: I wonder if it is possible to run monte carlo simulation with dse2 package(MonteCarloSimulations function) using ordinary differential equation. How do I define the model? Or if there are any functions which can run monte carlo simulation using ordinary differential equation.

Re: [R] bug checking

On Tue, 2005-10-11 at 16:07 -1000, Parlamis Franklin wrote: I have observed the following behavior, wondering if it is a bug before I submit a report. I am using the plot function with call: plot(X, Y, col=red, . . . ) where X is an object that inherits from classes 'dates' and

Re: [R] ML optimization question--unidimensional unfolding scaling

There may be a few problems where ML (or more generally Bayes) fails to give sensible answers, but they are relatively rare. What is your likelihood? How many parameters are you trying to estimate? Are you using constrained or unconstrained optimization? If

Re: [R] log4j

I just got 145 hits from RSiteSearch(debugger). Does this help? spencer graves Omar Lakkis wrote: Is there a log4j, or similar, package for R? __ R-help@stat.math.ethz.ch mailing list

Re: [R] bug checking

## Code was long, so I simplified it by creating vectors from scratch so it would run as is. Putative bug is still evidenced on the x axis discount.factors.dates - seq.dates(from=09/30/2005, to=09/30/2035) rates-seq(4.4, 5.2, by=0.0025); plot(discount.factors.dates[1:length(rates)], rates,

Re: [R] bug checking

Thanks for the code and the clarifications, including the PDF file. Yes, I can replicate the behavior here (R 2.2.0 on FC4) and I am cc:ing Kurt Hornik, who ported chron to R and is the chron package maintainer. It appears that the culprit is the argument 'col = red', which towards the end of

Re: [R] Under-dispersion - a stats question?

On Tue, 2005-10-11 at 17:16 -0400, Kjetil Holuerson wrote: Martin Henry H. Stevens wrote: Hello all: I frequently have glm models in which the residual variance is much lower than the residual degrees of freedom (e.g. Res.Dev=30.5, Res.DF = 82). Is it appropriate for me to use a

Re: [R] font=5 (Was: greek symbols using pch)

Earl, I don't think that's a bug. Try pdf(font5.pdf, onefile=FALSE) and similarly for postscript(). Ted. On 12/10/05 01:23, Earl F. Glynn wrote,: ecatchpole [EMAIL PROTECTED] wrote in message news:[EMAIL PROTECTED] Thanks for that. Very instructive, and much appreciated. And sorry, yes,

Re: [R] log4j

no, there isn't a general logging package. On 10/12/05, Spencer Graves [EMAIL PROTECTED] wrote: I just got 145 hits from RSiteSearch(debugger). Does this help? spencer graves Omar Lakkis wrote: Is there a log4j, or similar, package for R?