For MacOS we have
Binary packages, foo.tgz
Source packages, foo.tar.gz
Neither are `zip files' (things created by zip, usually with extension .zip).
It looks like you have not installed a source package before, and you
either do not have the development tools installed or they are not in your
On 9 Feb 2005 at 16:32, [EMAIL PROTECTED] wrote:
Petr,
It works perfectly! But I still have a question;
I have fit the following data;
x,y,z
1,10,11
2,11,15
3,12,21
4,13,29
5,14,39
6,15,51
7,16,65
8,17,81
9,18,99
10,19,119
dat.lm - lm(z~I(x^2)+y, data=dat)
dat.lm
Hi Bogdan
the behaviour is result of ct and lt POSIX class
try
str(test)
`POSIXlt', format: chr [1:4] 2005-02-08 18:49:15 2005-02-07
18:36:54 2005-02-04 18:37:03 2005-02-06 18:29:04
str(tst)
`POSIXct', format: chr [1:4] 2005-02-08 18:49:15 2005-02-07
18:36:54 2005-02-04 18:37:03 2005-02-06
Hello,
I ran my simulations on the Unix verson of R1.9.1 and the Windows version
of R2.0.1 on XP. I kept getting different values for the mean of the same
column of the same matrix, and am perplexed. I would appreciate if anyone
could help explain the difference?
Here is a sample code:
On Thu, 2005-02-10 at 08:01 +, Prof Brian Ripley wrote:
For MacOS we have
Binary packages, foo.tgz
Source packages, foo.tar.gz
Neither are `zip files' (things created by zip, usually with extension .zip).
It looks like you have not installed a source package before, and you
either
On Thu, 10 Feb 2005, Jari Oksanen wrote:
On Thu, 2005-02-10 at 08:01 +, Prof Brian Ripley wrote:
For MacOS we have
Binary packages, foo.tgz
Source packages, foo.tar.gz
Neither are `zip files' (things created by zip, usually with extension .zip).
It looks like you have not installed a source
Hi
Isn't it due to rounding and floating point precision of different
OSses?
Cheers
Petr
On 10 Feb 2005 at 3:44, Chang-Heok Soh wrote:
Hello,
I ran my simulations on the Unix verson of R1.9.1 and the Windows
version of R2.0.1 on XP. I kept getting different values for the mean
of the
Hello world,
I am actually transferring a course in data management for
students in biology, geography and agriculture
from statistica to R - it works
surprisingly well. If anyone is interested in my scratch/notepad
(in German language), please see
Hi Georg
If you know which value belongs to which year you can use
tapply, by, aggregate family like:
years-rep(c(2000,2001,2002,2003), each=10)
values-(1:40)
test-cbind(years,values)
test-data.frame(test)
aggregate(test$values,list(y=test$years),sum)
Cheers
Petr
On 10 Feb 2005 at 11:01,
Hello,
My problem is, that I have to build hundreds of GARCH models to obtain
volatility forecasts. I would like to run a loop, that would build those
forecasts for me. There is no problem, with writing only the results of
the forecasts, but I'd like to have stored results of the models in some
Dimitris Rizopoulos [EMAIL PROTECTED] writes:
maybe something like this:
dat - data.frame(years=rep(1994:2004, each=10), x=rnorm(110),
(I thought they abolished the 10 month year after the failure of the
French revolution. :-) )
y=rnorm(110))
lapply(split(dat[,-1], dat$years), cumsum)
or
On Thu, 10 Feb 2005, Wojtek Slusarski wrote:
Hello,
My problem is, that I have to build hundreds of GARCH models to obtain
volatility forecasts. I would like to run a loop, that would build those
forecasts for me. There is no problem, with writing only the results of
the forecasts, but
Can anyone explain why with latest version of R (2.0.1) on FC3, installed
from R-2.0.1-0.fdr.2.fc3.i386.rpm, update.packages() produces the message
/usr/lib/R/bin/Rcmd exec: INSTALL: not found.
Indeed /usr/lib/R/bin seems to lack various shell scripts (INSTALL,
REMOVE, etc).
This may or may not be related, but looking at the random number
generators in Numerical Recipies In C book it would appear that random
number generators with identical seeds will give different results
depending upon whther they are compliled with 64 bit or 32 bit
settings.
Is this the case
Roger Bivand wrote:
First, you should enclose the functions that output results in the loop in
print(), so that they do print, in loops this has to be done and is a FAQ
(often met with lattice graphics)
I am sorry for that, that I didn't find that. Now it works great. If
anyone needs that,
I M S White [EMAIL PROTECTED] writes:
Can anyone explain why with latest version of R (2.0.1) on FC3, installed
from R-2.0.1-0.fdr.2.fc3.i386.rpm, update.packages() produces the message
/usr/lib/R/bin/Rcmd exec: INSTALL: not found.
Indeed /usr/lib/R/bin seems to lack various shell scripts
That is correct, you have to specify with print or cat what you want to sink.
This should work for you:
print(summary(fit))
print(logLik(fit))
print(fcv)
I hope this helps.
On Thursday 10 February 2005 11:40, Wojtek Slusarski wrote:
Hello,
My problem is, that I have to build hundreds of
On Thu, 10 Feb 2005, Campbell wrote:
This may or may not be related,
It is completely unrelated: see earlier replies in this thread.
but looking at the random number
generators in Numerical Recipies In C book it would appear that random
number generators with identical seeds will give different
On Thu, 2005-02-10 at 13:52 +0100, Peter Dalgaard wrote:
I M S White [EMAIL PROTECTED] writes:
Can anyone explain why with latest version of R (2.0.1) on FC3, installed
from R-2.0.1-0.fdr.2.fc3.i386.rpm, update.packages() produces the message
/usr/lib/R/bin/Rcmd exec: INSTALL: not
Hi
Dumb question time again, for which I apologise.
I have a variable that contains the following numerical text 04010.
This is the name to access a list:
as.list(KEGGPATHID2NAME)$04010
[1] MAPK signaling pathway
Marvellous! Except I want to do that when 04010 is assigned to a
variable
The answer of course is parse()!
Thank you and good night!
M
-Original Message-
From: michael watson (IAH-C)
Sent: 10 February 2005 13:36
To: r-help@stat.math.ethz.ch
Subject: Using a number as a name to access a list
Hi
Dumb question time again, for which I apologise.
I have a
Quoting Jari Oksanen [EMAIL PROTECTED]:
On Thu, 2005-02-10 at 13:52 +0100, Peter Dalgaard wrote:
I M S White [EMAIL PROTECTED] writes:
Can anyone explain why with latest version of R (2.0.1) on FC3, installed
from R-2.0.1-0.fdr.2.fc3.i386.rpm, update.packages() produces the message
lists are like vectors, i.e.,
x - c(a = 1, b = 2)
path - b
x[path]
the same applies to lists:
lis - list(a = 1:3)
path - a
lis[path]
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address:
On Thu, Feb 10, 2005 at 01:35:56PM -, michael watson (IAH-C) wrote:
I have a variable that contains the following numerical text 04010.
This is the name to access a list:
as.list(KEGGPATHID2NAME)$04010
[1] MAPK signaling pathway
Marvellous! Except I want to do that when 04010 is
Hi all,
I would like to match two partitions. That is, if I have exactly the
same objects grouped
together for the two partitions, the labels may be arbitrarly permuted.
and so, i would like
to know the correspondances of the groups between the two clusterings.
In the e1701 pachage, it is
Dear R Users,
I have an xyplot() where different plotting symbols are used for
subgroups (originally used within S-Plus, but hopefully it´s also
applicable to R users).
How can I fit separate regression lines for every subgroup? So far, I
can only plot the overall fitted line.
The code looks
On Thursday 10 February 2005 09:34, Christoph Scherber wrote:
Dear R Users,
I have an xyplot() where different plotting symbols are used for
subgroups (originally used within S-Plus, but hopefully it´s also
applicable to R users).
How can I fit separate regression lines for every subgroup?
I have a working installation of RMySQL in Mac OS X (10.3.7). I have
libz in /usr/lib, as you do.
None of my personal installation notes mention doing anything special
for /usr/lib/libz*.
My R was built from a source distribution; it is not a binary
distribution. (version info below)
The only
Georg Hoermann [mailto:[EMAIL PROTECTED] wrote:
www.hydrology.uni-kiel.de/~schorsch/statistik/erle_stat.csv ...
contains a 10 year dataset.
We often need cumulative *annual* sums (sunshine, precipitation), i.e.
the sum must reset to 0
at the beginning of the year. I know of cumsum(), but I do
On Wed, 10 Feb 2005, Arin Basu wrote:
Hi All:
I am interested in tools available in R that will enable me to run
analysis of case crossover studies.
Conditional logistic regression, the basic tool, is the survival package.
-thomas
__
On Thu, 10 Feb 2005, michael watson (IAH-C) wrote:
The answer of course is parse()!
Nononono. If the answer is parse() you should usually rethink the
question.
The answer is as.list(WHATEVERITWASCALLED)[[path]]
-thomas
Thank you and good night!
M
-Original Message-
From: michael
Hi
I have some panel data for the 50 US states over about 25 years, and I
would like to test a simple model via OLS, using this data. I know how
to run OLS in R, and I think I can see how to create Panel Corrected
Standard Errors using
http://jackman.stanford.edu/classes/350C/pcse.r
What I
[EMAIL PROTECTED] writes:
Quoting Jari Oksanen [EMAIL PROTECTED]:
On Thu, 2005-02-10 at 13:52 +0100, Peter Dalgaard wrote:
I M S White [EMAIL PROTECTED] writes:
Can anyone explain why with latest version of R (2.0.1) on FC3,
installed
from R-2.0.1-0.fdr.2.fc3.i386.rpm,
In the nlme package you can find the gls() function to account for
autocorrelation over time using corAR1. Syntax might look something like
this:
fm1 - gls(response ~ IV, long, correlation=corAR1(form=~1|ID),
method='ML')
You can also use weights() for heteroscedasticity.
-Harold
-Original
Hi,
I was trying to compare differential expression data
from multiple studies ( 4 dataset downloaded from NCBI
GEO). The expression dataset are from Affy, as well as
cDNA (dual channel) platforms. Before I conduct
analysis, I have to map the gene IDs so that the
different genes are comparable.
Hei,
Olen lomalla ja palaan 16.02.2005.
Tilinavauksia hoitaa Annika Berglund 9 09-825 550.
Tilinavauksissa saattaa olla parin päivän viive.
T: Jaana Volotinen
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
Hi,
Thank you so much for the help. After setting the environment variables
I was able to install RMySQL. Thanks again.
Xander
On Thu, 2005-02-10 at 10:09, Don MacQueen wrote:
I have a working installation of RMySQL in Mac OS X (10.3.7). I have
libz in /usr/lib, as you do.
None of my
Pankaj Chopra wrote:
Hi,
I was trying to compare differential expression data
from multiple studies ( 4 dataset downloaded from NCBI
GEO). The expression dataset are from Affy, as well as
cDNA (dual channel) platforms. Before I conduct
analysis, I have to map the gene IDs so that the
different
In a multiple linear regression with two independent
variables is there any function in R to test for the
coefficients being different than some given values?
Example:
x1-rnorm(100)
x2-rnorm(100)
y-3+0.6*x1+0.3*x2
lm(y~x1+x2)
Obtain a test for the coefficients for x1 being
different than 0.6 and
You may also be interested in the following, which point out potentially
severe problems with such metanalyses.
http://www.cbs.dtu.dk/courses/norfa2004/Extras/5_1_Yves_TiG_meta-analysis_Ar
ticle.pdf
http://www.clinchem.org/cgi/content/full/47/8/1350
-- Bert Gunter
Genentech Non-Clinical
That's very helpful - I was on the point of giving up and going with
Stata! I will look into that in more detail. I assume that afterwards
it would be ok to apply the Beck and Katz procedure to get panel
corrected standard errors.
Cheers
David
On Thu, 10 Feb 2005 12:36:32 -0500, Doran, Harold
Hi,
I've got OSX 10.3 and R v2.0.1 and what I'd like to do is either take a
history file and run it as a source file (after I modify it in a text
editor) or simply output the command history to a plain text file that I
could then modify and later run as source.
The .history files are somewhat
Hi all,
We currently have 5 immediate openings in the Analytic Science/Product
Development organization of Fair Isaac in San Diego, CA, USA. Send me an
email if you are interested: [EMAIL PROTECTED]
Thanks much, and here's the blurb:
Dear Manuel,
Take a look at the linear.hypothesis function in the car package.
I hope this helps,
John
John Fox
Department of Sociology
McMaster University
Hamilton, Ontario
Canada L8S 4M4
905-525-9140x23604
http://socserv.mcmaster.ca/jfox
On 10 Feb 2005, at 19:26, Peter Dalgaard wrote:
[EMAIL PROTECTED] writes:
Quoting Jari Oksanen [EMAIL PROTECTED]:
On Thu, 2005-02-10 at 13:52 +0100, Peter Dalgaard wrote:
I M S White [EMAIL PROTECTED] writes:
Can anyone explain why with latest version of R (2.0.1) on FC3,
installed
from
Roger Bivand wrote:
On Thu, 10 Feb 2005, Georg Hoermann wrote:
Hello world,
based on the code of Roger I have now two solutions:
the first one (one line for the whole dataset)
- cut here --
erle -
read.csv(url(http://www.hydrology.uni-kiel.de/~schorsch/statistik/erle_stat.csv;))
jahre
Hi,
I run into a problem with R CMD Rd2dvi command: it gives me File `Rd.sty'
not found error (See the output message on the bottom). I get the same
error when running RCMD check. My system is: Windows 2000, R version
2.0.1, MiKTeX version 2.4. I do have a Rd.sty file in R_HOME/share/texm
Tuszynski, Jaroslaw W. JAROSLAW.W.TUSZYNSKI at saic.com writes:
:
: Hi,
:
: I run into a problem with R CMD Rd2dvi command: it gives me File `Rd.sty'
: not found error (See the output message on the bottom). I get the same
: error when running RCMD check. My system is: Windows 2000, R version
:
I rewrote scatter.smooth to handle missing values, but I have a question
about a move I had to make. Here's the code:
Mscatter.smooth-function (x, y, span = 2/3, degree = 1, family =
c(symmetric,
gaussian), xlab = deparse(substitute(x)), ylab =
deparse(substitute(y)),
ylim = range(y,
I really like these Trellis graphics but how do I get this code to
skip the missing?
logreg-read.csv(logreg.csv, header=TRUE, sep=,, na.string= )
attach(logreg)
bwplot(yesno~bc_pcb_tot |varlist, data=logreg, main=Box Cox PCB
transformation, auto.key=TRUE, fontfamily = HersheySans )
Lazy evaluation. See VR's S Programming for a good explanation. But, in
brief, the default values for xlab and ylab are not evaluated until they are
needed. If plot(..., xlab=...) appears after the assignments to x and y,
then x and y are no longer the original argument expressions but numeric
I am attempting to wrap the histogram function in my own custom
function, so that I can quickly generate some standard plots.
A part of what I want to do is to draw a normal curve over the histogram:
x - rnorm(1000)
hist(x, freq=F)
curve(dnorm(x), lty=3, add=T)
(for normal use, x would be a
Dear R-listers,
Does anybody know what is the correct source of rats dataset in survival
package?
The help gives the following information:
Rat data from survival5
Description:
48 rats were injected with a carcinogen, and then randomized to
either drug or placebo. The number of
Assuming I have years in YEAR and state ids in ID, I guess the
correlation ought to be
corAR1(form = ~ YEAR | ID)
?
Thanks a lot,
David
On Thu, 10 Feb 2005 12:36:32 -0500, Doran, Harold [EMAIL PROTECTED] wrote:
In the nlme package you can find the gls() function to account for
You need evaluate deparse(substitute(x)) _before_ you change x.
so have
xlab; ylab
early in the body of your function.
However, NEWS in R-devel says
o scatter.smooth() and loess.smooth() now handle missing values
in their inputs.
so why reinvent that wheel?
On Thu, 10 Feb
For which x do you think sd(x) is evaluated?
Hint: the help page shows
expr: an expression written as a function of 'x', or alternatively
the name of a function which will be plotted.
and you have written dnorm(x, mean=mean(x), sd=sd(x)) as function of x.
On Thu, 10 Feb 2005, Thomas
Thomas Hopper [EMAIL PROTECTED] writes:
I am attempting to wrap the histogram function in my own custom
function, so that I can quickly generate some standard plots.
A part of what I want to do is to draw a normal curve over the histogram:
x - rnorm(1000)
hist(x, freq=F)
Are the years equally spaced and in time order?
If so, it probably doesn't matter, and if not you may want corCAR1 not
corAR1.
On Thu, 10 Feb 2005, David Hugh-Jones wrote:
Assuming I have years in YEAR and state ids in ID, I guess the
correlation ought to be
corAR1(form = ~ YEAR | ID)
?
Thanks a
Just to explain my previous mail, here is the output I get.
dim(tissue.exp)
[1] 1532 20
pick-sample(tissue.exp,5,replace=TRUE)
dim(pick)
[1] 15325
tissue.exp.t-t(tissue.exp)
dim(tissue.exp.t)
[1] 20 1532
pick-sample(tissue.exp.t,5,replace=TRUE)
dim(pick)
NULL
Thanks
Prof Brian Ripley ripley at stats.ox.ac.uk writes:
:
: This is a problem with a recent-ish change to MikTeX to make it use
: non-standard paths: see the notes on Duncan Murdoch's site at
:
:http://www.murdoch-sutherland.com/Rtools/miktex.html
:
: Fptex is a standard distribution e.g.
Gonzalez Ruiz, Juan Ramon jrgonzalez at ico.scs.es writes:
:
: Dear R-listers,
:
: Does anybody know what is the correct source of rats dataset in survival
package?
:
: The help gives the following information:
:
: Rat data from survival5
: Description:
: 48 rats were injected with a
I am running R 2.01 on Mac OS 10.3.7. Whenever I save graphs, they are
saved as PDF files. Are there any other file formats to which I could save
my graphs. I would like to directly export my graphs to MS Word, if
possible.
Thank you in advance.
Sincerely,
Jim Milks
Graduate Student
Two time series questions:
FITTING TRANSFER FUNCTIONS WITH LAGS: Consider the following toy example:
dates - paste(11:21, /01/2005, sep=)
Dates - as.Date(dates, %d/%m/%Y)
set.seed(1)
DF - data.frame(date=Dates, y=rnorm(11), x=rnorm(11, 3))
arima(DF$y, c(1,0,0), xreg=lag(DF$x, 1))
Hello,
I want to build some tables in my project using R, does R have some tables form
that I can use? if i use the internal R table, how can initial a table before I
use it?
I want to my tables to have some columns, for example:
(trial, x, y, goal) I want to put these columns into my tables,
## The following corrects what I believe to have been an error in my
previous post:
Two time series questions:
I. FITTING TRANSFER FUNCTIONS WITH LAGS: Consider the following toy
example:
dates - paste(11:21, /01/2005, sep=)
Dates - as.Date(dates, %d/%m/%Y)
set.seed(1)
DF -
Howdy R gurus !
I am newbie to R
I use R 2.0.1 in Windows XP. When I run R
I got the follwoing memory error.
My physical memory size is 3 Gb.
My R got the memory problem when it reached to
about 2 Gb.
Thanks in advance,
library(spdep)
sfr.lagsarlm - lagsarlm(sfr.data$Bldgsqft ~
I'm having trouble using the read.mtp function in the foreign package
to import datasets from minitab (.mtw) format. Specifically, each file
I try to import fails to load any data beyond a header row stating the
version of Minitab that saved the dataset. I get this error:
incomplete final line
You need to read the rw-FAQ (as the posting guide asks) and make the R
executable capable of using more than 2Gb (assuming that you have already
enabled your OS).
On Fri, 11 Feb 2005, Kum-Hoe Hwang wrote:
Howdy R gurus !
I am newbie to R
I use R 2.0.1 in Windows XP. When I run R
I got the
Hi
On 10 Feb 2005 at 18:10, Cuichang Zhao wrote:
Hello,
I want to build some tables in my project using R, does R have some
tables form that I can use? if i use the internal R table, how can
initial a table before I use it?
How do you want to use a table???
I want to my tables to have
Hello,
I am encountering some problems when overloading the [ operator for a
new S4 class based on array. This is an example class definition:
setClass(foo,
representation(array),
prototype(array(NA, dim=c(3,3)),
dimnames=list(age=1:3, year=10:12))
)
And this the
Hi Peter
On 11 Feb 2005 at 1:06, Peter Jaffe wrote:
I'm having trouble using the read.mtp function in the foreign package
to import datasets from minitab (.mtw) format. Specifically, each file
I try to import fails to load any data beyond a header row stating the
version of Minitab that
Jim and Chana Milks wrote:
I am running R 2.01 on Mac OS 10.3.7. Whenever I save graphs, they are
saved as PDF files. Are there any other file formats to which I could save
my graphs. I would like to directly export my graphs to MS Word, if
possible.
Thank you in advance.
Sincerely,
Jim Milks
72 matches
Mail list logo