[R] R: to view the memory

2007-09-11 Thread Sumit.Gupta
Hello, I am wondering if it is possible to view what variables and vairable values are stored in the R memory. This to enable debugging of R-scripts I write. Sumit [[alternative HTML version deleted]] __ R-help@r-project.org mailing list http

[R] Implementation of the butter() function for filtering.

2007-09-11 Thread Remund, Todd
How do I use the Butterworth filter from the signal library to simulate an antialias filter, that was analog? I know basically how to use the butter() function, but am not sure how to apply it to the data. Does the filter() function work. Can I use the plane='s', i.e. analog version? Thanks for

[R] Help regading time series data reading

2007-09-11 Thread Tolga Uzuner
Dear R-Users, Have a question about reading in some data and manipulating dates. I have a data set in excel which looks like this: Date PEGETNHSPPTC 13/10/200441.9264.7529.869.27 14/10/200441.9361.79

[R] Creating Hmisc labels inside a function

2007-09-11 Thread Steve Powell
hello list members, I am wanting to write a label using the Hmisc package, but inside a function. Normally I can do: library(Hmisc) M=2 label(M)="lab" #But inside a function the "=" does not work: Test=function(obj,labe) { label(obj)=labe } Test(M,"new label") I usually use the "assign" functio

[R] Create a "local" repository

2007-09-11 Thread Scott Hyde
I'd like to create a small "local" repository that would be used to install a package for a class of students at their home. I don't want to upload it to CRAN, as I don't think it should be disseminated at that level. What I'd like to do is: > where="http://mysite.com/"; > install.packages("mypa

Re: [R] Missing data

2007-09-11 Thread David Kaplan
Excellent. Thanks! === David Kaplan, Ph.D. Professor Department of Educational Psychology University of Wisconsin - Madison Educational Sciences, Room 1061 1025 W. Johnson Street Madison, WI 53706 email: [EMAIL PROTECTED] Web:

Re: [R] Inferences for ratios of non-normal means

2007-09-11 Thread Stephen Weigand
Paul, On 9/11/07, JRG <[EMAIL PROTECTED]> wrote: > On 11 Sep 2007 at 22:10, Robert A LaBudde wrote: > > > I think a ratio of two normals has a Cauchy distribution, which > > doesn't have a variance (the singularity in the denominator), so the > > Central Limit theorem does not apply. > > > > The C

Re: [R] Missing data

2007-09-11 Thread James Reilly
The mice package might also be useful, especially the md.pattern function: http://finzi.psych.upenn.edu/R/library/mice/html/md.pattern.html James -- James Reilly Department of Statistics, University of Auckland Private Bag 92019, Auckland, New Zealand On 12/9/07 1:33 PM, Bill Pikounis wrote: > D

[R] Integrate() error message, I am at a loss

2007-09-11 Thread Sergey Goriatchev
Hello! I have a problem with integrate() in my function nctspa(). Integrate produces an error message "evaluation of function gave a result of wrong length". I don't know what that means. Could anyone suggest me what is wrong with my function? These are the examples of function calls that work OK

[R] Problem with integrate()

2007-09-11 Thread Sergey Goriatchev
Hello! I have a problem with integrate() in my function nctspa(). Integrate produces an error message "evaluation of function gave a result of wrong length". I don't know what that means. Could anyone suggest me what is wrong with my function? These are the examples of function calls that work OK

[R] install R packages [SEC=UNCLASSIFIED]

2007-09-11 Thread Jin.Li
Hi All, I installed R 2.5.1 recently on a PC (Windows XP Professional 2001) and tried to install some R packages. It took several minutes and gave me the following message. > utils:::menuInstallPkgs() --- Please select a CRAN mirror for use in this session --- Error in open.connection(fil

Re: [R] Inferences for ratios of non-normal means

2007-09-11 Thread JRG
On 11 Sep 2007 at 22:10, Robert A LaBudde wrote: > I think a ratio of two normals has a Cauchy distribution, which > doesn't have a variance (the singularity in the denominator), so the > Central Limit theorem does not apply. > The Cauchy results if the denominator normal distribution has mean

Re: [R] Inferences for ratios of non-normal means

2007-09-11 Thread Moshe Olshansky
It even does not have a mean, but this does not matter: both the numerator and denominator are asymptotically normal and so the ratio converges (in distribution) to the ratio of two normals (which indeed has a Cauchy distribution). --- Robert A LaBudde <[EMAIL PROTECTED]> wrote: > I think a ratio

Re: [R] off-topic: better OS for statistical computing

2007-09-11 Thread Scionforbai
> My question is what is the best OS on PC (laptop) for statistical > computing and why. A free/open source *nix operating system is the best you can have todate, for almost everything, noticeably stability, security, scalability, networking and development, given that your hardware is supported a

Re: [R] Inferences for ratios of non-normal means

2007-09-11 Thread Robert A LaBudde
I think a ratio of two normals has a Cauchy distribution, which doesn't have a variance (the singularity in the denominator), so the Central Limit theorem does not apply. I would suggest using bootstrap resampling to make inferences. At 08:10 PM 9/11/2007, Moshe wrote: >For large samples you ha

Re: [R] Missing data

2007-09-11 Thread Bill Pikounis
David, Frank Harrell's pair of packages Hmisc and Design has some functions for tabulating, visualizing, and accounting for missing data. I recommend them as one avenue to investigate. Frank's companion book "Regression Modeling Strategies" covers their use in-depth. Hope that helps, Bill __

Re: [R] question about non-linear least squares in R

2007-09-11 Thread Yu (Warren) Wang
Dear Nilsson, Thank you very much for your reply which gave me much instruction! I will have a try on the suggested way from you. In fact, to get the value of MA is the most important thing for me in the fitting by this model. And I have thought about the linear model fitting at first

Re: [R] Are the error messages of ConstrOptim() consisten with each other?

2007-09-11 Thread Yuchen Luo
Dear Professor Lumley. You are absolutely right! Please allow me to rephrase: Assuming Professor Murdoch's explanation that theta=c(0.5,0.3,0.5) violate the constrain while theta1=c(0.5,0.9,0.5) satisfies it. I did the following # > theta=c(0.5,0.3,0.5) > theta1=c(0.5,0.9,

Re: [R] Inferences for ratios of non-normal means

2007-09-11 Thread Moshe Olshansky
For large samples you have asymptotic normality! --- Paul Smith <[EMAIL PROTECTED]> wrote: > Dear All, > > The package mratios can perform inferences for > ratios of normal means. > Is there some other package to do the same but with > non-normal > populations. Since I have got large samples, an

Re: [R] Missing data

2007-09-11 Thread Rolf Turner
On 12/09/2007, at 11:33 AM, David Kaplan wrote: > I did. If you don't want to answer, then your unhelpful suggestions to > yourself. This is uncalled for. Bert Gunter's response was perfectly reasonable. Read the posting guide. cheers,

Re: [R] Missing data

2007-09-11 Thread David Kaplan
I did. If you don't want to answer, then your unhelpful suggestions to yourself. === David Kaplan, Ph.D. Professor Department of Educational Psychology University of Wisconsin - Madison Educational Sciences, Room 1061 1025 W. Jo

Re: [R] Missing data

2007-09-11 Thread Bert Gunter
Please use R's existing search tools before posting: RsiteSearch("imputation") RsiteSearch("missing data imputation") etc. Bert Gunter Genentech Nonclinical Statistics -Original Message- From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] On Behalf Of David Kaplan Sent: Tuesday, Septem

[R] Missing data

2007-09-11 Thread David Kaplan
Hi all, I'm looking for a contributed package that can provide a detailed account of missing data patterns and perhaps also provide imputation procedures, such as mean imputation or hot deck imputation and the like. Is there anything out there? Thanks in advance, David -- ===

[R] Online course - Introduction to Data Handling in R.

2007-09-11 Thread Peter Bruce
Dr. Paul Murrell will present his course "Introduction to R - Data Handling" online at statistics.com Sept. 28 - Oct. 26. Participants can ask questions and exchange comments with Dr. Murrell via a private discussion board throughout the period. This course will provide an easy introduction to R

Re: [R] Vector Size Error Message

2007-09-11 Thread jim holtman
What operating system are you running on? Do you know if you have run out of physical memory? If on Windows, what do you have '--max-mem-size' set to on the command line invoking the RGUI? On 9/11/07, Sasse, Mark <[EMAIL PROTECTED]> wrote: > I have been using R off an on for approximately 3 mont

Re: [R] Fitting Data to a Noncentral Chi-Squared Distribution using MLE

2007-09-11 Thread Charles C. Berry
This is what I might do: > y <- rchisq( 1000, df=10, ncp=2 ) > library( stats4 ) > res <- mle( function(x,z) -sum( dchisq(y, x, z , log=TRUE ) ), start=list( > x=5, z=5 ) ) > coef(res) x z 10.355711 1.586123 > > ## or just to keep clear of boundary constraints: > > res <- mle(

Re: [R] building with atlas version of blas and lapack

2007-09-11 Thread Yuelin Li
On my linux machine (Ubuntu Feisty on i-686) this works: ./configure --with-blas="lf77blas -latlas" Yuelin. -- K Vanw wrote --|Tue (Sep/11/2007)[03:56]|--: I'd like to build R using my optimized blas and lapack libraries. It seems know matter what I do, the configure script uses the b

Re: [R] building with atlas version of blas and lapack

2007-09-11 Thread K Vanw
Looks like the list-server didn't like whatever gmail does to the formating. Sorry. So here's my post in plain text: On 9/11/07, K Vanw <[EMAIL PROTECTED]> wrote: > I'd like to build R using my optimized blas and lapack libraries. It seems > know matter what I > do, the configure script uses the

[R] building with atlas version of blas and lapack

2007-09-11 Thread K Vanw
I'd like to build R using my optimized blas and lapack libraries. It seems know matter what I do, the configure script uses the blas supplied with the source. My blas and lapack libraries are in /usr/local/atlas/lib. How can I get configure to use these? [[alternative HTML version deleted]

[R] Fitting Data to a Noncentral Chi-Squared Distribution using MLE

2007-09-11 Thread Terence Broderick
Hi, I have written out the log-likelihood function to fit some data I have (called ONES20) to the non-central chi-squared distribution. >library(stats4) >ll<-function(lambda,k){x<-ONES20; 25573*0.5*lambda-25573*log(2)-sum(-x/2)-log((x/lambda)^(0.25*k-0.5))-log(besselI(sqrt(lambda*x),0.5*

[R] Vector Size Error Message

2007-09-11 Thread Sasse, Mark
I have been using R off an on for approximately 3 months. As such, I am not very knowledgeable about coding in R. Currently I am running into the following problem. I'm using scan to make data available in R. The data I'm loading consists of over 400,000 records with 5 data fields. After runnin

[R] Fitting data to chi-squared or noncentral chi-squared distributions

2007-09-11 Thread Terence Broderick
Does anybody have any experience fitting data to the non-central chi-squared or chi-squared distribution? I am trying to fit some data to this distribution but there is error after error. audaces fortuna iuvat - [[alternative HTML version deleted]

Re: [R] xyplot question

2007-09-11 Thread Deepayan Sarkar
On 9/11/07, Leeds, Mark (IED) <[EMAIL PROTECTED]> wrote: > Thanks Deepyan. I am doing it in source but I didn't realize that that > was a function. Thanks so much. You might try source(..., echo=TRUE) if you want to mimic console evaluation. -Deepayan

Re: [R] Percentiles in R

2007-09-11 Thread Jose Luis Aznarte M.
Ok, of course you were right. As Stefan pointed out, it was a problem of the type of the vector: it was a factor and hence the quantile was not what expected. Thank you all!! Jose Luis Aznarte M. escribió: > Hi there! Still struggling to translate Matlab code into R's tsDyn package. > Here i

Re: [R] xyplot question

2007-09-11 Thread Leeds, Mark (IED)
Thanks Deepyan. I am doing it in source but I didn't realize that that was a function. Thanks so much. -Original Message- From: Deepayan Sarkar [mailto:[EMAIL PROTECTED] Sent: Tuesday, September 11, 2007 12:38 PM To: Leeds, Mark (IED) Cc: [EMAIL PROTECTED] Subject: Re: [R] xyplot questio

Re: [R] POSIXct dates on x-axis using xyplot

2007-09-11 Thread Deepayan Sarkar
On 9/10/07, jim holtman <[EMAIL PROTECTED]> wrote: > I am using 'xyplot' in lattice to plot some data where the x-axis is a > POSIXct date. I have data which spans a 6 month period, but when I > plot it, only the last month is printed on the right hand side of the > axis. I would have expected th

Re: [R] Percentiles in R

2007-09-11 Thread Daniel Lakeland
On Tue, Sep 11, 2007 at 07:09:24PM +0200, Jose Luis Aznarte M. wrote: > Thank you all! But the problem is that quantile() returns "how many > data" are greater than p percent, and not a value in the domain of the > vector under scrutiny. For example, I have a vector not quite, quantile retur

Re: [R] xyplot question

2007-09-11 Thread Deepayan Sarkar
On 9/11/07, Leeds, Mark (IED) <[EMAIL PROTECTED]> wrote: > I have the code below and it works fine if I print the xyplot but if I > take the print out, then I just get a blank > pdf. The same holds if I just send the plot to the console without the > print ( I get nothing ). My question is whether

Re: [R] Percentiles in R

2007-09-11 Thread roger koenker
Try quantile(x, .1, type=1) or read ?quantile for other options. url:www.econ.uiuc.edu/~rogerRoger Koenker email[EMAIL PROTECTED]Department of Economics vox: 217-333-4558University of Illinois fax: 217-244-6678

Re: [R] Percentiles in R

2007-09-11 Thread Stefan Grosse
Original Message Subject: Re: [R] Percentiles in R From: Jose Luis Aznarte M. <[EMAIL PROTECTED]> To: [EMAIL PROTECTED] Date: 11.09.2007 19:09 > Thank you all! But the problem is that quantile() returns "how many > data" are greater than p percent, and not a value in the dom

Re: [R] persp() problem

2007-09-11 Thread Petr PIKAL
Hi He is actually plotting DM matrix against default values [0,1]. So what he needs to do is to change y to c(1,10) and x to seq(1,20,2) to enable his code work. And even in that case there is not possible AFAIK to control exact labeling of axes as they are internaly recalculated to [0,1] value

Re: [R] Percentiles in R

2007-09-11 Thread Jose Luis Aznarte M.
Thank you all! But the problem is that quantile() returns "how many data" are greater than p percent, and not a value in the domain of the vector under scrutiny. For example, I have a vector > x [1] "-0,39559" "1,1916""0,23214" "1,656" "-0,3439" "-0,50568" [7] "0,52281" "0,

Re: [R] Odp: storing text and decimal values in a matrix form

2007-09-11 Thread Petr PIKAL
Hi Gonçalo Ferraz <[EMAIL PROTECTED]> napsal dne 11.09.2007 16:35:57: > Thanks! What is the easiest way of declaring an empty data frame. Is there > anything similar to: > matrix(0,x,y) > In the help files, I am only finding longer commands that ask me to specify > each individual column. Yo

Re: [R] xyplot question

2007-09-11 Thread Stefan Grosse
Original Message Subject: [R] xyplot question From: Leeds, Mark (IED) <[EMAIL PROTECTED]> To: [EMAIL PROTECTED] Date: 11.09.2007 18:00 > > load("stocks.dat") > > pdf('temp.pdf', width=10, height=8) > does ist work with postscript? something like: postscript("temp.ps", width =

[R] Meeting this week

2007-09-11 Thread Brent Sapiro
This email may be an advertisement or solicitation. If you do not want to receive marketing messages from WebEx click here | view Privacy Policy   -WEBEX:enUS:sdFj9CMZvQ8KWIBz21EK7I- __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listi

[R] CLUE package consensus function cl_medoid behavior

2007-09-11 Thread Joseph Retzer
When using cl_medoid to create a consensus clustering I've noticed that the list element in the resulting object containing the cluster assignments by case will switch between the [1] and [3] positions. This has happened when re-running the code without changes. Has anyone else experienced this

[R] xyplot question

2007-09-11 Thread Leeds, Mark (IED)
I have the code below and it works fine if I print the xyplot but if I take the print out, then I just get a blank pdf. The same holds if I just send the plot to the console without the print ( I get nothing ). My question is whether this is always the case with xyplot or is there something wrong w

Re: [R] Percentiles in R

2007-09-11 Thread Julian Burgos
Hola Jose Luis, Sure, you can use quantile(). Use the probs argument. For example, to get the median do quantile(X,probs=0.5). Do ?quantile to learn about the different type of quantiles calculated by the funcion. Saludos, Julian Jose Luis Aznarte M. wrote: > Hi there! Still struggling t

Re: [R] Percentiles in R

2007-09-11 Thread Xavier Abulker
It looks like prctile(X,p) is the same as quantile(X,p) i.e x<-0:100 quantile(x,0.5) is the median(x) and quantile(x,0.1)=10 is the value that is greater than 10% percent of the values in X "José Luis Aznarte M." wrote: > > Hi there! Still struggling to translate Matlab code into R's tsDyn >

Re: [R] using bootstrap for tree selection step in rpart

2007-09-11 Thread Frank E Harrell Jr
Fiona Callaghan wrote: > Hi Chris. Thanks for replying. I want to take the pruned list of > subtrees, and select the final tree from this list using a bootstrap > technique (rather than cross validation). I think, and I could be wrong, The result will be very poor predictive accuracy on new sa

Re: [R] persp() problem

2007-09-11 Thread Xavier Abulker
There is an error in your code: in persp(x , y , z) the length of x is the number of rows of z, the length of y is the number of columns of z, You should also name the rows and columns in z with x and y Economics Guy wrote: > > I am having some trouble getting the persp() package to change th

[R] Percentiles in R

2007-09-11 Thread Jose Luis Aznarte M.
Hi there! Still struggling to translate Matlab code into R's tsDyn package. Here is my question: Is there in R an equivalent function to Matlab's prctile()? To the moment I thought it was quantile(), but I just realized I was wrong. The definition of the Matlab function: prctile Percentiles of a

Re: [R] Editor for R under Fedora 7

2007-09-11 Thread Josué Polanco
Hi all, Another good alternative is VIM www.vim.org, vim Rules :D Too, check it in the historic list: http://tolstoy.newcastle.edu.au/R/help/06/05/26548.html cheerios ..-.. Josue On 9/11/07, Stefan Grosse <[EMAIL PROTECTED]> wrote: > I use emacs with the ess extension on Fedora 7. Emacs is

Re: [R] Editor for R under Fedora 7

2007-09-11 Thread Gavin Simpson
On Tue, 2007-09-11 at 08:36 -0500, [EMAIL PROTECTED] wrote: > Hello everyone, > > I wanted to ask what will be a good editor to write R scripts in Fedora 7. Emacs + ESS is my editor of choice. There is more info here: http://ess.r-project.org/ and a list of editors and the like can be found he

Re: [R] Editor for R under Fedora 7

2007-09-11 Thread Stefan Grosse
I use emacs with the ess extension on Fedora 7. Emacs is easily installable via yum and for ESS (emacs speaks statistics) you could use the Moertel rpms: http://community.moertel.com/ss/space/RPMs To edit: save a "buffer" as *.r, then some buttons appear which enable you to start an R process and

Re: [R] Editor for R under Fedora 7

2007-09-11 Thread Vladimir Eremeev
Did you consider R GUI project? http://www.sciviews.org/_rgui/ I used kate on Fedora, my general impression is not bad. But one should be careful with running R session inside kate. I killed it accidentally several times without saving of the state, when exited the editor. Antonio_Paredes wro

Re: [R] Editor for R under Fedora 7

2007-09-11 Thread Gustavo Carvalho
Fedora 7 is a linux distribution On 9/11/07, Paul Smith <[EMAIL PROTECTED]> wrote: > On 9/11/07, [EMAIL PROTECTED] > <[EMAIL PROTECTED]> wrote: > > I wanted to ask what will be a good editor to write R scripts in Fedora 7. > > Tell us first what is your operating system. > > Paul > > _

Re: [R] Editor for R under Fedora 7

2007-09-11 Thread Paul Smith
On 9/11/07, [EMAIL PROTECTED] <[EMAIL PROTECTED]> wrote: > > Fedora 7 What about kate then? Paul > "Paul Smith" <[EMAIL PROTECTED]> > Sent by: [EMAIL PROTECTED] > > 09/11/2007 09:25 AM > > To [EMAIL PROTECTED] > > cc > > > Subject Re: [R] Editor for R under Fedora 7 > > > > > > > > > On 9/11/07

Re: [R] Are the error messages of ConstrOptim() consisten with each other?

2007-09-11 Thread Thomas Lumley
On Mon, 10 Sep 2007, Yuchen Luo wrote: > Dear Professor Murdoch. > Thank you for your help! > 1. I believe c(0.5,0.3,0.5) satisfies the constrain because I did the > following experiment > ui=-1*ui > ci=-1*ci > constrOptim(c(0.5,0.3,0.5), f=fit.error, gr=fit.error.grr, ui=ui,ci=ci) > > The same er

Re: [R] Editor for R under Fedora 7

2007-09-11 Thread Paul Smith
On 9/11/07, [EMAIL PROTECTED] <[EMAIL PROTECTED]> wrote: > I wanted to ask what will be a good editor to write R scripts in Fedora 7. Tell us first what is your operating system. Paul __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/li

Re: [R] as.POSIXlt, starting at 0 or 18000?

2007-09-11 Thread Gabor Grothendieck
This is what I get: > unclass(strptime("1970-01-01", "%Y-%m-%d", tz="GMT") + 0) [1] 0 attr(,"tzone") [1] "GMT" > unclass(strptime("1970-01-01", "%Y-%m-%d", tz="EST5EDT") + 0) [1] 18000 attr(,"tzone") [1] "EST5EDT" > R.version.string # Windows Vista [1] "R version 2.6.0 alpha (2007-09-06 r42791)"

[R] as.POSIXlt, starting at 0 or 18000?

2007-09-11 Thread yoooooo
Hi all, I tried to do this: unclass(strptime("1970-01-01", "%Y-%m-%d", tz="") + 0) [1] 18000 I thought the starting time is zero? Same result returns if I do: unclass(strptime("1970-01-01", "%Y-%m-%d", tz="EDT") + 0) unclass(strptime("1970-01-01", "%Y-%m-%d", tz="EST") + 0) unclass(st

[R] Odp: storing text and decimal values in a matrix form

2007-09-11 Thread Petr PIKAL
Hi Although not specifically stated in help page, matrix is a vector with dim attribute. Therefore it can have only one type of values either numeric or character or logical or factor. So if you want textual and numeric values you probably have to stick with data frames. Regards Petr [EMAIL PR

[R] storing text and decimal values in a matrix form

2007-09-11 Thread Gonçalo Ferraz
Hi, I would like to store some text and decimal values in the following form: namedec.val dec.val ... dec.val namedec.val dec.val ... dec.val ... namedec. valdec.val ... dec.val To do so, I created a matrix of the necessary size (x row per y cols) with the com

Re: [R] using bootstrap for tree selection step in rpart

2007-09-11 Thread Fiona Callaghan
Hi Chris. Thanks for replying. I want to take the pruned list of subtrees, and select the final tree from this list using a bootstrap technique (rather than cross validation). I think, and I could be wrong, that bagging takes a bunch of bootstrap samples, grows one tree per boostrap sample and

[R] Editor for R under Fedora 7

2007-09-11 Thread Antonio_Paredes
Hello everyone, I wanted to ask what will be a good editor to write R scripts in Fedora 7. Tony [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

[R] why I cannot change the font and use hot-key in Editor of JGR 1.5-6

2007-09-11 Thread Xingwang Ye
Dear R users, Yesterday, I updated R 2.4.1 to R 2.5.1 on Unbuntu 7.0.4 successfully. I also update the JGR to 1.5-6 according to http://rosuda.org/JGR/linux.shtml (after "sudo update-java-alternatives -s java-6-sun", I cannot do "sudo update-java-alternatives -s java-1.6.0-sun"). I can use the

Re: [R] alternative way to loop

2007-09-11 Thread Duncan Murdoch
On 9/11/2007 8:54 AM, kevinchang wrote: > I heard that if-loop may result in errors if we are going to loop for a lot > of times. I don't know what you mean by an "if-loop", but if doing it lots of times caused an error, that would likely be a bug in R. And I know of no such bug. > And sum()

Re: [R] Generating Replicate Datasets (using loops or other means)

2007-09-11 Thread VTLT1999
Hello Moshe, Replication in this sense means that I want to use the parameters specified in a1 b1, and c1 to produce a set of simulated responses, coded as a 0 or a 1, and do this many times. So a1, b1, and c1 do not change, and I would like to produce 30 different sets of simulated responses, a

[R] alternative way to loop

2007-09-11 Thread kevinchang
I heard that if-loop may result in errors if we are going to loop for a lot of times. And sum() can be the alternative way to do this kind of hugbe looping. But I haven't figure out how . Can someone please give me an concrete example of using sum() for this purpose? -- View this message in con

[R] how to run a mixed design ANOVA with repeated measures

2007-09-11 Thread Gilbert G
Hi, I would like to run an ANOVA with repeated measures with both a between subjects variable and within subjects (repeated measures) variables. I know how to do this when I have equal numbers of subjects in both groups (i.e., balanced design). That has been very well documented by Baron, for exa

Re: [R] Function to get a sequence of months

2007-09-11 Thread Gabor Grothendieck
Here is a solution using character manipulation. Noting that month.name is built into R we paste together the character string: "January Februrary ... December January ... December" Then we use perl style ungreedy matching to get the shortest substring matching the indicated expression then splitt

Re: [R] Multiple contrast of ANOVA

2007-09-11 Thread Richard M. Heiberger
Download the multcomp package for multiple comparisons. You might also want to use the MMC functions in the HH package to display the contrasts. __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting

Re: [R] install packages automatically

2007-09-11 Thread john seers (IFR)
Hi Vladimir You may well be right and it is better, or at least better practice. But I have a feeling I tried it and it was not as successful, but I cannot remember why. Perhaps it was because when installing new packages they were not installed in that (USER) directory? JS -Original Mess

Re: [R] [solved] Re: I can't do it again on an other PC : R+RMySQL ->error loading dll

2007-09-11 Thread Peter Dalgaard
Ptit_Bleu wrote: > Hi, > > just to say that I "solved" my problem. I wrote "solved" because all I did > was to uninstall R and the DBI and RMySQL libraries and then install > everything again. > The only change is that I used another RMySQL.zip file found at this address > : > http://www.bioconduct

Re: [R] Function to get a sequence of months

2007-09-11 Thread Vladimir Eremeev
Duncan Murdoch-2 wrote: > > On 11/09/2007 6:36 AM, Arun Kumar Saha wrote: >> Hi all, >> >> I am looking for a function for following calculation. >> >> start.month = "July" >> end.month = "January" >> >> months = f(start.month, end.month, by=1) >> >> * f is the function that I am looking for

Re: [R] Function to get a sequence of months

2007-09-11 Thread Duncan Murdoch
On 11/09/2007 6:36 AM, Arun Kumar Saha wrote: > Hi all, > > I am looking for a function for following calculation. > > start.month = "July" > end.month = "January" > > months = f(start.month, end.month, by=1) > > * f is the function that I am looking for. > > Actually I want to get months = c(

Re: [R] off-topic: better OS for statistical computing

2007-09-11 Thread Patrick Connolly
On Mon, 10-Sep-2007 at 07:26PM -0400, Gabor Grothendieck wrote: |> My sense is that R users are even split between UNIX and Windows |> users so either will do in terms of the larger community. |> |> Some R packages may not be avaliable on every platform or will |> be available on one platform bef

Re: [R] [solved] Re: I can't do it again on an other PC : R+RMySQL ->error loading dll

2007-09-11 Thread Gustaf Rydevik
On 9/11/07, Peter Dalgaard <[EMAIL PROTECTED]> wrote: > Ptit_Bleu wrote: > > Hi Peter, > > > > I really have no idea about your comment. I was so happy to be able to > > connect to my database ... > > > > In fact I found the address in the following post : > > https://stat.ethz.ch/pipermail/r-help/

[R] Function to get a sequence of months

2007-09-11 Thread Arun Kumar Saha
Hi all, I am looking for a function for following calculation. start.month = "July" end.month = "January" months = f(start.month, end.month, by=1) * f is the function that I am looking for. Actually I want to get months = c("July", "August",.."January") If start.month = 6 and end.

Re: [R] question about non-linear least squares in R

2007-09-11 Thread Nilsson Fredrik X
Dear Warren, I had similar problems, this is (roughly) how I solved it translated to your problem: x <- c(-0.06,-0.04,-0.025,-0.015,-0.005,0.005,0.015,0.025,0.04,0.06) y <- c(1866760,1457870,1314960,1250560,1184850,1144920, 1158850,1199910,1263850,1452520) dafa<-data.frame(x,y) plot(x,y)

Re: [R] [solved] Re: I can't do it again on an other PC : R+RMySQL ->error loading dll

2007-09-11 Thread Peter Dalgaard
Ptit_Bleu wrote: > Hi Peter, > > I really have no idea about your comment. I was so happy to be able to > connect to my database ... > > In fact I found the address in the following post : > https://stat.ethz.ch/pipermail/r-help/2007-August/138142.html > > Do you think I can't use it freely ? :( >

[R] Inferences for ratios of non-normal means

2007-09-11 Thread Paul Smith
Dear All, The package mratios can perform inferences for ratios of normal means. Is there some other package to do the same but with non-normal populations. Since I have got large samples, an asymptotic procedure would be fine. Thanks in advance, Paul ___

Re: [R] [solved] Re: I can't do it again on an other PC : R+RMySQL ->error loading dll

2007-09-11 Thread Vladimir Eremeev
Do you need precisely RMySQL and DBI? RODBC can successfully replace this. I used RMySQL some time ago, but then, after the next version upgrade it stopped working. Switching to RODBC was almost without pain, since function names were quite similar. I simply replaced dbGetQuery with sqlGetQuery

Re: [R] [solved] Re: I can't do it again on an other PC : R+RMySQL ->error loading dll

2007-09-11 Thread Ptit_Bleu
Hi Peter, I really have no idea about your comment. I was so happy to be able to connect to my database ... In fact I found the address in the following post : https://stat.ethz.ch/pipermail/r-help/2007-August/138142.html Do you think I can't use it freely ? :( Thanks for your comments, Ptit

Re: [R] install packages automatically

2007-09-11 Thread Vladimir Eremeev
I think using R_LIBS_USER is better. >From ?.libPaths: The library search path is initialized at startup from the environment variable R_LIBS (which should be a semicolon-separated list of directories at which R library trees are rooted) followed by those in environment variable R_LIBS_USER.

[R] Multiple contrast of ANOVA

2007-09-11 Thread Tony zeng
Hi,all guy, I am doing ANOVA of some data with one factor,and the result shows factor's effect is not significant,but through some plots ,I know maybe some levels of the factor is significant,so I need to do multiple contrast to tell me which levels ,What I can do? thanks a lot!

Re: [R] install packages automatically

2007-09-11 Thread john seers (IFR)
Hi I use the following method which works well for me. I create a directory to hold the installed packages e.g. C:\Program Files\R\mylibrary Then in my C:\Program Files\R\R-2.5.1\.Renviron file I have a line as follows: R_LIBS=C:/PROGRA~1/R/mylibrary (I think you have to use Windows short

Re: [R] finding the minimum positive value of some data

2007-09-11 Thread Vladimir Eremeev
min(x[x>0]) dxc13 wrote: > > useRs, > > I am looking to find the minimum positive value of some data I have. > Currently, I am able to find the minimum of data after I apply some other > functions to it: > >> x > [1] 1 0 1 2 3 3 4 5 5 5 6 7 8 8 9 9 10 10 > >> sort(x) > [1]

[R] [solved] Re: I can't do it again on an other PC : R+RMySQL ->error loading dll

2007-09-11 Thread Ptit_Bleu
Hi, just to say that I "solved" my problem. I wrote "solved" because all I did was to uninstall R and the DBI and RMySQL libraries and then install everything again. The only change is that I used another RMySQL.zip file found at this address : http://www.bioconductor.org/packages/2.0/extra/bin/w

[R] composition of a matrix by matches between two data frames

2007-09-11 Thread Patrick Zimmermann
Dear R-users, I want to compose a new table (a presence/absence matrix of species in regions) by matches between two data frames. The first data frame has the species as header and the habitats in which they occure as elements. The second data frame has the regions as header and the occuring habit

Re: [R] install packages automatically

2007-09-11 Thread Vladimir Eremeev
This was discussed here sometime ago. Several variants are possible. You can install the new version of R over the old one, in the same folder, this preserves the installed packages. At least, this works fine for me on Windows. You can also store your list of packages in a simple R script and ex