Re: [R] The hidden costs of GPL software? - None

2004-11-24 Thread Uwe Ligges
Martin, what about setting up a new mailing list R-hcgs? (acronym for R - The hidden costs of GPL software?) Seems to be worth given the amount of messages in this thread(s). ;-) Uwe __ [EMAIL PROTECTED] mailing list

[R] T-test syntax question

2004-11-24 Thread Vito Ricci
Hi, You did not specify if data are paired or not, as data are paired you should use option paired=TRUE in t.test(). Variances of the two samples have to be not significatevely different, (see ? var.test) to use t.test, if not you should specify var.equal=FALSE. var.equal: a logical variable

[R] Ks.test (was (no subject))

2004-11-24 Thread Vito Ricci
Hi Angela, I believe you should introduce only df as parameters; t distribution as by default mean=0; see this example. x-rt(100,10) ks.test(x, pt,10) One-sample Kolmogorov-Smirnov test data: x D = 0.1414, p-value = 0.03671 alternative hypothesis: two.sided Ciao Vito you wrote:

Re: [R] (no subject)

2004-11-24 Thread Uwe Ligges
[EMAIL PROTECTED] wrote: Good morning, I have to apply the Ks test with the the t distribution. I know I have to write ks.test(data_name,distribution_name, parameters..) but I don't know what is the name fot t distribution and which parameters to introduce? may be mean=0 and freedom degrees in my

[R] Is there a package in R that lets me fit a robust linear mixed model?

2004-11-24 Thread Frank Duan
Dear R people, Happy Thanksgiving! I just wonder if there is a R package that can supply some kind of robust way to fit a linear mixed model. I mean assigning small weights to those observations with large residuals, like iteratively-reweighted-least-squares approach. Many thanks, Frank

Re: [R] The hidden costs of GPL software? - None

2004-11-24 Thread Ted Harding
On 24-Nov-04 John wrote: Off hand, the costs of GPL'd software are not hidden at all. R for instance demands that a would be user sit down and learn the language. This in turn pushes a user into learning more about statistics than the simple overview that Stat 1 presents a student. I'd see

[R] Posting to 2 mailing lists {was How to extract data?}

2004-11-24 Thread Martin Maechler
You posted the referenced message to both R-help and R-sig-finance. Please do *not* post to more than one R-list! Please do *not* post to more than one R-list! Please do *not* post to more than one R-list! Please decide if something is specific for an R-SIG-foo list or if it belongs to

[R] problem with anova and glmmPQL

2004-11-24 Thread Andrew R. Criswell
Hello: I am getting an error message when appplying anova() to two equations estimated using glmmPQL. I did look through the archives but didn't finding anything relevant to my problem. The R-code and results follow. Hope someone can help. ANDREW fm1 -

[R] LDA with previous PCA for dimensionality reduction

2004-11-24 Thread Christoph Lehmann
Dear all, not really a R question but: If I want to check for the classification accuracy of a LDA with previous PCA for dimensionality reduction by means of the LOOCV method: Is it ok to do the PCA on the WHOLE dataset ONCE and then run the LDA with the CV option set to TRUE (runs LOOCV) --

[R] Grumble ...

2004-11-24 Thread Ted Harding
Hi Folks, A Grumble ... The message I just sent to R-help about The hidden costs of GPL ... has evoked a Challenge response: Hi, You´ve just sent a message to [EMAIL PROTECTED] In order to confirm the sent message, please click here This confirmation is necessary because [EMAIL

Re: [R] T-test syntax question

2004-11-24 Thread Adaikalavan Ramasamy
As Vito Ricci has already pointed out, the Welsh test is for two group unpaired data with unequal variance assumption. If you have the original data, say x and y, then you can simply do t.test( x, y, paired=FALSE, var.equal=FALSE ). If you do not have the original data, you can still calculate

[R] Automatic file reading

2004-11-24 Thread Anders Malmberg
Hi, I want to do automatic reading of a number of tables (files) stored in ascii format without having to specify the variable name in R each time. Below is an example of how I would like to use it (I assume files pair1,...,pair8 exist in spec. dire.) for (i in 1:8){ name -

[R] Re: T-test syntax question

2004-11-24 Thread Vito Ricci
Hi, In case of paired data, if you have only differencies and not original data you can get this t test based on differencies: Say d is the vector with differencies data and suppose you wish to test if the mean of differency is equal to zero: md-mean(d) ## sample mean of differencies sdd-sd(d)

[R] Re: CRLF-terminated Fortran files

2004-11-24 Thread Prof Brian Ripley
I added a check for CRLF termination of Fortran and C++ source files to R CMD check and found potential problems in packages BsMD MCMCpack (C++) asypow aws bayesSurv (C++) eha fBasics/fOptions/fSeries gam mclust ncomplete noverlap pan rrcov subselect (C++) survrec I'd be interested to know if C++

Re: [R] Grumble ...

2004-11-24 Thread Duncan Murdoch
On Wed, 24 Nov 2004 10:36:35 - (GMT), (Ted Harding) [EMAIL PROTECTED] wrote: Hi Folks, A Grumble ... The message I just sent to R-help about The hidden costs of GPL ... has evoked a Challenge response: Hi, You´ve just sent a message to [EMAIL PROTECTED] In order to confirm the sent

Re: [R] Automatic file reading

2004-11-24 Thread Adaikalavan Ramasamy
for(i in 1:10){ assign( paste(data, i), i ) } data1 [1] 1 data5 [1] 5 data8 + data5 [1] 13 See help(assign) for more details and examples. On Wed, 2004-11-24 at 11:10, Anders Malmberg wrote: Hi, I want to do automatic reading of a number of tables (files) stored in ascii format

Re: [R] Automatic file reading

2004-11-24 Thread Uwe Ligges
Anders Malmberg wrote: Hi, I want to do automatic reading of a number of tables (files) stored in ascii format without having to specify the variable name in R each time. Below is an example of how I would like to use it (I assume files pair1,...,pair8 exist in spec. dire.) for (i in 1:8){

Re: [R] Automatic file reading

2004-11-24 Thread Arne Henningsen
Hi Andreas, what's about: pair - list() for (i in 1:8){ name - paste(pair,i,sep=) pair[[ i ]] - read.table(paste(/home/andersm/tmp/,name,sep=)) } Arne On Wednesday 24 November 2004 12:10, Anders Malmberg wrote: Hi, I want to do automatic reading of a number of tables (files) stored in

Re: [R] Automatic file reading

2004-11-24 Thread Sean Davis
If you simply want read all files in a given directory, you can do something like: fullpath = /home/andersm/tmp filenames - dir(fullpath,pattern=*) pair - sapply(filenames,function(x) {read.table(paste(fullpath,'/',x,sep=))}) Sorry, untested. But the point is that you can use dir to get all

Re: [R] data.frame into vector

2004-11-24 Thread Jan Goebel
Hi, as other already pointed out as.matrix is what you need. Just one comment: as.matrix(x[1,]) should be much faster for larger data frames compared to as.matrix(x)[1,] Best jan On Tue, 23 Nov 2004, Tiago R Magalhaes wrote: Hi I want to extract a row from a data.frame but I want

Re: [R] timeDate

2004-11-24 Thread Gabor Grothendieck
Gabor Grothendieck ggrothendieck at myway.com writes: : : Yasser El-Zein abu3ammar at gmail.com writes: : : : : : I am looking for up to the millisecond resolution. Is there a package : : that has that? : : : : On Mon, 22 Nov 2004 21:48:20 + (UTC), Gabor Grothendieck : : ggrothendieck at

Re: [R] Grumble ...

2004-11-24 Thread Martin Maechler
Ted == Ted Harding [EMAIL PROTECTED] on Wed, 24 Nov 2004 10:36:35 - (GMT) writes: Ted Hi Folks, A Grumble ... Ted The message I just sent to R-help about The hidden Ted costs of GPL ... has evoked a Challenge response: Ted Hi, You´ve just sent a message to Ted

Re: [R] Grumble ...

2004-11-24 Thread Peter Dalgaard
Duncan Murdoch [EMAIL PROTECTED] writes: On Wed, 24 Nov 2004 10:36:35 - (GMT), (Ted Harding) [EMAIL PROTECTED] wrote: Hi Folks, A Grumble ... The message I just sent to R-help about The hidden costs of GPL ... has evoked a Challenge response: Hi, You´ve just sent a message

Re: [R] LDA with previous PCA for dimensionality reduction

2004-11-24 Thread Ramon Diaz-Uriarte
Dear Cristoph, I guess you want to assess the error rate of a LDA that has been fitted to a set of currently existing training data, and that in the future you will get some new observation(s) for which you want to make a prediction. Then, I'd say that you want to use the second approach. You

Re: [R] Two factor ANOVA in lme

2004-11-24 Thread Fernando Henrique Ferraz P. da Rosa
nat writes: I want to specify a two-factor model in lme, which should be easy? Here's what I have: factor 1 - treatment FIXED (two levels) factor 2 - genotype RANDOM (160 genotypes in total) I need a model that tells me whether the treatment, genotype and interaction terms are

Re: [R] LDA with previous PCA for dimensionality reduction

2004-11-24 Thread Torsten Hothorn
On Wed, 24 Nov 2004, Ramon Diaz-Uriarte wrote: Dear Cristoph, I guess you want to assess the error rate of a LDA that has been fitted to a set of currently existing training data, and that in the future you will get some new observation(s) for which you want to make a prediction. Then, I'd

RE: [R] an R function to search on Prof. Baron's site

2004-11-24 Thread Andy Bunn
Using this function with 2.0.0 XP and Firefox 1.0 (I've rediscovered the internet) produces a curious result. myString - RSiteSearch(string = 'Ripley') myString [1] http://finzi.psych.upenn.edu/cgi-bin/htsearch?config=htdigrun1;restrict=Rhe

Re: [R] LDA with previous PCA for dimensionality reduction

2004-11-24 Thread Christoph Lehmann
Thank you, Torsten; that's what I thought, as long as one does not use the 'class label' as a constraint in the dimension reduction, the procedure is ok. Of course it is computationally more demanding, since for each new (unknown in respect of the class label) observation one has to compute a

[R] Respuesta Automatica CorreoDirect.

2004-11-24 Thread infocd
### Este email esta generado de manera automatica ### Ha escrito usted a una cuenta de correo generica de CorreoDirect. Si desea conocer nuestra política de privacidad, pulse aquí.

[R] scatterplot of 100000 points and pdf file format

2004-11-24 Thread Witold Eryk Wolski
Hi, I want to draw a scatter plot with 1M and more points and save it as pdf. This makes the pdf file large. So i tried to save the file first as png and than convert it to pdf. This looks OK if printed but if viewed e.g. with acrobat as document figure the quality is bad. Anyone knows a way

[R] 2GB dataset

2004-11-24 Thread apollo wong
Hi, do any one have experience with loading dataset that is larger than 2GB into R. My organization is a SAS oriented shop and I'm in the process of switching it to R. One of the complain about R has always been it's inability to handle large dataset (GB) efficiently. I would like some comments

[R] tsdiag for ar?

2004-11-24 Thread Dr Carbon
Is there a way to have the ar function work with tsdiag for on-the-fly visualization of ar fits? I have to fit a great many models of varying order and would like to save the diagnostic graphs. For instance, tsdiag(ar(lh)) tsdiag(arima(lh, order = c(3,0,0))) Thanks...

Re: [R] 2GB dataset

2004-11-24 Thread Prof Brian Ripley
Absolutely no problem on 64-bit OSes with enough memory. Many 32-bit OSes have problems with 2Gb files. Please do read the posting guide and tell us basic facts like which OS you are running on, so we don't have to speculate to answer your question. Also, what you want to do with the

Re: [R] tsdiag for ar?

2004-11-24 Thread Prof Brian Ripley
On Wed, 24 Nov 2004, Dr Carbon wrote: Is there a way to have the ar function work with tsdiag for on-the-fly visualization of ar fits? I have to fit a great many models of varying order and would like to save the diagnostic graphs. First you have to produce them, surely? For instance,

Re: [R] scatterplot of 100000 points and pdf file format

2004-11-24 Thread Marc Schwartz
On Wed, 2004-11-24 at 16:34 +0100, Witold Eryk Wolski wrote: Hi, I want to draw a scatter plot with 1M and more points and save it as pdf. This makes the pdf file large. So i tried to save the file first as png and than convert it to pdf. This looks OK if printed but if viewed e.g. with

RE: [R] scatterplot of 100000 points and pdf file format

2004-11-24 Thread Ted Harding
On 24-Nov-04 Witold Eryk Wolski wrote: Hi, I want to draw a scatter plot with 1M and more points and save it as pdf. This makes the pdf file large. So i tried to save the file first as png and than convert it to pdf. This looks OK if printed but if viewed e.g. with acrobat as document

[R] OOT: frailty-multinivel

2004-11-24 Thread Kjetil Brinchmann Halvorsen
Hola! I started to search for information about multilevel survival models, and found frailty in R. This seems to be something of the same, is it the same? Then: why the name frailty (weekness?) -- Kjetil Halvorsen. Peace is the most effective weapon of mass construction. -- Mahdi

RE: [R] scatterplot of 100000 points and pdf file format

2004-11-24 Thread Liaw, Andy
Marc/Eryk, I have no experience with it, but I believe the hexbin package in BioC was there for this purpose: avoid heavy over-plotting lots of points. You might want to look into that, if you have not done so yet. Best, Andy From: Marc Schwartz On Wed, 2004-11-24 at 16:34 +0100, Witold

Re: [R] scatterplot of 100000 points and pdf file format

2004-11-24 Thread Witold Eryk Wolski
Hi, I tried the ps idea. But I am using pdflatex. You get a even larger size reduction if you convert the ps into a pdf using ps2pdf. But unfortunately there is a quality loss. I have found almost a working solution: a) Save the scatterplot without axes and with par(mar=c(0,0,0,0)) as png . b)

Re: [R] scatterplot of 100000 points and pdf file format

2004-11-24 Thread Thomas Lumley
On Wed, 24 Nov 2004, Witold Eryk Wolski wrote: Hi, I want to draw a scatter plot with 1M and more points and save it as pdf. Try the hexbin Bioconductor package, which gives hexagonally-binned density scatterplots. Even for tens of thousands of points this is often much better than a

Re: [R] scatterplot of 100000 points and pdf file format

2004-11-24 Thread Witold Eryk Wolski
Hi, Yes, indeed the hexbin package generates very cool pix. They look great. I was using it already. But this time I am interested in visualizing exactly the _scatter_ of some extreme points. Eryk Liaw, Andy wrote: Marc/Eryk, I have no experience with it, but I believe the hexbin package in

RE: [R] scatterplot of 100000 points and pdf file format

2004-11-24 Thread Prof Brian Ripley
On Wed, 24 Nov 2004 [EMAIL PROTECTED] wrote: On 24-Nov-04 Witold Eryk Wolski wrote: Hi, I want to draw a scatter plot with 1M and more points and save it as pdf. This makes the pdf file large. So i tried to save the file first as png and than convert it to pdf. This looks OK if printed but if

RE: [R] scatterplot of 100000 points and pdf file format

2004-11-24 Thread Matt Nelson
Witold, I have found that plotting more than a few thousand data points at a time quickly becomes a loosing proposition. That is, the dense overlap of data points tends to obscure the patterns of interest, with only outliers distinctly visible. I typically deal with this in two ways. The

[R] coplot =? gannt chart + bargraph

2004-11-24 Thread Jeff D. Hamann
I would like to display some results from simulations in the form of a Gantt chart (progress) with a barchart (production) of another variable below (something very similar to coplot charts). I'm not sure if I should attempt to build this from scratch (using grid or some of the basic graphics

Re: [R] OOT: frailty-multinivel

2004-11-24 Thread Thomas Lumley
On Wed, 24 Nov 2004, Kjetil Brinchmann Halvorsen wrote: Hola! I started to search for information about multilevel survival models, and found frailty in R. This seems to be something of the same, is it the same? More or less. Shared frailty models are the same as hierarchical/mixed survival

[R] what does order() stand for in an lme formula?

2004-11-24 Thread Harry Athanassiou
I'm a beginner in R, and trying to fit linear models with different intercepts per group, of the type y ~ A*x1 + B, where x1 is a numerical variable. I cannot understand whether I should use y1 ~ x1 +1 or y1 ~ order(x1) + 1 Although in the toy example included it makes a

Re: [R] scatterplot of 100000 points and pdf file format

2004-11-24 Thread John
On Wednesday 24 November 2004 07:34, Witold Eryk Wolski wrote: Hi, I want to draw a scatter plot with 1M and more points and save it as pdf. This makes the pdf file large. So i tried to save the file first as png and than convert it to pdf. This looks OK if printed but if viewed e.g. with

Re: [R] scatterplot of 100000 points and pdf file format

2004-11-24 Thread Sean Davis
Do you have a measures of scatter or can you pick outliers that could allow you to produce a mixed plot using either density or hexbinned data with only outliers placed after-the-fact using points()? Sean -Original Message- From: Witold Eryk Wolski [mailto:[EMAIL PROTECTED] Sent:

Re: [R] scatterplot of 100000 points and pdf file format

2004-11-24 Thread james . holtman
Have you tried plot(...,pch='.') This will use the period as the plotting character instead of the 'circle' which is drawn. This should reduce the size of the PDF file. I have done scatter plots with 2M points and they are typically meaningless with that many points overlaid. Check out

RE: [R] Re: T-test syntax question

2004-11-24 Thread Jessie F
Actually, you can still use t.test with one vector of data. Say, the differences is d ( a vector (or an array of numbers), you can use t.test(d), then by default, it testing whether mu=0, you can also specify confidence level by adding conf.level = 0.95 etc. You can also type ?t.test in R

Re: [R] scatterplot of 100000 points and pdf file format

2004-11-24 Thread Barry Rowlingson
I would strongly suggest a different method to present the data such as a contour plot or 3D bar plot. An XY plot with a million points is unlikely to be readable unless it is produced as a large format print. At 200 DPI printed, 1,000,000 discrete points requires a minimum of a 5 inch (12.7

[R] Searching for antilog function

2004-11-24 Thread Heather J. Branton
Dear R-users, I have a basic question about how to determine the antilog of a variable. Say I have some number, x, which is a factor of 2 such that x = 2^y. I want to figure out what y is, i.e. I am looking for the antilog base 2 of x. I have found log2 in the Reference Manual. But I am

Re: [R] what does order() stand for in an lme formula?

2004-11-24 Thread Peter Dalgaard
Harry Athanassiou [EMAIL PROTECTED] writes: I'm a beginner in R, and trying to fit linear models with different intercepts per group, of the type y ~ A*x1 + B, where x1 is a numerical variable. I cannot understand whether I should use y1 ~ x1 +1 or y1 ~ order(x1) + 1

RE: [R] SAS or R software

2004-11-24 Thread bogdan romocea
neela v writes: Hi all there Can some one clarify me on this issue, features wise which is better R or SAS, leaving the commerical aspect associated with it. I suppose there are few people who have worked on both R and SAS and wish they would be able to help me in deciding on this. THank

[R] SMVs

2004-11-24 Thread stephenc
Hi Everyone I am struggling to get going with support vector machines in R - smv() and predict() etc. Does anyone know of a good tutorial covering R and these things? Stephen [[alternative HTML version deleted]] __ [EMAIL PROTECTED]

RE: [R] Searching for antilog function

2004-11-24 Thread Liaw, Andy
What's wrong with log2()? log2(16) [1] 4 Isn't that exactly what you asked for? Andy From: Heather J. Branton Dear R-users, I have a basic question about how to determine the antilog of a variable. Say I have some number, x, which is a factor of 2 such that x = 2^y. I want to

Re: [R] Searching for antilog function

2004-11-24 Thread Spencer Graves
Consider: exp(log(1:11)) [1] 1 2 3 4 5 6 7 8 9 10 11 2^log(1:11, 2) [1] 1 2 3 4 5 6 7 8 9 10 11 2^logb(1:11, 2) [1] 1 2 3 4 5 6 7 8 9 10 11 10^log10(1:11) [1] 1 2 3 4 5 6 7 8 9 10 11 2^log2(1:11) [1] 1 2 3 4 5 6 7 8 9 10 11 Does this

Re: [R] Searching for antilog function

2004-11-24 Thread Duncan Murdoch
On Wed, 24 Nov 2004 12:26:46 -0500, Heather J. Branton [EMAIL PROTECTED] wrote : Dear R-users, I have a basic question about how to determine the antilog of a variable. Say I have some number, x, which is a factor of 2 such that x = 2^y. I want to figure out what y is, i.e. I am looking for

Re: [R] Searching for antilog function

2004-11-24 Thread Heather J. Branton
Yes! Somehow I must have made an entry error when I tried that before as I was getting something completely different! Thank you. ...heather Liaw, Andy wrote: What's wrong with log2()? log2(16) [1] 4 Isn't that exactly what you asked for? Andy From: Heather J. Branton Dear R-users, I

[R] seriesMerge

2004-11-24 Thread Yasser El-Zein
Is there a function in R that is equivalent to S-PLUS's seriesMerge(x1, x2, pos=union) where x1, and x2 are of class timeSeries seriesMerge is in S-PLUS's finmetrics. I looked into R's mergeSeries (in fSeries part of Rmetrics) but I could not make it behave quite the same. In R it expected a

Re: [R] scatterplot of 100000 points and pdf file format

2004-11-24 Thread Greg Snow
How about the following to plot only the 1,000 or so most extreem points (the outliers): x - rnorm(1e6) y - 2*x+rnorm(1e6) plot(x,y,pch='.') tmp - chull(x,y) while( length(tmp) 1000 ){ tmp - c(tmp, seq(along=x)[-tmp][ chull(x[-tmp],y[-tmp]) ] ) } points(x[tmp],y[tmp], col='red') now

Re: [R] Searching for antilog function

2004-11-24 Thread Heather J. Branton
Thank you so much for each of your responses. But to make sure I am clear (in my own mind), is this correct? If x = 2^y Then y = log2(x) Thanks again. I know this is basic. ...heather Duncan Murdoch wrote: On Wed, 24 Nov 2004 12:26:46 -0500, Heather J. Branton [EMAIL PROTECTED] wrote : Dear

Re: [R] reshaping of data for barplot2

2004-11-24 Thread Marc Schwartz
On Wed, 2004-11-24 at 19:24 +0100, Jean-Louis Abitbol wrote: Dear All, I have the following data coming out from s - with(final, summarize(norm, llist(gtt,fdiab), function(norm) { n - sum(!is.na(norm)) s -

Re: [R] seriesMerge

2004-11-24 Thread Dirk Eddelbuettel
On Wed, Nov 24, 2004 at 03:29:53PM -0500, Yasser El-Zein wrote: Is there a function in R that is equivalent to S-PLUS's seriesMerge(x1, x2, pos=union) where x1, and x2 are of class timeSeries seriesMerge is in S-PLUS's finmetrics. I looked into R's mergeSeries (in fSeries part of Rmetrics)

[R] confidence interval of a average...

2004-11-24 Thread Duncan Harris
I have a sample of lung capacities from a population measured against height. I need to know the 95% CI of the lung capacity of a person of average height. I have fitted a regression line. How do I get a minimum and maximum values of the 95% CI? My thinking was that this has something to do

[R] Installing gregmisc under windows 2000

2004-11-24 Thread Robert W. Baer, Ph.D.
I went through the following steps using RGUI menus to install gregmisc from CRAN. It appears to install but at the end R does not seem to be able to find it. Any idea what I'm doing wrong? Thankjs, Rob local({a - CRAN.packages() +

[R] Modeling censored binomial / Poisson data

2004-11-24 Thread Greg Tarpinian
I have some count data (0 - 1 at each time point for each test subject) that I would like to model. Since the 1's are rather sparse, the Poisson distribution comes to mind but I would also consider the binomial. The data are censored as the data come from a clinical trial and subjects were able

[R] how to remove time series trend in R?

2004-11-24 Thread Terry Mu
I got a set of data which has seasonal trend in form of sinx, cosx, I don't have any idea on how to deal with it. Can you give me a starting point? Thanks, Terry __ [EMAIL PROTECTED] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do

Re: [R] Installing gregmisc under windows 2000

2004-11-24 Thread Prof Brian Ripley
gregmisc is a bundle, not a package. Its description on CRAN is gregmiscBundle of gtools, gdata, gmodels, gplots so try one of those packages. On Wed, 24 Nov 2004, Robert W. Baer, Ph.D. wrote: I went through the following steps using RGUI menus to install gregmisc from CRAN. It appears

Re: [R] Installing gregmisc under windows 2000

2004-11-24 Thread Peter Dalgaard
Robert W. Baer, Ph.D. [EMAIL PROTECTED] writes: I went through the following steps using RGUI menus to install gregmisc from CRAN. It appears to install but at the end R does not seem to be able to find it. Any idea what I'm doing wrong? It's a bundle nowadays, so you need to load one of

[R] RE: RODBC and Table views

2004-11-24 Thread Niels Steen Krogh
channel2- odbcConnectAccess(C:\\Documents and Settings\\Fælles\\Journal\\DATASUPERMARKED\\DANBIONOVEMBER2004.mdb, uid=) sqlQuery(channel2,select * from Afdelinger_output_tabel1B order by antal desc) Does take views and tables! Niels Steen Krogh Konsulent ZiteLab Mail: -- [EMAIL

Re: [R] Installing gregmisc under windows 2000

2004-11-24 Thread Robert W. Baer, Ph.D.
Thanks for the clarification. Pursuant to the recent dicussion of GUI promoting ignornce among users, I plead guilty for CRAN installs but they have generally saved so much time.g. This does raise the question as to whether gregmisc and other bundles should appear on the install packages from

Re: [R] confidence interval of a average...

2004-11-24 Thread Robert W. Baer, Ph.D.
It depends on whether you want to do 95% ocnfidence intervals on the predicition or the mean vital capacity. Try the following and see if it gets you started: #Simulate data height=48:72 vc=height*10+20*rnorm(72-48+1) # Do regression lm.vc=lm(vc~height) # Confidence interval on mean vc

Re: [R] confidence interval of a average...

2004-11-24 Thread Robert W. Baer, Ph.D.
Sorry. The last code line got destroyed by my emailer and should read: matlines(height,predict.lm(lm.vc,interval=p), + lty=c(1,3,3),col=c('black','red','red')) - Original Message - From: Robert W. Baer, Ph.D. [EMAIL PROTECTED] To: [EMAIL PROTECTED] Sent: Wednesday, November 24, 2004

Re: [R] Installing gregmisc under windows 2000

2004-11-24 Thread Yuandan Zhang
That seems not the case under linux in term of installation. You can install this bundle in the same way as installing an individul package. eg R CMD INSTALL CRAN/contrib/main/gregmisc_2.0.0.tar.gz to get all constituent packages installed. On Wed, 24 Nov 2004 22:14:14 + (GMT) Prof

Re: [R] scatterplot of 100000 points and pdf file format

2004-11-24 Thread Patrick Connolly
On Wed, 24-Nov-2004 at 10:22AM -0600, Marc Schwartz wrote: | On Wed, 2004-11-24 at 16:34 +0100, Witold Eryk Wolski wrote: | Hi, | | I want to draw a scatter plot with 1M and more points and save it as pdf. | This makes the pdf file large. | So i tried to save the file first as png and than

RE: [R] Installing gregmisc under windows 2000

2004-11-24 Thread Liaw, Andy
If I'm not mistaken, bundle is really only useful as a concept for distribution and installation. You distribute and install a bundle, but load the individual packages when you want to use them. Once you install the bundle, you won't see the name of the bundle in the list of installed packages,

RE: [R] scatterplot of 100000 points and pdf file format

2004-11-24 Thread Ted Harding
On 24-Nov-04 Prof Brian Ripley wrote: On Wed, 24 Nov 2004 [EMAIL PROTECTED] wrote: 1. Multiply the data by some factor and then round the results to an integer (to avoid problems in step 2). Factor chosen so that the result of (4) below is satisfactory. 2. Eliminate duplicates in

RE: [R] scatterplot of 100000 points and pdf file format

2004-11-24 Thread Austin, Matt
-Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Behalf Of [EMAIL PROTECTED] Sent: Wednesday, November 24, 2004 16:37 PM To: R Help Mailing List Subject: RE: [R] scatterplot of 10 points and pdf file format On 24-Nov-04 Prof Brian Ripley wrote: On

RE: [R] scatterplot of 100000 points and pdf file format

2004-11-24 Thread Ted Harding
On 25-Nov-04 Ted Harding wrote: 'unique' will eat x for breakfast, indeed, but will have some trouble chewing (x,y). I still can't think of a neat way of doing that. Best wishes, Ted. Sorry, I don't want to be misunderstood. I didn't mean that 'unique' won't work for arrays. What I meant

Re: [R] confidence interval of a average...

2004-11-24 Thread Robert W. Baer, Ph.D.
Sorry if this was not clear. This is more of a theoreticla question rather than a R-coding question. I need to calculate The predicted response and 95% prediction interval for a man of average height So I need to predict the average response, which is easily done by taking the mean height

[R] Re:Hi!

2004-11-24 Thread APHA Registration Services
Thank you for your interest in APHA 2004. This is an automated reply confirming that we have received your email inquiry. General questions, requests, modifications and/or new registrations received by email, fax or mail will be processed within 7 business days. At that time, you will receive

RE: [R] scatterplot of 100000 points and pdf file format

2004-11-24 Thread Ted Harding
On 25-Nov-04 Austin, Matt wrote: -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Behalf Of [EMAIL PROTECTED] Sent: Wednesday, November 24, 2004 16:37 PM To: R Help Mailing List Subject: RE: [R] scatterplot of 10 points and pdf file format On 24-Nov-04

[R] Error in anova(): objects must inherit from classes

2004-11-24 Thread Andrew Criswell
Hello: Let me rephrase my question to attract interest in the problem I'm having. When I appply anova() to two equations estimated using glmmPQL, I get a complaint, anova(fm1, fm2) Error in anova.lme(fm1, fm2) : Objects must inherit from classes gls, gnls lm,lmList, lme,nlme,nlsList, or nls The

RE: [R] scatterplot of 100000 points and pdf file format

2004-11-24 Thread Liaw, Andy
From: [EMAIL PROTECTED] On 25-Nov-04 Ted Harding wrote: 'unique' will eat x for breakfast, indeed, but will have some trouble chewing (x,y). I still can't think of a neat way of doing that. Best wishes, Ted. Sorry, I don't want to be misunderstood. I didn't mean that

Re: [R] scatterplot of 100000 points and pdf file format

2004-11-24 Thread hadley wickham
Another possibility might be to use a 2d kernel density estimate (eg. kde2d from library(MASS). Then for the high density areas plot the density contours, for the low density areas plot the individual points. Hadley __ [EMAIL PROTECTED] mailing list

[R] logistic regression and 3PL model

2004-11-24 Thread Michael Lau
Hello colleagues, I am a novice with R and am stuck with an analysis I am trying to conduct. Any suggestions or feedback would be very much appreciated. I am analyzing a data set of psi (ESP) ganzfeld trials. The response variable is binary (correct/incorrect), with a 25% base rate. I've

RE: [R] Error in anova(): objects must inherit from classes

2004-11-24 Thread Austin, Matt
The lme method for anova() checks the inheritance of the object when a single object is supplied, which is why there is no error when you use one object at a time. When two objects are supplied, the method uses the class of the object by invoking the data.class function (which does not list

Re: [R] seriesMerge

2004-11-24 Thread Gabor Grothendieck
Yasser El-Zein abu3ammar at gmail.com writes: : : Is there a function in R that is equivalent to S-PLUS's : seriesMerge(x1, x2, pos=union) : where x1, and x2 are of class timeSeries : : seriesMerge is in S-PLUS's finmetrics. I looked into R's mergeSeries : (in fSeries part of Rmetrics) but I

Re: [R] Searching for antilog function

2004-11-24 Thread Gabor Grothendieck
Heather J. Branton hjb at pdq.com writes: : : Thank you so much for each of your responses. But to make sure I am : clear (in my own mind), is this correct? : : If x = 2^y : Then y = log2(x) : : Thanks again. I know this is basic. Although its not a proof, you can still use R to help you

Re: [R] Installing gregmisc under windows 2000

2004-11-24 Thread Prof Brian Ripley
On Thu, 25 Nov 2004, Yuandan Zhang wrote: That seems not the case under linux in term of installation. You can install this bundle in the same way as installing an individul package. eg R CMD INSTALL CRAN/contrib/main/gregmisc_2.0.0.tar.gz to get all constituent packages installed. And the same