On Fri, 29 Apr 2005, Patrick Connolly wrote:
On Thu, 28-Apr-2005 at 05:25PM -0400, Xiang-Jun Lu wrote:
|> Hi,
|>
|> I have just noticed the following problem with R2.1.0 running on SuSE 9.1,
|> [However, version 2.0.1 (2004-11-15) on the same machine works Okay]:
|>
|> -
Do you mean that you have one file where the data is ordered by
group? If so, using scan and a loop is an easy option. If you have
the data in separate files by group, just make a character vector
with your filenames as elements and use a loop to read the data using
read.table at each iteration
Hello
Instead of reading in group1.txt I want to read in groups1 for the first
iteration of i, then groups2 for the second and so on. Obviously I can't use
groups(i) but assume there is a way to do this.
group<-read.table("C:/Data/April 2005/group1.txt",header=T)
thanks in advance
Meredith
_
I have snippets of code that I have either taken from examples or off of the
list. I apologise to those I have stolen it from but I didn't keep the proper
references.
n <- 100
xx <- c(0:n, n:0)
yy <- c(c(0,cumsum(rnorm(n))), rev(c(0,cumsum(rnorm(n)
plot (xx, yy, type="n", xlab="Time", yla
Just a little bit of trivia.
The 2001 Census in Australia had a significant group of people who responded to
the question of religion with the answer Jedi or Jedi Knight. Unfortunately the
Australian Bureau of Statistics is a bit fuddy duddy about the issue as they
see it as a trivialisation of
On Thu, Apr 28, 2005 at 05:11:11PM +0100, V Nyirongo wrote:
> Dear all,
>
> Am new on this list but need help:
>
> I have a fortran code which I compile into a dynamic library (on Solaris).
>
> My probem is that when I call this code soon after starting R, it runs ok.
> But it doesn't for the s
On 4/28/05, Jonathan Q. <[EMAIL PROTECTED]> wrote:
>
> just started using r. a few questions I can't figure out.
>
> first, when loading in a text file, I can do it from the web but not
> from hard drive.
> Fil<- 'http://www.test.com/file.txt' -- works
> Fil<- 'c:\program files\file.txt' --does
Dear Ed,
I think the discrepancy you are seeing is due to a change in
parameterization of the inverse Wishart distribution in MCMCpack from
MCMCpack version 0.6-1 to 0.6-2 and later. Starting in 0.6-2 we've
parameterized the inverse Wishart in the same way as in Appendix A of
Gelman, Carlin, St
just started using r. a few questions I can't figure out.
first, when loading in a text file, I can do it from the web but not
from hard drive.
Fil<- 'http://www.test.com/file.txt'-- works
Fil<- 'c:\program files\file.txt' --does not. what am i doing wrong??
also, how do you plot a graph of
On Thu, 2005-04-28 at 16:00 -0500, Ghosh, Sandeep wrote:
> Hi,
>
> I'm having problems getting the resulting R image generated as an png
> file when trying to feed the R cmds through a java prog to R. I was
> facing the same problem once before and Deepayan had suggested to use
> "print" along wit
Dear R
I'm trying to plot the lda means onto a 2 D plot of discriminant scores.
Preferably I'd like these to be in a larger font compared to the
discriminant scores.
I tried
skull.mean.pred <- predict(skulls.lda, as.data.frame(skulls.lda$means),
dimen=2)
from which I got
skull.mean.pred
$clas
Xiang-Jun Lu <[EMAIL PROTECTED]> writes:
> Hi,
>
> I have just noticed the following problem with R2.1.0 running on SuSE
> 9.1, [However, version 2.0.1 (2004-11-15) on the same machine works
> Okay]:
> -
> > hist(rnorm(100))
On Thu, 28-Apr-2005 at 05:25PM -0400, Xiang-Jun Lu wrote:
|> Hi,
|>
|> I have just noticed the following problem with R2.1.0 running on SuSE 9.1,
|> [However, version 2.0.1 (2004-11-15) on the same machine works Okay]:
|>
|> --
At last, simple addition in R is now fast and easy!!!
They said it could never be done, but they've never heard of Fortran.
In order to add two numbers in R before you had to type:
> 7 + 3
It was tedious to type and several hours could pass as you waited for
your return value.
So, imagine the
You are passing just a string to subset(). At the very least you need
to parse it (but still this does not work easily with subset() -- see
below). But are you sure you need to do this? subset() for dataframes
already accepts subset expressions involving the columns of the
dataframe, e.g.:
"Pascal Boisson" <[EMAIL PROTECTED]> writes:
> Hello all,
>
> I have some trouble in reconstructing a valid expression within a
> function,
> here is my question.
>
> I am building a function :
>
> SUB<-function(DF,subset=TRUE) {
> #where DF is a data frame, with Var1, Var2, Fact1, Fact2, Fact3
Hi,
I have just noticed the following problem with R2.1.0 running on SuSE 9.1,
[However, version 2.0.1 (2004-11-15) on the same machine works Okay]:
-
hist(rnorm(100))
Error in title(main = main, sub = sub, xlab = xlab, ylab
Hello all,
I have some trouble in reconstructing a valid expression within a
function,
here is my question.
I am building a function :
SUB<-function(DF,subset=TRUE) {
#where DF is a data frame, with Var1, Var2, Fact1, Fact2, Fact3
#and subset would be an expression, eg. Fact3 == 1
#in a first
Hi,
I'm having problems getting the resulting R image generated as an png file when
trying to feed the R cmds through a java prog to R. I was facing the same
problem once before and Deepayan had suggested to use "print" along with the
chart cmd like print(barchart(...)), but the print cmd was n
Hi all,
I'm writing a vignette for my package, and I would like to include some of
the package help files in there as well. Is there an easy way of doing so?
I tried using R CMD Rdconv to generate latex files from .Rd files but I am
not sure how to include these into a .Rnw file (ie the vignette so
Here is summary of
l<-qqnorm(kk) # kk is my sample
l$y (which is my sample)
l$x (which is therotical quantile)
diff<-l$y-l$x
and
> summary(l$y)
Min. 1st Qu. MedianMean 3rd Qu.Max.
0.9007 0.9942 0.9998 0. 1.0060 1.1070
> summary(l$x)
Min.1st Qu. Median Me
Hi Alex,
The binary, from the Mac OS website, works on my system (R-2.1.0):
library(sna)
help(package=sna)
example(bbnam)
I've not tried to install the source version. Let me know if this works
for you.
Rob
On Apr 28, 2005, at 1:03 PM, Alex Bach wrote:
Hi folks,
I am trying to install the sna pac
Hi folks,
I am trying to install the sna package on a MacOS X 10.3.9 and
10.3.4. In neither case I have been able to get the R CMD Install
command to succeed from the Terminal window.
With 10.3.9 I get an error with the Make command, and with 10.3.4 I
get an error on the Description(???)
Has
Hua Li <[EMAIL PROTECTED]> writes:
> In R, is there any way that I can let R not printing out the "Read 4
> items " message?
>
> > a<-scan("probes.txt")
> Read 4 items
Yes, and ?scan tells you how to make it quiet almost immediately!
--
O__ Peter Dalgaard Blegdamsvej 3
Well, you could always wrap it in
sink(...)
scan(...)
sink()
I suppose.
-- Bert Gunter
Genentech Non-Clinical Statistics
South San Francisco, CA
"The business of the statistician is to catalyze the scientific learning
process." - George E. P. Box
> -Original Message-
> From: [EMA
use quiet=T
help("scan")
norm
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Hua Li
Sent: Thursday, April 28, 2005 12:32 PM
To: r-help@stat.math.ethz.ch
Subject: [R] environment variables
In R, is there any way that I can let R not printing out the "
Totally understood. If you remain "on vacation" (we wish eh?!) and
get further messages like that from the list s'ware, just accept my
standing apologies and delete them.
Very welcome re letter!
Very best,
Chris
--- Forwarded message follows ---
Send reply to: "Alana O'C" <[
In R, is there any way that I can let R not printing out the "Read 4
items " message?
> a<-scan("probes.txt")
Read 4 items
Thanks
Hua LI
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting g
Wayne Fuller's Measurement Error Models is a good reference.
url:www.econ.uiuc.edu/~rogerRoger Koenker
email [EMAIL PROTECTED] Department of Economics
vox:217-333-4558University of Illinois
fax:217-244-6678
Bert wrote:
>>You probably have noticed that
>> I'm quite new
>> to statistics, but I'm working on that...
>>
>>
> And you want to use Bayesian methods?!
>
I was always under the impression that it's mostly a matter of mindset if
you go Bayesian or frequentist, not of your statistical skills.
[..
Stock returns and other financial data have often found to be heavy-tailed.
Even Cauchy distributions (without even a first absolute moment) have been
entertained as models.
Your qq function subtracts numbers on the scale of a normal (0,1)
distribution from the input data. When the input data are
On Thu, 2005-04-28 at 13:53 -0400, Andy Bunn wrote:
> > PS 6 days to the big Jedi holiday!
>
> I wonder if anybody who gave the matter any thought would be surprised that
> the R-Help list is populated by ubergeeks.
Perhaps, but we all managed to miss Pi Day last month
> Pi.Day <- as.POSIXct
- Original Message -
From: "Berton Gunter" <[EMAIL PROTECTED]>
To: "'Pieter Provoost'" <[EMAIL PROTECTED]>;
Sent: Thursday, April 28, 2005 6:26 PM
Subject: RE: [R] normality test
>
>
> > -Original Message-
> > From: [EMAIL PROTECTED]
> > [mailto:[EMAIL PROTECTED] On Behalf Of Pi
Could someone please help me by giving me a reference to how one computes
standard errors for the coefficients in an orthogonal linear regression, or
perhaps someone has some R code? (I would accept a derivation or formula, but
as a former teacher, I know how that can rankle.) I tried to imitate
If you're trying to find the textual form of an actual argument, here's
one way:
> foo <- function(x) {
+ xn <- substitute(x)
+ if (is.name(xn) && !exists(as.character(xn)))
+ as.character(xn)
+ else
+ x
+ }
> foo(x)
[1] 3
> foo(xx)
[1] "xx"
> foo(list(xx))
Error in fo
> PS 6 days to the big Jedi holiday!
I wonder if anybody who gave the matter any thought would be surprised that
the R-Help list is populated by ubergeeks.
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do r
Dear R-helpers:
I pointed out my question last time but it is only partially solved.
So I would like to point it out again since I think it is very
interesting, at least to me.
It is a question not about how to use R, instead it is a kind of
therotical plus practical question, represented by R.
I
On Thu, 28 Apr 2005, Omar Lakkis wrote:
r
[1] open settle
<0 rows> (or 0-length row.names)
class(r)
[1] "data.frame"
this is an empty data.frame I get back from a sql statement that
returns an empty result set. How can I create such an empty data.frame
using the data.frame() constructor?
I want t
On Thu, 2005-04-28 at 19:02 +0200, Achim Zeileis wrote:
> On Thu, 28 Apr 2005 11:50:29 -0500 Marc Schwartz wrote:
>
> > On Thu, 2005-04-28 at 18:33 +0200, Achim Zeileis wrote:
> > > On Thu, 28 Apr 2005 12:22:51 -0400 Rajarshi Guha wrote:
> > >
> > > > On Thu, 2005-04-28 at 17:01 +0200, Henric Nil
On 4/28/05, Omar Lakkis <[EMAIL PROTECTED]> wrote:
>
> > r
> [1] open settle
> <0 rows> (or 0-length row.names)
> > class(r)
> [1] "data.frame"
>
> this is an empty data.frame I get back from a sql statement that
> returns an empty result set. How can I create such an empty data.frame
> using th
There is a force function in the Base package!
Does that counteract the evil of the try function?
Erin
__
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/post
Marc Schwartz wrote:
Shouldn't that be more like:
"Much to learn you have still, young apprentice mine"
R offends my Jedi religious sensibilities by having a 'try' function.
[see: http://www.starwars.com/databank/character/yoda/ about four paras
from bottom]
Baz
PS 6 days to the big Jedi holida
> r
[1] open settle
<0 rows> (or 0-length row.names)
> class(r)
[1] "data.frame"
this is an empty data.frame I get back from a sql statement that
returns an empty result set. How can I create such an empty data.frame
using the data.frame() constructor?
I want to have a data.frame with 0 rows but
Achim Zeileis wrote:
On Thu, 28 Apr 2005 08:52:33 -0500 roger koenker wrote:
For my money, Frank's comment should go into fortunes. It seems a
rather Sisyphean battle to keep the lessons of robustness on the
statistical table but nevertheless well worthwhile.
Added.
On more comment: maybe it's
On Thu, 28 Apr 2005 11:50:29 -0500 Marc Schwartz wrote:
> On Thu, 2005-04-28 at 18:33 +0200, Achim Zeileis wrote:
> > On Thu, 28 Apr 2005 12:22:51 -0400 Rajarshi Guha wrote:
> >
> > > On Thu, 2005-04-28 at 17:01 +0200, Henric Nilsson wrote:
> > > > Heather Joan Lynch said the following on 2005-04
On Thu, 2005-04-28 at 18:33 +0200, Achim Zeileis wrote:
> On Thu, 28 Apr 2005 12:22:51 -0400 Rajarshi Guha wrote:
>
> > On Thu, 2005-04-28 at 17:01 +0200, Henric Nilsson wrote:
> > > Heather Joan Lynch said the following on 2005-04-28 16:14:
> >
> > > First, citing Simon `Yoda' Blomberg, "This is
R users-
In moving from R 2.0.0 to R 2.1.0 in Windows, I have encountered a problem
with the "riwish" command in the package "MCMCpack". I've searched the
documentation and can't seem to figure it out. For example:
Define a matrix:
> lam <- matrix(c(.00233,-.00057,-.00057,.00190),2,2)
and then
On Thu, 28 Apr 2005 12:22:51 -0400 Rajarshi Guha wrote:
> On Thu, 2005-04-28 at 17:01 +0200, Henric Nilsson wrote:
> > Heather Joan Lynch said the following on 2005-04-28 16:14:
>
> > First, citing Simon `Yoda' Blomberg, "This is R. There is no if.
> > Only how."
>
> I hope this quote makes into
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] On Behalf Of Pieter Provoost
> Sent: Thursday, April 28, 2005 7:52 AM
> To: R-help@stat.math.ethz.ch
> Subject: Re: [R] normality test
>
> Thanks all for your comments and hints. I will try to keep
> them in mi
On Thu, 2005-04-28 at 17:01 +0200, Henric Nilsson wrote:
> Heather Joan Lynch said the following on 2005-04-28 16:14:
> First, citing Simon `Yoda' Blomberg, "This is R. There is no if. Only how."
I hope this quote makes into the fortunes package :)
--
Dear all,
Am new on this list but need help:
I have a fortran code which I compile into a dynamic library (on Solaris).
My probem is that when I call this code soon after starting R, it runs ok. But
it doesn't for the second time without exiting R first. In this code I have to
generate some rando
Below.
>
> Usually (but not always) doing tests of normality reflect a lack of
> understanding of the power of rank tests, and an assumption of high
> power for the tests (qq plots don't always help with that because of
> their subjectivity). When possible it's good to choose a
> robust m
Hi
Do you know how I can do a mixed factorial manova with R ?
Thanks
Julien Martin
Laboratoire de Denis Réale
Chaire de Recherche du Canada en écologie comportementale
Canadian Research Chair in behavioural ecology
Département des sciences biologiques
Université du Québec à Montréal
Case postal
On Thu, 2005-04-28 at 10:14 -0400, Heather Joan Lynch wrote:
> Hello,
>
> I cannot figure out how to shade between two lines in a plot. For
> example, if I am trying to plot a confidence envelope and I would like to
> shade the interior of the envelope grey. Is this possible in R?
>
> Thanks in
Heather Joan Lynch said the following on 2005-04-28 16:14:
I cannot figure out how to shade between two lines in a plot. For
example, if I am trying to plot a confidence envelope and I would like to
shade the interior of the envelope grey. Is this possible in R?
First, citing Simon `Yoda' Blomber
Thanks all for your comments and hints. I will try to keep them in mind.
Since a number of people asked me what I'm trying to do: I want to apply
Bayesian inference to a simple ecological model I wrote, and therefore I
need to fit (uniform, normal or lognormal) distributions to sets of observed
dat
Hello,
I cannot figure out how to shade between two lines in a plot. For
example, if I am trying to plot a confidence envelope and I would like to
shade the interior of the envelope grey. Is this possible in R?
Thanks in advance.
Heather Lynch
__
R-
On Thu, 28 Apr 2005 08:52:33 -0500 roger koenker wrote:
> For my money, Frank's comment should go into fortunes. It seems a
> rather Sisyphean battle to keep the lessons of robustness on the
> statistical table but nevertheless well worthwhile.
Added.
On more comment: maybe it's also worth no
Dear Sir,
My name is Javier Bussi and I´m a professor at the University of Rosario,
Argentina. The National Household Survey has changed the periodicity of the
information provided on unemployment from biannual to quarterly. It is not
the usual situation where we have two series at the same period
Dear listers,
After activating the name space for my bioconductor package (prada) I
successfully ran R CMD check. However when loading the package in R and
running the examples the imported function brewer.pal from package
RColorBrewer is not found. I can directly call brewer.pal from the
RColorBre
For my money, Frank's comment should go into fortunes. It seems a
rather Sisyphean battle to keep the lessons of robustness on the
statistical table
but nevertheless well worthwhile.
url:www.econ.uiuc.edu/~rogerRoger Koenker
email [EMAIL PROTECTED] De
Thanks for the advice!
For those interested in what the final solution was to installing RMySQL (R
2.1.0, Linux SuSE 9.3):
The following packages were not installed:
php4-mysql
php5-mysql
php5-mysqli
mysql-devel
at least one of them was obviously missing
On Thursday 28 April 2005 10:37, yo
Romain Francois wrote:
Le 28.04.2005 13:16, Pieter Provoost a écrit :
Hi,
I have a small set of data on which I have tried some normality tests.
When I make a histogram of the data the distribution doesn't seem to
be normal at all (rather lognormal), but still no matter what test I
use (Shapiro,
- Original Message -
From: "Romain Francois" <[EMAIL PROTECTED]>
To: "Pieter Provoost" <[EMAIL PROTECTED]>; "RHELP"
Sent: Thursday, April 28, 2005 2:03 PM
Subject: Re: [R] normality test
> Le 28.04.2005 13:16, Pieter Provoost a écrit :
>
> >Hi,
> >
> >I have a small set of data on which
Le 28.04.2005 13:16, Pieter Provoost a écrit :
Hi,
I have a small set of data on which I have tried some normality tests. When I make a histogram of the data the distribution doesn't seem to be normal at all (rather lognormal), but still no matter what test I use (Shapiro, Anderson-Darling,...) it
On Thu, 28 Apr 2005, Crispin Miller wrote:
Im' using matrix() and layout() to set up a fairly complex plot...
I was wondering if anyone can point me to something that describes the
appropriate incantations to determine, given a text string, the maximum
font size that will allow it to fit into the c
Hello, Sebastian.
It looks like you are also missing MySQL API installed.
It consists of several C headers and libraries with definitions of
functions, working with mysql server daemon.
In your case rackage installation process needs to compile some C sources, and
a compiler must know w
Hi,
I have a small set of data on which I have tried some normality tests. When I
make a histogram of the data the distribution doesn't seem to be normal at all
(rather lognormal), but still no matter what test I use (Shapiro,
Anderson-Darling,...) it returns a very small p value (which as far
check this:
mat <- matrix(rnorm(5*5), 5, 5)
mat[1, 3] <- 1; mat[4, 2] <- 2
which(mat==1, TRUE)
which(mat==2, TRUE)
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuven, Belgium
Tel:
Dear useRs,
a new version of the lmtest package is available from CRAN. Thanks to
Giovanni Millo who provided the initial versions of several new
functions (and many helpful discussions), there is new functionality for
the comparison of nested and non-nested linear models.
For non-nested model co
Dear useRs,
recently, there were several discussions on R-help about how to
conveniently fit dynamic linear models and time series regressions. The
package dynlm tries to address this problem by 1. providing some more
functions like lags L() and differences d() and season() in the formula
specific
Joao:
1) The error message you get when setting nu=0 is due to the fact that
no support vectors can be found with that extreme restriction, and this
confuses the predict function (try svm(, fitted = false): the model
returned is empty). In fact, the C++ code interfaced by svm() clearly
allows
Dear useRs,
a new version of the zoo package for totally ordered observations is
available from CRAN. It contains many new features, the most important
of which are: support of regular time series, an improved merge method
and functions for carrying out rolling analyses of time series. A more
deta
Hi
Ive got big matrixes with just few non missing data, and
would like to use or do a funtion to find the index of the row and col of one
specified value. For example I just have one "1" value in the matrix, and one
"2" value in the same matrix.
Has anyone an idea about how to extract the rowind
Dear all,
Im' using matrix() and layout() to set up a fairly complex plot...
I was wondering if anyone can point me to something that describes the
appropriate incantations to determine, given a text string, the maximum
font size that will allow it to fit into the current plotting area
without be
75 matches
Mail list logo
