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 u
Heather J. Branton 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
verify
Yasser El-Zein 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 could not m
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
glmmPQ
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
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
ht
> 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
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",
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
>>
>>
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 e
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
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
> -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 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 dup
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,
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 pn
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 B
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
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
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
predict.lm
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.. This does raise the question as to whether gregmisc and other
bundles should appear on the "install packages from CR
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 PRO
"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
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 to
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 r
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 t
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()
+ install.packages(select.list(a[,
On Wed, 24 Nov 2004 14:38:19 -0500, "Heather J. Branton" <[EMAIL PROTECTED]>
wrote :
>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)
Yes, that's right.
Duncan Murdoch
__
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 wit
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 Rm
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))
>
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 :
D
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')
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 ti
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 h
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 looki
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 thi
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
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] mai
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.
>
> TH
"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 ~ orde
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 strugglin
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 <- sum(norm, na.rm=T)
binconf(s, n)
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
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 comma
On Wed, 2004-11-24 at 17:43 +0100, Witold Eryk Wolski wrote:
> 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:
>
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 'he
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: W
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. w
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 sma
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 mode
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 featu
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 most
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 view
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 BioC
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 scatterpl
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) co
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
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 El
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 documen
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. w
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,
tsdiag(ar(lh)
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 dataset
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...
__
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 fr
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 to
On 24 Nov 2004, at 10:16, Christoph Lehmann wrote:
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
###
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í.
http://www.correodir
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 n
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
lp00/archive|Rhelp01/archive|Rhelp02a/archive;forma
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
On 24 Nov 2004 13:59:41 +0100, Peter Dalgaard
<[EMAIL PROTECTED]> wrote :
>Duncan Murdoch <[EMAIL PROTECTED]> writes:
>> Here's a strategy that I hope subverts this irritating mechanism:
>> Every now and then I get a challenge about a message that I didn't
>> send, because someone (or some virus)
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 s
On 24-Nov-04 Martin Maechler wrote:
>> "Ted" == Ted Harding <[EMAIL PROTECTED]>
> [...]
> Ted> [...] The "Challenge"
> Ted> mechanism would destroy the mailing-list community if
> Ted> it became widely adopted.
>
> Exactly.
> I've received such a message myself from the same "mach
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 mi
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,
> >
> "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
Gabor Grothendieck myway.com> writes:
:
: Yasser El-Zein 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
: : myway.com> wrote:
: : > Yasser El-Zein gma
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
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 ge
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 (file
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){
na
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 form
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 con
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++ g
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)
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 <- paste("pair",i,se
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 th
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 PROT
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)
-- OR
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 <- glmmPQL(ch
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- list
or if it belongs to R-h
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 s
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
___
[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
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 wro
[EMAIL PROTECTED] wrote:
By 'ignore', can we delete those from the list of data? I would then
assume that if you have a sequence of +0+0+ that you would want the last
"+" for the increase of three.
If that is the case, then do a 'diff' and delete the entries that are 0.
Then create a new 'diff' a
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 indi
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
https://stat.ethz.ch/mailman/listinf
99 matches
Mail list logo