Re: [R] iplots problem

2007-05-25 Thread ryestone
You have to use Sun's java on your PC and not microsoft's. I just loaded it on mine and then restarted computer and it seemed to work. Hope this helps! mister_bluesman wrote: > > Hi > > How did u 'load' sun's java in R? > > Many thanks > > > > ryestone wrote: >> >> Did you load Sun's java

Re: [R] Estimation of Dispersion parameter in GLM for Gamma Dist.

2007-05-25 Thread Prof Brian Ripley
This is discussed in the book the MASS package (sic) supports, and/or its online material (depending on the edition). On Fri, 25 May 2007, fredrik odegaard wrote: > Hi All, > could someone shed some light on what the difference between the > estimated dispersion parameter that is supplied with t

Re: [R] Problem with rpart

2007-05-25 Thread Prof Brian Ripley
You only have 43 cases. After one split, the groups are too small to split again with the default settings. See ?rpart.control. On Fri, 25 May 2007, Silvia Lomascolo wrote: > > I work on Windows, R version 2.4.1. I'm very new with R! > > I am trying to build a classification tree using rpa

Re: [R] How to get the "Naive SE" of coefficients from the zelig output

2007-05-25 Thread Ferdinand Alimadhi
Dear Abdus, Can you try: summary(il6w.out)$table[,"(Naive SE)"] If you want to know more details of where that's coming from, run the following command in your R prompt and have a look at the outputted code > survival:::summary.survreg You might also consider subscribing at the zelig mailing

Re: [R] Interactive plots?

2007-05-25 Thread Michael Lawrence
You could load it into GGobi using the rggobi package and then use the identification mode, as well as all the other interactive features... http://www.ggobi.org Michael On 5/25/07, mister_bluesman <[EMAIL PROTECTED]> wrote: > > > Hi there. > > I have a matrix that provides place names and the d

Re: [R] How to get the "Naive SE" of coefficients from the zelig output

2007-05-25 Thread Abdus Sattar
Dear R-user: After the fitting the Tobit model using zelig, if I use the following command then I can get the regression coefficents: beta=coefficients(il6.out) > beta (Intercept) apache 4.7826 0.9655 How may I extract the "Naive SE" from the following output please? > summary

Re: [R] 3D plots with data.frame

2007-05-25 Thread Duncan Murdoch
On 25/05/2007 8:22 PM, [EMAIL PROTECTED] wrote: > > You could try the function 'plot3d', in package 'rgl': > > library(rgl) > ?plot3d > x<-data.frame(a=rnorm(100),b=rnorm(100),c=rnorm(100)) > plot3d(x$a,x$b,x$c) Or more simply, plot3d(x) (which plots the 1st three columns). Duncan Murdoch > >

Re: [R] 3D plots with data.frame

2007-05-25 Thread J . delasHeras
You could try the function 'plot3d', in package 'rgl': library(rgl) ?plot3d x<-data.frame(a=rnorm(100),b=rnorm(100),c=rnorm(100)) plot3d(x$a,x$b,x$c) Jose Quoting "H. Paul Benton" <[EMAIL PROTECTED]>: > Dear all, > > Thank you for any help. I have a data.frame and would like to plot > it

Re: [R] normality tests [Broadcast]

2007-05-25 Thread Frank E Harrell Jr
[EMAIL PROTECTED] wrote: > Following up on Frank's thought, why is it that parametric tests are so > much more popular than their non-parametric counterparts? As > non-parametric tests require fewer assumptions, why aren't they the > default? The relative efficiency of the Wilcoxon test as compar

Re: [R] normality tests [Broadcast]

2007-05-25 Thread Cody_Hamilton
Following up on Frank's thought, why is it that parametric tests are so much more popular than their non-parametric counterparts? As non-parametric tests require fewer assumptions, why aren't they the default? The relative efficiency of the Wilcoxon test as compared to the t-test is 0.955, and y

[R] Estimation of Dispersion parameter in GLM for Gamma Dist.

2007-05-25 Thread fredrik odegaard
Hi All, could someone shed some light on what the difference between the estimated dispersion parameter that is supplied with the GLM function and the one that the 'gamma.dispersion( )' function in the MASS library gives? And is there consensus for which estimated value to use? It seems that the

Re: [R] normality tests [Broadcast]

2007-05-25 Thread Cody_Hamilton
You can also try validating your regression model via the bootstrap (the validate() function in the Design library is very helpful). To my mind that would be much more reassuring than normality tests performed on twenty residuals. By the way, be careful with the correlation test - it's only good

Re: [R] Interactive plots?

2007-05-25 Thread Tony Plate
The package RSVGTipsDevice allows you to do just it just -- you create a plot in an SVG file that can be viewed in a browser like FireFox, and the points (or shapes) in that plot can have pop-up tooltips. -- Tony Plate mister_bluesman wrote: > Hi there. > > I have a matrix that provides place

Re: [R] normality tests [Broadcast]

2007-05-25 Thread wssecn
The normality of the residuals is important in the inference procedures for the classical linear regression model, and normality is very important in correlation analysis (second moment)... Washington S. Silva > Thank you all for your replies they have been more useful... well > in my case

[R] Interactive plots?

2007-05-25 Thread mister_bluesman
Hi there. I have a matrix that provides place names and the distances between them: Chelt Exeter London Birm Chelt 0 118 96 50 Exeter 1180 118 163 London 96 118 0 118 Birm 50 163 118 0 After performing multi

Re: [R] normality tests [Broadcast]

2007-05-25 Thread Frank E Harrell Jr
[EMAIL PROTECTED] wrote: > Thank you all for your replies they have been more useful... well > in my case I have chosen to do some parametric tests (more precisely > correlation and linear regressions among some variables)... so it > would be nice if I had an extra bit of support on my decision

Re: [R] normality tests [Broadcast]

2007-05-25 Thread Frank E Harrell Jr
Lucke, Joseph F wrote: > Most standard tests, such as t-tests and ANOVA, are fairly resistant to > non-normalilty for significance testing. It's the sample means that have > to be normal, not the data. The CLT kicks in fairly quickly. Testing > for normality prior to choosing a test statistic is

[R] 3D plots with data.frame

2007-05-25 Thread H. Paul Benton
Dear all, Thank you for any help. I have a data.frame and would like to plot it in 3D. I have tried wireframe() and cloud(), I got scatterplot3d(xs) Error: could not find function "scatterplot3d" > wireframe(xs) Error in wireframe(xs) : no applicable method for "wireframe" > persp(x=x, y=

[R] Problem with rpart

2007-05-25 Thread Silvia Lomascolo
I work on Windows, R version 2.4.1. I'm very new with R! I am trying to build a classification tree using rpart but, although the matrix has 108 variables, the program builds a tree with only one split using one variable! I know it is probable that only one variable is informative, but I think

Re: [R] iplots problem

2007-05-25 Thread mister_bluesman
Hi How did u 'load' sun's java in R? Many thanks ryestone wrote: > > Did you load Sun's java? Try that and also try rebooting your machine. > This worked for me. > > > mister_bluesman wrote: >> >> Hi. I try to load iplots using the following commands >> >>> library(rJava) >>> library(ipl

Re: [R] covariance question which has nothing to do with R

2007-05-25 Thread toby909
while my other program is running. The reference I mentioned previously addresses exactly this. Snijders and Bosker's Multilevel Analysis book on page 31 and 33, section 3.6.2 and 363 discuss this. When you say that the Xs are correlated then you would need to say according to which struct

[R] Scale mixture of normals

2007-05-25 Thread Anup Nandialath
Dear Friends, Is there an R package which implements regression models with error distributions following a scale mixture of normals? Thanks Anup __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the

Re: [R] email the silly fuckers instead.

2007-05-25 Thread mister_bluesman
What? I'm not sure about anyone else but I have absolutely no idea what you're going on about. foxy-3 wrote: > > An alle finanzinvestoren! Diese aktie wird durchstarten! Dienstag 29. > Mai startet die hausse! > > Firma: Talktech Telemedia > Kurzel: OYQ / OYQ.F / OYQ.DE > WKN: 278104 > ISIN: US

[R] email the silly fuckers instead.

2007-05-25 Thread foxy
An alle finanzinvestoren! Diese aktie wird durchstarten! Dienstag 29. Mai startet die hausse! Firma: Talktech Telemedia Kurzel: OYQ / OYQ.F / OYQ.DE WKN: 278104 ISIN: US8742781045 Preis: 0.54 5-T Prognose: 1.75 Realisierter kursgweinn von 450% in 5 tagen! OYQ Wird wie eine rakete durchstarten! D

Re: [R] File path expansion

2007-05-25 Thread Duncan Murdoch
On 5/25/2007 1:09 PM, Prof Brian Ripley wrote: > On Fri, 25 May 2007, Martin Maechler wrote: > >> >>> path.expand("~") >> [1] "/home/maechler" > > Yes, but beware that may not do what you want on Windows in R <= 2.5.0, > since someone changed the definition of 'home' but not path.expand. A more

Re: [R] In which package is the errbar command located?

2007-05-25 Thread Charles C. Berry
help.search("errbar") shows the one installed on my system (and might on yours), but RSiteSearch("errbar") shows that more than one package contains a function by that name On Fri, 25 May 2007, Judith Flores wrote: > > __ > R-help@

Re: [R] Calculation of ratio distribution properties

2007-05-25 Thread Ravi Varadhan
Mike, Attached is an R function to do this, along with an example that will reproduce the MathCad plot shown in your attached paper. I haven't checked it thoroughly, but it seems to reproduce the MathCad example well. Ravi.

Re: [R] lme with corAR1 errors - can't find AR coefficient in output

2007-05-25 Thread Stephen Weigand
Millo, On 5/24/07, Millo Giovanni <[EMAIL PROTECTED]> wrote: > Dear List, > > I am using the output of a ML estimation on a random effects model with > first-order autocorrelation to make a further conditional test. My model > is much like this (which reproduces the method on the famous Grunfeld

Re: [R] iplots problem

2007-05-25 Thread ryestone
Did you load Sun's java? Try that and also try rebooting your machine. This worked for me. mister_bluesman wrote: > > Hi. I try to load iplots using the following commands > >> library(rJava) >> library(iplots) > > but then I get the following error: > > Error in .jinit(cp, parameters = "-Xm

[R] /tmp/ gets filled up fast

2007-05-25 Thread Alessandro Gagliardi
Dear useRs, I'm running some pretty big R scripts using a PBS that calls upon the RMySQL library and it's filling up the /tmp/ directory with Rtmp* files. I thought the problem might have come from scripts crashing and therefore not getting around to executing the dbDisconnect() function, but I j

Re: [R] normality tests [Broadcast]

2007-05-25 Thread gatemaze
Thank you all for your replies they have been more useful... well in my case I have chosen to do some parametric tests (more precisely correlation and linear regressions among some variables)... so it would be nice if I had an extra bit of support on my decisions... If I understood well from al

[R] iPlots package

2007-05-25 Thread ryestone
I am having trouble connecting two points in iplot. In the normal plot command I would use segments(). I know there is a function ilines() but can you just enter coordinates of 2 points? -- View this message in context: http://www.nabble.com/iPlots-package-tf3817683.html#a10808180 Sent from the

Re: [R] In which package is the errbar command located?

2007-05-25 Thread John Kane
Hmisc --- Judith Flores <[EMAIL PROTECTED]> wrote: __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-con

Re: [R] Calculation of ratio distribution properties

2007-05-25 Thread Mike Lawrence
According to the paper I cited, there is controversy over the sufficiency of Hinkley's solution, hence their proposed more complete solution. On 25-May-07, at 2:45 PM, Lucke, Joseph F wrote: > The exact ratio is given in > > On the Ratio of Two Correlated Normal Random Variables, D. V. > Hi

Re: [R] xyplot: different scales accross rows, same scales within rows

2007-05-25 Thread Deepayan Sarkar
On 5/25/07, Marta Rufino <[EMAIL PROTECTED]> wrote: > Dear list members, > > > I would like to set up a multiple panel in xyplots, with the same scale > for all colunms in each row, but different accross rows. > relation="free" would set up all x or y scales free... which is not what > I want :-( >

Re: [R] xyplot: different scales accross rows, same scales within rows

2007-05-25 Thread Gabor Grothendieck
xlim= can take a list: # CO2 is built into R library(lattice) xlim <- rep(list(c(0, 1000), c(0, 2000)), each = 2) xyplot(uptake ~ conc | Type * Treatment, data = CO2, scales = list(relation = "free"), xlim = xlim) On 5/25/07, Marta Rufino <[EMAIL PROTECTED]> wrote: > Dear list members, >

Re: [R] normality tests [Broadcast]

2007-05-25 Thread Lucke, Joseph F
Most standard tests, such as t-tests and ANOVA, are fairly resistant to non-normalilty for significance testing. It's the sample means that have to be normal, not the data. The CLT kicks in fairly quickly. Testing for normality prior to choosing a test statistic is generally not a good idea. -

[R] In which package is the errbar command located?

2007-05-25 Thread Judith Flores
__ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

Re: [R] Calculation of ratio distribution properties

2007-05-25 Thread Mike Lawrence
MathCad code failed to attach last time. Here it is. On 25-May-07, at 2:24 PM, Mike Lawrence wrote: I came across this reference: http://www.informaworld.com/smpp/content? content=10.1080/03610920600683689 The authors sent me code (attached with permission) in MathCad to perform the calcu

Re: [R] File path expansion

2007-05-25 Thread Prof Brian Ripley
On Fri, 25 May 2007, Martin Maechler wrote: > >> path.expand("~") > [1] "/home/maechler" Yes, but beware that may not do what you want on Windows in R <= 2.5.0, since someone changed the definition of 'home' but not path.expand. > >> "RobMcG" == McGehee, Robert <[EMAIL PROTECTED]> >>

Re: [R] normality tests [Broadcast]

2007-05-25 Thread Liaw, Andy
From: [EMAIL PROTECTED] > > On 25/05/07, Frank E Harrell Jr <[EMAIL PROTECTED]> wrote: > > [EMAIL PROTECTED] wrote: > > > Hi all, > > > > > > apologies for seeking advice on a general stats question. I ve run > > > normality tests using 8 different methods: > > > - Lilliefors > > > - Shapiro-Wilk

[R] xyplot: different scales accross rows, same scales within rows

2007-05-25 Thread Marta Rufino
Dear list members, I would like to set up a multiple panel in xyplots, with the same scale for all colunms in each row, but different accross rows. relation="free" would set up all x or y scales free... which is not what I want :-( Is this possible? Thank you in advance, Best wishes, Marta

Re: [R] windows to unix

2007-05-25 Thread Barry Rowlingson
Martin Maechler wrote: >> "Erin" == Erin Hodgess <[EMAIL PROTECTED]> >> on Fri, 25 May 2007 06:10:10 -0500 writes: > > Erin> Dear R People: > Erin> Is there any way to take a Windows version of R, compiled from > source, > Erin> compress it, and put it on a Unix-like envi

Re: [R] normality tests

2007-05-25 Thread gatemaze
On 25/05/07, Frank E Harrell Jr <[EMAIL PROTECTED]> wrote: > [EMAIL PROTECTED] wrote: > > Hi all, > > > > apologies for seeking advice on a general stats question. I ve run > > normality tests using 8 different methods: > > - Lilliefors > > - Shapiro-Wilk > > - Robust Jarque Bera > > - Jarque Bera

Re: [R] File path expansion

2007-05-25 Thread Martin Maechler
> path.expand("~") [1] "/home/maechler" > "RobMcG" == McGehee, Robert <[EMAIL PROTECTED]> > on Fri, 25 May 2007 11:44:27 -0400 writes: RobMcG> R-Help, RobMcG> I discovered a "mis-feature" is ghostscript, which is used by the bitmap RobMcG> function. It seems that specify

Re: [R] normality tests

2007-05-25 Thread Frank E Harrell Jr
[EMAIL PROTECTED] wrote: > Hi all, > > apologies for seeking advice on a general stats question. I ve run > normality tests using 8 different methods: > - Lilliefors > - Shapiro-Wilk > - Robust Jarque Bera > - Jarque Bera > - Anderson-Darling > - Pearson chi-square > - Cramer-von Mises > - Shapiro

Re: [R] windows to unix

2007-05-25 Thread Martin Maechler
> "Erin" == Erin Hodgess <[EMAIL PROTECTED]> > on Fri, 25 May 2007 06:10:10 -0500 writes: Erin> Dear R People: Erin> Is there any way to take a Windows version of R, compiled from source, Erin> compress it, and put it on a Unix-like environment, please? Since nobody has

Re: [R] trouble with snow and Rmpi

2007-05-25 Thread Ramon Diaz-Uriarte
Dear Erin, What operating system are you trying this on? Windows? In Linux you definitely don't need MPICH2 but, rather, LAM/MPI. Best, R. On 5/25/07, Erin Hodgess <[EMAIL PROTECTED]> wrote: > Dear R People: > > I am having some trouble with the snow package. > > It requires MPICH2 and Rmpi. >

Re: [R] File path expansion

2007-05-25 Thread Gabor Grothendieck
Try ?path.expand On 5/25/07, McGehee, Robert <[EMAIL PROTECTED]> wrote: > R-Help, > I discovered a "mis-feature" is ghostscript, which is used by the bitmap > function. It seems that specifying file names in the form "~/abc.png" > rather than "/home/directory/abc.png" causes my GS to crash when I

[R] File path expansion

2007-05-25 Thread McGehee, Robert
R-Help, I discovered a "mis-feature" is ghostscript, which is used by the bitmap function. It seems that specifying file names in the form "~/abc.png" rather than "/home/directory/abc.png" causes my GS to crash when I open the bitmap device on my Linux box. The easiest solution would seem to be to

Re: [R] Read in 250K snp chips

2007-05-25 Thread James W. MacDonald
bhiggs wrote: > I'm having trouble getting summaries out of the 250K snp chips in R. I'm > using the oligo package and when I attempt to create the necessary SnpQSet > object (to get genotype calls and intensities) using snprma, I encounter > memory issues. > > Anyone have an alternative package

[R] Cochran-Armitage

2007-05-25 Thread [EMAIL PROTECTED]
Hi, I have to calculate a Cochran-Armitage test on data. I don't know much on this test and then, I have few questions on the independence_test function. I have read the help and search the archives but this doesn't help me... - They are two ways: either with a formula (independence_test(tum

[R] Read in 250K snp chips

2007-05-25 Thread bhiggs
I'm having trouble getting summaries out of the 250K snp chips in R. I'm using the oligo package and when I attempt to create the necessary SnpQSet object (to get genotype calls and intensities) using snprma, I encounter memory issues. Anyone have an alternative package or workaround for these l

[R] Help with complex lme model fit

2007-05-25 Thread Colin Beale
Hi R helpers, I'm trying to fit a rather complex model to some simulated data using lme and am not getting the correct results. It seems there might be some identifiability issues that could possibly be dealt with by specifying starting parameters - but I can't see how to do this. I'm compari

Re: [R] Make check failure for R-2.4.1

2007-05-25 Thread Adam Witney
> Adam> Here is the results: > >>> xB <- c(2000,1e6,1e50,Inf) >>> for(df in c(0.1, 1, 10)) > Adam> +for(ncp in c(0, 1, 10, 100)) > Adam> +print(pchisq(xB, df=df, ncp=ncp), digits == 15) > Adam> Error in print.default(pchisq(xB, df = df, ncp = ncp), digits ==

Re: [R] Make check failure for R-2.4.1

2007-05-25 Thread Martin Maechler
> "Adam" == Adam Witney <[EMAIL PROTECTED]> > on Fri, 25 May 2007 14:48:18 +0100 writes: >>> ##-- non central Chi^2 : >>> xB <- c(2000,1e6,1e50,Inf) >>> for(df in c(0.1, 1, 10)) >> + for(ncp in c(0, 1, 10, 100)) stopifnot(pchisq(xB, df=df, ncp=ncp) ==1) >> Erro

Re: [R] Problem with numerical integration and optimization with BFGS

2007-05-25 Thread Deepankar Basu
Ravi, Thanks a lot for your detailed suggestions. I will certainly look at the links that you have sent and the package "mnormt". For the moment, I have managed to analytically integrate the expression using "pnorm" along the lines suggested by Prof. Ripley yesterday. For instance, my first inte

Re: [R] Problem with numerical integration and optimization with BFGS

2007-05-25 Thread Ravi Varadhan
Deepankar, If the problem seems to be in the evaluation of numerical quadrature part, you might want to try quadrature methods that are better suited to integrands with strong peaks. The traditional Gaussian quadrature methods, even their adaptive versions such as Gauss-Kronrod, are not best suit

Re: [R] Make check failure for R-2.4.1

2007-05-25 Thread Adam Witney
>> ##-- non central Chi^2 : >> xB <- c(2000,1e6,1e50,Inf) >> for(df in c(0.1, 1, 10)) > + for(ncp in c(0, 1, 10, 100)) stopifnot(pchisq(xB, df=df, ncp=ncp) ==1) > Error: pchisq(xB, df = df, ncp = ncp) == 1 is not all TRUE > Execution halted > > Ok, thanks; > so, if we want to learn mo

[R] normality tests

2007-05-25 Thread gatemaze
Hi all, apologies for seeking advice on a general stats question. I ve run normality tests using 8 different methods: - Lilliefors - Shapiro-Wilk - Robust Jarque Bera - Jarque Bera - Anderson-Darling - Pearson chi-square - Cramer-von Mises - Shapiro-Francia All show that the null hypothesis that

[R] Competing Risks Analysis

2007-05-25 Thread Kevin E. Thorpe
I am working on a competing risks problem, specifically an analysis of cause-specific mortality. I am familiar with the cmprsk package and have used it before to create cumulative incidence plots. I also came across an old (1998) s-news post from Dr. Terry Therneau describing a way to use coxph t

Re: [R] Make check failure for R-2.4.1

2007-05-25 Thread Martin Maechler
> "Adam" == Adam Witney <[EMAIL PROTECTED]> > on Fri, 25 May 2007 09:38:29 +0100 writes: Adam> Thanks for your replies Details inline below: Adam> On 24/5/07 17:12, "Martin Maechler" <[EMAIL PROTECTED]> wrote: >>> "UweL" == Uwe Ligges <[EMAIL PROTECTED]> >

Re: [R] R-About PLSR

2007-05-25 Thread Gavin Simpson
On Fri, 2007-05-25 at 17:25 +0530, Nitish Kumar Mishra wrote: > hi R help group, > I have installed PLS package in R and use it for princomp & prcomp > commands for calculating PCA using its example file(USArrests example). > But How I can use PLS for Partial least square, R square, mvrCv one more

[R] R-About PLSR

2007-05-25 Thread Nitish Kumar Mishra
hi R help group, I have installed PLS package in R and use it for princomp & prcomp commands for calculating PCA using its example file(USArrests example). But How I can use PLS for Partial least square, R square, mvrCv one more think how i can import external file in R. When I use plsr, R2, RMSEP

[R] lmer and scale parameter in glmm model

2007-05-25 Thread Olivier MARTIN
Hi all, I try to fit a glmm model with binomial distribution and I would to verify that the scale parameter is close to 1... the lmer function gives the following result : Estimated scale (compare to 1 ) 0.766783 But I would like to know how this estimation (0.766783) is performed, and I woul

[R] iplots problem

2007-05-25 Thread mister_bluesman
Hi. I try to load iplots using the following commands > library(rJava) > library(iplots) but then I get the following error: Error in .jinit(cp, parameters = "-Xmx512m", silent = TRUE) : Cannot create Java Virtual Machine Error in library(iplots) : .First.lib failed for 'iplots' What

[R] windows to unix

2007-05-25 Thread Erin Hodgess
Dear R People: Is there any way to take a Windows version of R, compiled from source, compress it, and put it on a Unix-like environment, please? thanks in advance, Sincerely, Erin Hodgess Associate Professor Department of Computer and Mathematical Sciences University of Houston - Downtown mailt

[R] trouble with snow and Rmpi

2007-05-25 Thread Erin Hodgess
Dear R People: I am having some trouble with the snow package. It requires MPICH2 and Rmpi. Rmpi is fine. However, I downloaded the MPICH2 package, and installed. There is no mpicc, mpirun, etc. Does anyone have any suggestions, please? Thanks in advance! Sincerely, Erin Hodgess Associate P

Re: [R] in unix opening data object created under win

2007-05-25 Thread John Kane
--- [EMAIL PROTECTED] wrote: > On Unix R's version is 2.3.1 and on PC its 2.4.1. > > I dont have the rights to install newer version of R > on Unix. > > I tried different upload methods. No one worked. > > On Unix it looks as follows (dots to hide my > userid): > > > > load("/afs/ir/users/..

Re: [R] how to mimic plot=F for truehist?

2007-05-25 Thread Vladimir Eremeev
By defining your own function. You can get the function body by typing its name in the R command line and pressing Enter. Copy-paste the function body in ascii file (source R code), redefine it as you like, for example, by adding desired argument and code for processing it, then source that file

Re: [R] L1-SVM

2007-05-25 Thread Gorden T Jemwa
> Is there a package in R to find L1-SVM. I did search and found svmpath > but not sure if it is what I need. Try kernlab __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www

[R] testing difference (or similarities) between two distance matrices (not independent)

2007-05-25 Thread Stephane . Buhler
Hi, i'm looking to test if two distance matrices are statistically different from each others. These two matrices have been computed on the same set of data (several population samples) 1. using a particular genetic distance 2. weighting that genetic distance with an extra factor (we can lo

[R] how to mimic plot=F for truehist?

2007-05-25 Thread Johannes Graumann
Dear Rologists, In order to combine plots I need to get access to the some "par"s specific to my plot prior to replot it with modified parameters. I have not found any option like "plot=F" associated with truehist and would like to know whether someone can point out how to overcome this problem.

Re: [R] Speeding up resampling of rows from a large matrix

2007-05-25 Thread Juan Pablo Lewinger
That's beautiful. For the full 120 x 65,000 matrix your approach took 85 seconds. A truly remarkable improvement over my 80 minutes! Thank you! __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the pos

Re: [R] Question concerning "pastecs" package

2007-05-25 Thread Rainer M. Krug
Dear Philippe Thanks a lot for your information - it is extremely usefull Rainer Philippe Grosjean wrote: > Hello, > > I already answered privately to your question. No, there is no > translation of pastecs.pdf. The English documentation is accessible, as > usual, by: > > ?turnpoints > > R

[R] Off topic: S.E. for cross validation

2007-05-25 Thread Gad Abraham
Hi, I'm performing (blocked) 10-fold cross-validation of a several time series forecasting methods, measuring their mean squared error (MSE). I know that the MSE_cv is the average over the 10 MSEs. Is there a way to calculate the standard error as well? The usual SD/sqrt(n) formula probably do

Re: [R] Question concerning "pastecs" package

2007-05-25 Thread Philippe Grosjean
Hello, I already answered privately to your question. No, there is no translation of pastecs.pdf. The English documentation is accessible, as usual, by: ?turnpoints Regarding your specific question, 'info' is the quantity of information I associated with the turning points: I = -log2 P(t) w

Re: [R] Why might X11() not be found?

2007-05-25 Thread Patrick Connolly
On Fri, 25-May-2007 at 08:25AM +1200, Patrick Connolly wrote: |> > sessionInfo() |> R version 2.5.0 (2007-04-23) |> x86_64-unknown-linux-gnu |> |> locale: |> |> LC_CTYPE=en_US.UTF-8;LC_NUMERIC=C;LC_TIME=en_US.UTF-8;LC_COLLATE=en_US.UTF-8;LC_MONETARY=en_US.UTF-8;LC_MESSAGES=en_US.UTF-8;LC_PAPE

Re: [R] Make check failure for R-2.4.1

2007-05-25 Thread Adam Witney
Thanks for your replies Details inline below: On 24/5/07 17:12, "Martin Maechler" <[EMAIL PROTECTED]> wrote: >> "UweL" == Uwe Ligges <[EMAIL PROTECTED]> >> on Thu, 24 May 2007 17:34:16 +0200 writes: > > UweL> Some of these test are expected from time to time, since they are

Re: [R] Speeding up resampling of rows from a large matrix

2007-05-25 Thread Bill.Venables
Here is a possibility. The only catch is that if a pair of rows is selected twice you will get the results in a block, not scattered at random throughout the columns of G. I can't see that as a problem. ### --- start code excerpt --- nSNPs <- 1000 H <- matrix(sample(0:1, 120*nSNPs , replace=T),

Re: [R] Running R in Bash and R GUI

2007-05-25 Thread michael watson \(IAH-C\)
There are two things that occur. Firstly, I normally have to unset no_proxy: %> unset no_proxy; R Secondly, if for some reason http_proxy isn't being seen in R, you can use the Sys.putenv() function within R to manipulate the environment -Original Message- From: [EMAIL PROTECTED] [mailt