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
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
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
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
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
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
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
>
>
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
[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
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
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
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
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
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
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
[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
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
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=
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
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
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
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
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
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
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
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@
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.
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
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
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
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
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
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
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
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 :-(
>
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,
>
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-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.
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
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]>
>>
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
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
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
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
> 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
[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
> "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
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.
>
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-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
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
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
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
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
> 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 ==
> "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
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
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
>> ##-- 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
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
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
> "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]>
>
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
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
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
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
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
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
--- [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/..
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
> 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
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
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.
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
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
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
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
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
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
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),
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
81 matches
Mail list logo