Re: [R] fitdistr, mle's and gamma distribution

2003-09-30 Thread Ben Bolker
Spencer Graves's suggestion of using shape and scale parameters on a log scale is a good one. To do specifically what you want (check values for which the objective function is called and see what happens) you can do the following (untested!), which makes a local copy of dgamma that you

Re: [R] fitdistr, mle's and gamma distribution

2003-09-30 Thread Ben Bolker
PS. 11 MILLION entries?? On Tue, 30 Sep 2003, Ben Bolker wrote: Spencer Graves's suggestion of using shape and scale parameters on a log scale is a good one. To do specifically what you want (check values for which the objective function is called and see what happens) you can

Re: [R] truncated multivariate normal

2003-09-30 Thread Ben Bolker
Do you mean you want to generate deviates from this distribution or that you want the density function (or distribution function)? mvrnorm in package MASS (library(MASS)) will generate multivariate normal deviates. A brute-force approach would just keep picking values until you get ones in

Re: [R] Help needed: plotting with no device

2003-09-29 Thread Ben Bolker
Can you use save.image() to rescue your results? I would try save.image(file=salvage.RData) and see if the file appears. Otherwise I would say you're probably out of luck. On Mon, 29 Sep 2003, Ross Boylan wrote: I have a hung session I would very much like to recover, since it has

Re: [R] Std. errors of intercept and slope

2003-09-26 Thread Ben Bolker
textbook or reference material ... Ben Bolker On Fri, 26 Sep 2003, Yao, Minghua wrote: Dear all, I have the following output generated by linear regression. Since there is only one regression intercept and one slope for one set of data, what is the meaning of std. error for intercept

RE: [R] Std. errors of intercept and slope

2003-09-26 Thread Ben Bolker
I'm afraid you're going to have to look it up in a basic statistics textbook. Ben Bolker On Fri, 26 Sep 2003, Yao, Minghua wrote: Thanks, Ben. Could you tell me the formula for calculating this sd., given (x_i, y_i) (i=1,2,...,N)? We only have one intercept and slope for them

Re: [R] discretization method

2003-09-23 Thread Ben Bolker
Consider: ?cut ?round With more detail we might be able to help more ... On Tue, 23 Sep 2003, Jaime Lopez Carvajal wrote: Hi R users I need to apply discretization to my continuous data. Is there a method in R to do this? Thanks in advance, Jaime -- 620B Bartram Hall

Re: [R] How to extract data from Excel

2003-09-23 Thread Ben Bolker
Your best bet is saving as a comma-separated value file (.csv) and using read.csv to get the data into R. Ben On Tue, 23 Sep 2003 [EMAIL PROTECTED] wrote: Hi, I would like to know how to extract the data from Excel Spreadsheet. Thank you very much. Melissa

Re: [R] 3D plot/surface rotation

2003-09-17 Thread Ben Bolker
in a web browser. Not as good as x/ggobi for interaction with data, but possibly OK for display. Ben Bolker On Wed, 17 Sep 2003, Brian D Ripley wrote: On Wed, 17 Sep 2003, Thomas W Blackwell wrote: Perhaps MZodet wants the interactive, mouse controlled rotation capability offered by ggobi

RE: [R] Retrieve ... argument values

2003-09-16 Thread Ben Bolker
Yes, although this becomes tedious if (e.g.) you have a function that calls two different functions, each of which has many arguments (e.g. plot() and barplot(); then you have to set up a whole lot of arguments that default to NULL and, more annoyingly, you have to document them all in

Re: [R] convert a Character-string to a number

2003-09-12 Thread Ben Bolker
?as.numeric On Fri, 12 Sep 2003, Michael Kirschbaum wrote: Hi Everyone. I have a simple problem but don't know, how to get along. how can I convert the vector a-c(0,01,1,00) in a vector b-c(0.01,1.00) Thank you for suggestions M.Kirschbaum [[alternative HTML

Re: [R] insert eps into microsft word

2003-09-10 Thread Ben Bolker
If you're working in windows, the WMF (windows metafile) format is probably your best bet; it's a vector format like PostScript. On Wed, 10 Sep 2003, Karim Elsawy wrote: it seems that word can not read encapsupalted postscripts generated by R I used this command

Re: [R] Combining the components of a character vector

2003-04-04 Thread Ben Bolker
paste(x,collapse= ) On Wed, 2 Apr 2003, John Miyamoto wrote: Dear Help, Suppose I have a character vector. x - c(Bob, loves, Sally) I want to combine it into a single string: Bob loves Sally . paste(x) yields: paste(x) [1] Bob loves Sally The following function combines

RE: [R] Does R have an inverse wishart distribution?

2003-04-01 Thread Ben Bolker
be able to do some QA work for you. :) -Kevin -Original Message- From: Ben Bolker [mailto:[EMAIL PROTECTED] Sent: Monday, March 31, 2003 3:25 PM To: Kevin Karty Cc: [EMAIL PROTECTED] Subject: Re: [R] Does R have an inverse wishart distribution? Do you want

Re: [R] 2 scatter plots in the same graphic window....

2003-03-26 Thread Ben Bolker
Either n1 - length(x1) n2 - length(x2) plot(c(x1,x2),c(y1,y2),col=rep(c(1,2),c(n1,n2)), pch=rep(c(1,2),c(n1,n2)) OR plot(x1,y1,xlim=range(c(x1,x2)),ylim=range(c(y1,y2)),col=1,pch=1) points(x2,y2,col=2,pch=2) will work. Ben Bolker On Wed, 26 Mar 2003, Vasudevan, Geetha wrote: Hello

Re: [R] Problem with read.table

2003-03-21 Thread Ben Bolker
Should that be sep=\t ? On Fri, 21 Mar 2003, sabrina servanty wrote: Dear all, I was used to work on R1.6 and I have now passed on R1.6.2 but I can't read my file (and that is a big problem!!). I made a data sheet with some spreadsheet in Excell, and save it as separeted by tab .txt.

Re: [R] Error bars for interaction plots

2003-02-28 Thread Ben Bolker
Put your own solution together with arrows(...,angle=90), or use plotCI within gregmisc, or get plotCI from my website (http://www.zoo.ufl.edu/bolker/R/windows -- you want the bbmisc bundle) Ben On Fri, 28 Feb 2003, Katrina Grech wrote: How do I add error bars to an interaction plot of

Re: [R] outliers/interval data extraction

2003-02-20 Thread Ben Bolker
II. How to extract only those values from vector which fulfill the condition of interval (higher than A, and lower than B)? x[xA xB] Rado Bonk __ [EMAIL PROTECTED] mailing list http://www.stat.math.ethz.ch/mailman/listinfo/r-help --

Re: [R] Help in separate window under X11

2003-02-19 Thread Ben Bolker
Your best bet is help.start(). On Wed, 19 Feb 2003, Alexey Shipunov wrote: Dear R users, Is there the possibitily in R under X11 to get (after typing help(...) command) separate help window, as it is in Windows version? Best wishes, = Dr. Alexey

Re: [R] glm and overdispersion

2003-02-18 Thread Ben Bolker
Use family=quasibinomial in your glm() statement. R does what you say rather than doing what you mean: it specifically says something like (Dispersion for binomial family taken to be 1), which is an indication that R is not using a dispersion factor 1 in this case. (There is a hint of this

Re: [R] My remark on libraries

2003-02-07 Thread Ben Bolker
person for updating (and for keeping it coherent ...) Ben Bolker On Fri, 7 Feb 2003, Eric Lecoutre wrote: Hi R community, These days, I am writing some functions to work with 2-ways frequency tables ; you know all this tuff about measures of association: Chisq and derived (phi

Re: [R] Plotting

2003-01-30 Thread Ben Bolker
Depending on what you mean (same x-axis and y-axis), look at ?points or ?lines OR (same x-axis, different y-axes) plot(...); par(new=TRUE); plot(...,ann=FALSE). (The latter is asked frequently on this list.) Ben Bolker On Thu, 30 Jan 2003, Francisco do Nascimento Junior wrote: Hi

Re: [R] Double Exponential

2003-01-30 Thread Ben Bolker
What do you mean by double exponential? Can you give us another more familiar name (e.g. Laplacian), or a reference, or a more detailed description, or a formula for the cumulative or probability density function? Ben Bolker On Thu, 30 Jan 2003, Jason Parcon wrote: To whom it may concern

Re: [R] reading non-existent files

2003-01-28 Thread Ben Bolker
?file.exists On Tue, 28 Jan 2003 [EMAIL PROTECTED] wrote: Dear R-experts I would like to read all files from a directory, the files have names myname0001.txt etc. I paste the directory plus file names and use read.delim(). My problem is that some file names are missing, so I get an

Re: [R] Profile on optim/nlm

2003-01-23 Thread Ben Bolker
I have some code to do this; I keep meaning to clean it up and put it on CRAN, but it's fairly hard work to do it right. I also thought I heard a rumor that there would be something like that in MASS4. There should be a version on http://www.zoo.ufl.edu/bolker/R (subdirectories src and

Re: [R] How to change the label position in axis() ?

2003-01-14 Thread Ben Bolker
See ?par, the source of all wisdom: par(las=1) On Tue, 14 Jan 2003 [EMAIL PROTECTED] wrote: Dear R People, I'm working on a plot function to produce the graph shown below. One of the features my supervisor does not like is that the labels in the Y- axix are shown vertically. Is there

Re: [R] Help on R commands

2002-12-20 Thread Ben Bolker
of these questions are pretty basic; you really need to try to get some help off-list. You may not get too much more help here unless you explain what your situation is (are you doing this for a class, independent project, etc.) and give some indication that you've exhausted local resources ... cheers, Ben

<    1   2   3