Re: [R] Making R packages (Unix)

2003-09-10 Thread Michael Mader
Gattuso, Jean-Pierre wrote: Hi: I have have taken over from a colleague who prepared an R package and failed to build it on Windows. I am doing this with unix as I am a mac user. Below is the output I get when I use the build command: [gattuso:unix/R/CO2.Rcheck] gattuso% R CMD build

[R] Need your help-SOS

2003-09-10 Thread WeiQiang . Li
Hello, I am a newbie in R project and trying to call prcomp(x) of R function using (D)COM server communicate with R in ASP, and encountering the error Runtime error -2147221493(8004000b). Automation Error, Object is static, operation not allowed. Source code is shown as below:

[R] C code for KalmnaLike

2003-09-10 Thread orasi
Hi it is possible to see the C code for the KalmanLike and Kalmansmooth functions with R? Otherwise, without using R, how can I get the code? Thank arianna __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help

[R] geoR variogram problem

2003-09-10 Thread Dave Nys
Dear GeoR-er, If I use the variog function in the latest release of geoR, the first lag is always ignored. For instance, if you read in geodata, calculates the variogram using the variog function and give in a uvec like uvec=seq(0,max,by=2.44), it only starts giving results from distance=4.88 and

[R] Need your help with SJava package on W2K

2003-09-10 Thread Raimondas B.
Dear R expert, I have the problems with running R from Java on Windows 2000. This is my what i get when i run the program: Loading RInterpreter library java.lang.UnsatisfiedLinkError: no RInterpreter in java.library.path I set all variable (environment). I'd like to notice,that Java from R

Re: [R] Need your help-SOS

2003-09-10 Thread Thomas W Blackwell
WeiQiang - As I read it, both difficulties arise on the DCOM side, not in the R syntax. Problem 1, and I'm just guessing, could arise if you are not allowed to overwrite the value of Result in the DCOM environment. Try again, using two different variable names in the two successive lines.

[R] Off Topic: Good reference for sample size calculations

2003-09-10 Thread Warnes, Gregory R
Hi All, This is off topic, but we're drawing a blank here.. In a presentation I'll be giving next week, I want to include a reference to a good general text on computing sample sizes for standard experiments. Can anyone recommend a good book to use for this purpose? Thanks, -Greg

RE: [R] Off Topic: Good reference for sample size calculations

2003-09-10 Thread Harold Doran
Jacob Cohen's book Statistical Power Analysis for the Behavioral Sciences is one. -- Harold C. Doran Director of Research and Evaluation New American Schools 675 N. Washington Street, Suite 220 Alexandria, Virginia 22314 703.647.1628 http://www.edperform.net -Original

Re: [R] Off Topic: Good reference for sample size calculations

2003-09-10 Thread Marc Schwartz
On Wed, 2003-09-10 at 07:37, Warnes, Gregory R wrote: Hi All, This is off topic, but we're drawing a blank here.. In a presentation I'll be giving next week, I want to include a reference to a good general text on computing sample sizes for standard experiments. Can anyone recommend a

Re: [R] C code for KalmnaLike

2003-09-10 Thread Spencer Graves
Did you try www.r-project.org - Download CRAN - {select a local mirror} - Source code for all platforms? From what I hear, your prayers should be answered there (though I have no personnally built links to compiled code since S-Plus 3.3). hope this helps. spencer graves [EMAIL

Re: [R] C code for KalmnaLike

2003-09-10 Thread Paul Gilbert
[EMAIL PROTECTED] wrote: Hi it is possible to see the C code for the KalmanLike and Kalmansmooth functions with R? It is possible to see all the code distributed with R and the R packages on CRAN. (This is why it is sometimes called an open source project.) Kalman smoothing and Kalman

[R] Plot survey data

2003-09-10 Thread TyagiAnupam
I am trying to make plots that take into account survey weights. This a survey of the US population. To start with I want to explore the data using pairs, plot, coplots and lattice. Are there specialized methods that handle survey weights for plotting? Any pointers? Anupam.

Re: [R] geoR variogram problem

2003-09-10 Thread Paulo Justiniano Ribeiro Jr
Dave GHard to tell withoiut see your data. Pleas send me and example code (and data if necessary) and I will chack P.J. On Wed, 10 Sep 2003, Dave Nys wrote: Dear GeoR-er, If I use the variog function in the latest release of geoR, the first lag is always ignored. For instance, if you

[R] sort a matrix on just one column

2003-09-10 Thread Paul Green
How can I sort(decreasing) a matrix on just the first column? For example, I can I get 8 2 7 5 4 1 from 7 5 4 1 8 2 Thanks __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help

Re: [R] Making R packages (Unix)

2003-09-10 Thread Thomas Lumley
On Tue, 9 Sep 2003, Gattuso, Jean-Pierre wrote: Hi: I have have taken over from a colleague who prepared an R package and failed to build it on Windows. I am doing this with unix as I am a mac user. Below is the output I get when I use the build command: This looks like what happens under

Re: [R] sort a matrix on just one column

2003-09-10 Thread Petr Pikal
Hi On 10 Sep 2003 at 10:50, Paul Green wrote: How can I sort(decreasing) a matrix on just the first column? For example, I can I get 8 2 7 5 4 1 from 7 5 4 1 8 2 I am sure in help pages for sort() is a link to order() mat x y [1,] 7 5 [2,] 4 1 [3,]

[R] dataframe subsetting

2003-09-10 Thread Ryan Thomas Moore
I can create a small dataset, x below, and subset out rows based on values of a certain variable. However, on the dataset I'm working on now, latdata below, I get a subscript error. Any advice is appreciated! Ryan Successful: is.data.frame(x) [1] TRUE x X1 X2 X3 1 1 3 5 2 2 4 6

Re: [R] Plot survey data

2003-09-10 Thread John Fox
Dear Anupam, I may be wrong, but I don't think that there's any standard method to use in plotting with case weights. I can think of two approaches, however: (1) If you have a large sample, and if the range of the weights isn't too large, you could sample your observations with probability of

[R] what is set.fit in function predict.lm

2003-09-10 Thread solares
Hi, what es the parameter set.fit in function predict.lm, is set.fit True then i need the standard error How i cant calculate it?. It is the different what? i see the code of predict.lm How i cant see the matemathics formula for the calculation of standard error.

Re: [R] C code for KalmnaLike

2003-09-10 Thread Paul Gilbert
My apologies. I now realize this question was more specifically about the KalmanLike and KalmanSmooth functions distributed in package ts with base, and not generally about code for calculating the Kalman smoother or likelihood from the Kalman filter. Of course, the sources for these functions

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] insert eps into microsft word

2003-09-10 Thread Marc Schwartz
On Wed, 2003-09-10 at 12:12, Karim Elsawy wrote: it seems that word can not read encapsupalted postscripts generated by R I used this command postscript(output.eps,horizontal=F,onefile=TRUE) since onefile=TRUE produces an encapsualted postscript actually what I'm trying to do is to insert

[R] scan() problem

2003-09-10 Thread Paul Bayer
Dear R-helpers, I have to read some large csv-files into R (30 - 100MB). Since reading with read.csv leads to memory exhausted, I tried with scan(), skipping not needed columns by NULL-elements in what. When these skipped elements are quoted strings with commata inside, R interprets each such

[R] PLS LDA

2003-09-10 Thread Christoph Lehmann
Dear R experts I saw and downloaded the fresh pls package for R. Is there any way of using this pls package for PLS discriminant analysis? If not, is there any other package available. I need a way of classifying objects into e.g. two groups, where nbr_observations nbr_variables many thanks for

Re: [R] insert eps into microsft word

2003-09-10 Thread Martin Biuw
You can also add a Windows metafile preview to your eps image using for instance Ghostscript. Martin On Wed, 10 Sep 2003 18:12:00 +0100, Karim Elsawy [EMAIL PROTECTED] wrote: it seems that word can not read encapsupalted postscripts generated by R I used this command

[R] logistic regression for a data set with perfect separation

2003-09-10 Thread Christoph Lehmann
Dear R experts I have the follwoing data V1 V2 1 -5.800 0 2 -4.800 0 3 -2.867 0 4 -0.867 0 5 -0.733 0 6 -1.667 0 7 -0.133 1 8 1.200 1 9 1.333 1 and I want to know, whether V1 can predict V2: of course it can, since there is a perfect

Re: [R] logistic regression for a data set with perfect separation

2003-09-10 Thread David Firth
On Wednesday, Sep 10, 2003, at 18:50 Europe/London, Christoph Lehmann wrote: Dear R experts I have the follwoing data V1 V2 1 -5.800 0 2 -4.800 0 3 -2.867 0 4 -0.867 0 5 -0.733 0 6 -1.667 0 7 -0.133 1 8 1.200 1 9 1.333 1 and I want to

[R] how to calculate Spearman correlation with missing values

2003-09-10 Thread szhan
Hello, there: I got data matix with missing values. I want to calculate any possible pairwise Spearman correlation rho for each column. Is there a function just like cor(x, y, use=complete.obs) for Pearson correlation? Thanks in advance! Josh __

Re: [R] how to calculate Spearman correlation with missing values

2003-09-10 Thread Torsten Hothorn
On Wed, 10 Sep 2003 [EMAIL PROTECTED] wrote: Hello, there: I got data matix with missing values. I want to calculate any possible pairwise Spearman correlation rho for each column. Is there a function just like cor(x, y, use=complete.obs) for Pearson correlation? Thanks in advance!

Re: [R] insert eps into microsft word

2003-09-10 Thread Gavin Simpson
Depending on your version of Word that simply is *not* true! I can't test this here (I now have Word XP/2002) but my colleagues have Word 97 on their college-supplied service (yes, really!) and despite not being able to see a preview, they can quite happily import eps files I have produced

Re: [R] insert eps into microsft word

2003-09-10 Thread Peter Dalgaard BSA
Rafael A. Irizarry [EMAIL PROTECTED] writes: for word documents submitted to picky journals i usually use seomthing like this: bitmap(plot_1.png,width=6,height=6,res=600,pointsize=12,family=Times) on my computer this resutls in quality just as good (to mu eye) as with postscript. you

Re: [R] insert eps into microsft word

2003-09-10 Thread Brett Melbourne
Your problem may also be that you can't get the figure to print from Word? To get Word to print anything other than a blank box for the eps, you will need to install and use a postscript printer driver for your printer. cheers Brett Brett Melbourne, Postdoctoral Fellow Biological Invasions

RE: [R] logistic regression for a data set with perfect separati

2003-09-10 Thread Ted Harding
On 10-Sep-03 Christoph Lehmann wrote: I have the follwoing data V1 V2 1 -5.800 0 2 -4.800 0 3 -2.867 0 4 -0.867 0 5 -0.733 0 6 -1.667 0 7 -0.133 1 8 1.200 1 9 1.333 1 and I want to know, whether V1 can predict V2: of course it

RE: [R] PLS LDA

2003-09-10 Thread Liaw, Andy
Do you mean the pls.pcr package by Prof. Wehrens? This is what I do: o Code the two groups as 0s and 1s (numeric, not factor). o Run PLS as usual. Cases with predicted values 0.5 get classified as 1s, otherwise as 0s. o Note that you need to modify the code inside the mvr()

RE: [R] PLS LDA

2003-09-10 Thread Christoph Lehmann
Hi Andy Great and thanks a lot! Yes, it is the package from Prof. Wehrens. So I just run the PLS like a Logistic Regression, coding the endogenous variable as binary. So no need of specifying a binary-link function (as we have to when using glm)? And yes of course: I need the LVs which give the

[R] coef names in lm

2003-09-10 Thread Jean Eid
Dear all, I am interested in finding out how to change the names of coefficients in the lm function. I have a design matrix which I called design where each variate has its own name. However when I issue the command: lm.1-lm(response~design-1, weights=some.weights) and follow it with:

RE: [R] coef names in lm

2003-09-10 Thread Liaw, Andy
If you coerce design into a data frame and then do lm.1 - lm(response ~ . - 1, data=design, ...) that should work. Andy -Original Message- From: Jean Eid [mailto:[EMAIL PROTECTED] Sent: Wednesday, September 10, 2003 4:28 PM To: [EMAIL PROTECTED] Subject: [R] coef names in lm

Re: [R] coef names in lm

2003-09-10 Thread Roger D. Peng
Maybe this will work for you: df - as.data.frame(design) lm.2 - lm(response ~ ., df) -roger Jean Eid wrote: Dear all, I am interested in finding out how to change the names of coefficients in the lm function. I have a design matrix which I called design where each variate has its own name.

[R] regression questions

2003-09-10 Thread Paul, David A
I have been puzzling over how to fit some fixed effects models to a set of data. My response function is response - function(a, b, c, alpha1, alpha2, indicator, t, t2) { z = a + b * (t) * exp(-alpha1 * t) + indicator *c * (t2) * exp(-alpha2 * t2) } where

Re: [R] Computing a CDF or many quantiles

2003-09-10 Thread Kevin S. Van Horn
Your method looks like a naive reimplementation of integration, and won't work so well for distributions that have the great majority of the probability mass concentrated in a small fraction of the sample space. I was hoping for something that would retain the adaptability of integrate().

Re: [R] Computing a CDF or many quantiles

2003-09-10 Thread Jerome Asselin
On September 10, 2003 04:03 pm, Kevin S. Van Horn wrote: Your method looks like a naive reimplementation of integration, and won't work so well for distributions that have the great majority of the probability mass concentrated in a small fraction of the sample space. I was hoping for

Re: [R] Computing a CDF or many quantiles

2003-09-10 Thread Jerome Asselin
Also look at ecdf() from package stepfun. HTH, Jerome Asselin __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help

[R] Customised legend in lattice

2003-09-10 Thread Alexander . Herr
Hi List, Am trying to customize a legend in trellis: Draws 2x5 lines in 5 colors and 2 linetypes. I would like to add two more items to the legend showing the key for the line types above the colored legend. Any suggestions welcome - thanks Herry # #Following example

Re: [R] scan() problem

2003-09-10 Thread Gabor Grothendieck
If the records are always of the form: number,...,...,number where ... may contain commas but not double quotes then here is a kludgy solution. Perhaps its sufficient? # scan in data using as the delimiter and keep first and last fields s -