Hi Laurent,
You should look the R windows FAQ web page:
(a) Use the alternative |internet2.dll| by starting R with the flag
--internet2 (see How do I install R for Windows?
http://cran.r-project.org/bin/windows/base/rw-FAQ.html#How-do-I-install-R-for-Windows_003f)
which uses the Internet
Please do read the R posting guide and supply essential information such
as your OS (and not send HTML mail).
If you are using Windows, see the rw-FAQ Q2.19. Note that R never asks
you for `an active directory login/password': it may *if you ask it to*
ask for proxy identification. Quite
At 12:55 07.05.2006 +0100, Prof Brian Ripley wrote:
On Sun, 7 May 2006, Heinz Tuechler wrote:
Hello Xiaochun Li!
Thank you for submitting the function. At the time I had that problem I
solved it in a somewhat different way.
I changed a few lines in the print.survfit method. I introduced a
Hi Jinsong,
you can use the package drc on CRAN to fit logit and Weibull models
(but not the probit model) with natural mortality/response and/or
natural immunity to binomial data (maximum likelihood estimation).
To get an idea try:
library(drc)
?earthworms
Christian
I want to label groups of points in
the plot of residuals against predicted values
generated by plot.lm().
Labelling the points by color
x-1:15
y-rnorm(15)
f-factor(rep(c(1:3),each=5))
g-glm(y~x)
plot(g,which=1,col=c(1,2,3)[f])
I receive the error message
Error in title(main, sub, xlab,
David Foreman wrote:
Following upgrading to R-2.3.0 SciViews generates 2 warning messages (I
suspect 1 for each dock) saying 'null environment deprecated', and the R
Commander menu simply returns a 'file-name not found' error. I've had
to go back to 2.2.1 for the moment (where everything
The solution is to use a separate function to do this (and call it from the
print method).
On Mon, 8 May 2006, Heinz Tuechler wrote:
At 12:55 07.05.2006 +0100, Prof Brian Ripley wrote:
On Sun, 7 May 2006, Heinz Tuechler wrote:
Hello Xiaochun Li!
Thank you for submitting the function. At
Hello Anne,
the function 'ur.pp' contained in the package 'urca' has been ported
into the package 'fSeries' (see the documentaion of ?urppTest).
help(ur.pp-class, package=urca) will tell you that p-values are not
part of this class, i.e., these are not computed but critical values are
returned.
thank you Martin , but the problem now is that when I run R CMD INSTALL -c
SJava.tar .gz I have the message perl n'est pas reconnue en tant que command
interne ou externe
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch
On Mon, 8 May 2006, Ulrich Halekoh wrote:
I want to label groups of points in
the plot of residuals against predicted values
generated by plot.lm().
Labelling the points by color
x-1:15
y-rnorm(15)
f-factor(rep(c(1:3),each=5))
g-glm(y~x)
plot(g,which=1,col=c(1,2,3)[f])
I receive the
hello
I would like to know which package contains stepAIC
function ?
regards
--
Ahmet Temiz
Jeoloji Müh.
Afet İşleri Genel Müdürlüğü
Deprem Araştırma Dairesi
Tel: (312) 287 89 51 veya (312) 287 26 80/1547
Faks: (312) 287 89 51
E. Posta: [EMAIL PROTECTED]
www.deprem.gov.tr
Ahmet Temiz
miniar mansouri wrote:
thank you Martin , but the problem now is that when I run R CMD INSTALL -c
SJava.tar .gz I have the message perl n'est pas reconnue en tant que
command interne ou externe
[[alternative HTML version deleted]]
So what about installing perl on your machine - or
help.search(stepAIC) and RSiteSearch(stepAIC, restrict=functions)
show it is in package MASS.
orkun wrote:
hello
I would like to know which package contains stepAIC
function ?
regards
--
Chuck Cleland, Ph.D.
NDRI, Inc.
71 West 23rd Street, 8th floor
New York, NY 10010
tel: (212)
Hi All,
I,m trying to do a GARCH simulation in R 2.3.0 release
in Windows XP. I've seen garchsim function but that is
for garch (1,1) and ?garch gives an example for ARCH
simulation. Can anyone help me how can i extend the
help shown in ?garch to GARCH simulation? Please help
me in this regard.
Hello,
Can someone point me to documentation or ideas on how to calculate the
centroids of clusters identified with hclust ?
I would like to be able to chose the number of clusters (in the style of
cutree) and then get the centroids of these clusters.
This seems like a quite obvious task to
I use the R-2.3.0 with sciViews in windows XP homw
now. Mostly it is fine. But every booting, sciView
failed for the first time and said that it could not
establish connection with R. After closing
everything, it works fine since. I hope these
information is helpful. It is a great software
Steve Buysek Wrote
Using only informative incomplete trios and dropping others turns out to
bias the TDT (Curtis Sham 95). There are a variety of modifications to
handle case-parent trios, starting (I think) with Weinberg's 1999 paper that
uses the EM algorithm and a likelihood-ratio test. I
UweL == Uwe Ligges [EMAIL PROTECTED]
on Fri, 05 May 2006 17:50:22 +0200 writes:
UweL r user wrote:
I want to check if an object (dataset, vector, etc) is
?present?. If it is present, I will do nothing.
If it is not present, I will load it from my hard drive.
Dear R Users,
Here is another probelm/question.
I would like to run some panel regressions with R. Therefore I have
combined several time periods of data for different individuals in my
database.
I have already run pooled OLS but I would need to calculate a Fixed
Effects Estimator (within
I have a problem with the export of a simulated matrix (nrow=1000, ncol=492). I
want to obtain a .txt file with fixed format data (five character for example).
I've tried to use the instructions: options(digits=5) followed by
write.table(...), but the .txt file reports always the data with
Dear R-help list members
I gathered longitudinal data on fish behaviour which I try to analyse using
a multi level model for change. Mostly, I am following Singer Willett
(2003), who provide also the S/R code for their examples in the book (e.g.
One way is to use the sprintf-function which returns strings, and then write
these strings in a file. But it is kludgy and there must be better ways...
Best
Søren
-Oprindelig meddelelse-
Fra: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] På vegne af Andrea Toreti
Sendt: 8. maj 2006
I may be out of my statistical depth here, but isn't it the case that if one
has an experimental design with random effects, one has to include the
random effects, even if they appear to be non-significant?
AFAIK there are two reasons: one is the possibility of 'restriction errors'
that arise
Dear R-users
Can anyone inform me of a library or more specifically functions that can
maximise (or calculate) a Pairwsie likelihood from a data.
Better still, i would like to know if there is a function (library) that
fits regression models based on pairwise likelihoods.
Some time ago,I wrote some simple functions to do the
fixed/between/random effects regression.If you want the code,pls drop
me a line.
2006/5/8, David STADELMANN [EMAIL PROTECTED]:
Dear R Users,
Here is another probelm/question.
I would like to run some panel regressions with R. Therefore I
Hello,
I am trying to use the biplot function on an object (called 'acp')
resulting from a prcomp call. However, the scale I get when calling the
biplot function differs according to how I call:
biplot(acp) gives me a scale of -0.10 - 0.10 for the variable loadings
and a scale of -15 - 15 for
Sachin J wrote:
Hi,
How can I get the values of mean and median (not only points but values
too) on the boxplot. I am using boxplot function from graphics package.
Following is my data set
df
[1] 5 1 1 0 0 10 38 47 2 5 0 28 5 8 81 21 12 9 1 12 2 4 22 3
When I run the code : R CMD INSTALL -c SJava_0.69-0.tar.gz
i have:
gzip:stdin:invalid compressed data--crc error
tar:child returned status 1
cannot untar the package
what i do?
[[alternative HTML version deleted]]
__
Hi,
how about round() ?
best
On Monday 08 May 2006 12:36 Søren Højsgaard wrote:
One way is to use the sprintf-function which returns strings, and then
write these strings in a file. But it is kludgy and there must be better
ways... Best
Søren
-Oprindelig meddelelse-
Fra: [EMAIL
Hi all,
I would like to know if there is a way to have the diplay the breaks
of the hist (as myhist$breaks here in the example) as the labels,
instead of the myhist$counts:
input = rpois(1000,5)
myhist = hist(input, breaks = 15, labels=TRUE, axes=FALSE)
So that, instead of having:
94 140 187
try something along these lines:
input - rpois(1000, 5)
myhist - hist(input, breaks = 15, plot = FALSE)
plot(myhist, labels = as.character(myhist$breaks[-1]), axes = FALSE)
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
hello
I built logistic regression model.
To model check I used stepAIC. But I don't know how it
is interpreted . I could not any find any explanation about it
For instance which model is preferable ? What are the critarias
to choose beter model
I will appreciate if you give me an explanation ?
On Mon, 2006-05-08 at 14:22 +0200, Dimitris Rizopoulos wrote:
try something along these lines:
input - rpois(1000, 5)
myhist - hist(input, breaks = 15, plot = FALSE)
plot(myhist, labels = as.character(myhist$breaks[-1]), axes = FALSE)
it seems to give me the labels slided by one, but it
I am hoping for some assistance with a problem that has puzzled me.
Immediately below is the error messages I obtained when I tried to perform
two functions.
(A)
tapply(outcome,na.rm=T,grp,mean)
Error in tapply(outcome, na.rm = T, grp, mean) :
arguments must have same length
On 5/8/2006 4:52 AM, Bob Green wrote:
I am hoping for some assistance with a problem that has puzzled me.
Immediately below is the error messages I obtained when I tried to perform
two functions.
(A)
tapply(outcome,na.rm=T,grp,mean)
Error in tapply(outcome, na.rm = T, grp, mean) :
Do you know the RSQLite package? It uses the DBI package which gives a
common interface to various DB engines.
With it, you can explicitely treat data.frames as tables, and execute SQL
querys on them.
Antonio, Fabio Di Narzo.
2006/5/5, Robert Citek [EMAIL PROTECTED]:
Is there a cheat-sheet
There is a difference in the p- value from 0.000 and 0.012 when I am
using SPSS.
0.000 when I am using the independent variable as scaled 0.012 if I am
using the variable as ordinal.
The independent variable is ordinal but it seems that R is using the
variable as an scaled, because the P-
On Fri, 5 May 2006, Robert Citek wrote:
On May 5, 2006, at 11:30 AM, Thomas Lumley wrote:
In addition to Uwe's message it is worth pointing out that gc()
reports
the maximum memory that your program has used (the rightmost two
columns).
You will probably see that this is large.
Reloading
On Fri, 5 May 2006, Berton Gunter wrote:
?try
as in
result- try (some R expression...)
if (inherits(result,'try-error')) ...do something
else ...do something else
Or tryCatch(), as described in the FAQ
-thomas
Thomas Lumley Assoc. Professor, Biostatistics
On Sat, 6 May 2006, Chad Reyhan Bhatti wrote:
Hello,
My question concerns model selection, stepAIC(), add1(), and coxph().
In Venables and Ripley (3rd Ed) pp389-390 there is an example of using
stepAIC() for the automated selection of a coxph model for VA lung cancer
data.
A statistics
Thomas Kaliwe wrote:
Hi,
I'm sorry for this not being related to R but I think this is a good
place to ask. I'm looking for DOE examples(experiments) that can be done
at home or in class, such as Paper Helicopter, Paper Towel etc.. I'm
thankful for any comment.
Thomas
We at the KULeuven
stepAIC is support software for a book. Please do as the posting guide
asks and consult the book, which explains both AIC and how hte function
works.
The book is even a reference on the help page.
On Mon, 8 May 2006, orkun wrote:
hello
I built logistic regression model.
To model check I
On Mon, 8 May 2006, Thomas Lumley wrote:
On Sat, 6 May 2006, Chad Reyhan Bhatti wrote:
Hello,
My question concerns model selection, stepAIC(), add1(), and coxph().
In Venables and Ripley (3rd Ed) pp389-390 there is an example of using
stepAIC() for the automated selection of a coxph model
Ordinal variables should be stored as ordered factors in R. See ?ordered.
Andy
From: Knut Krueger
There is a difference in the p- value from 0.000 and 0.012
when I am
using SPSS.
0.000 when I am using the independent variable as scaled
0.012 if I am using the variable as ordinal.
I guess that you are on Windows, and that you need to install the
tools required to build R packages. The instructions for doing this
are at
http://www.r-project.org/
and then visiting the Manuals and 'R Installation and
Administration links. See especially section 3 and appendix F. This
can be
I am trying to locate this package but can not find it on the website,
and a search turns up nothing either. Does anybody know where this
package is available?
Also, has anybody heard of a library called CarbamateData? [as in:
library(CarbamateData) ]
Thanks.
Melanie
[[alternative
Hello all,
Is there a way to export C++ unsigned integer to R? (I am trying to parse
files in BPMAP format, and some variables are of type unsigned int (first
declared in C++) ).
I know that to export signed integer to R,
I can do the following in C++:
int Some_INTEGER = 5;
int *p_myInt;
On 5/8/2006 2:44 PM, Boom 2k1 wrote:
Hello all,
Is there a way to export C++ unsigned integer to R? (I am trying to parse
files in BPMAP format, and some variables are of type unsigned int (first
declared in C++) ).
I know that to export signed integer to R,
I can do the following in
It is a trivial matter to adjust the linewidth of the axes and plot
line in a 2D plot:
x - 1:10;
y - 1:10;
par(lwd=2); # axis (and others) linewidth
plot(x,y,l,lwd=6); # plot linewidth
Is there anything comparable to this for persp or wireframe? I've
tried a number of ways, and the
-Original Message-
From: [EMAIL PROTECTED] [mailto:r-help-
[EMAIL PROTECTED] On Behalf Of Melanie Edwards
Sent: Monday, May 08, 2006 10:37 AM
To: r-help@stat.math.ethz.ch
Subject: [R] DRUtils package
I am trying to locate this package but can not find it on the website,
and a
On 5/8/2006 3:02 PM, Matthew Sundling wrote:
It is a trivial matter to adjust the linewidth of the axes and plot
line in a 2D plot:
x - 1:10;
y - 1:10;
par(lwd=2); # axis (and others) linewidth
plot(x,y,l,lwd=6); # plot linewidth
Is there anything comparable to this for persp or
Hello,
i would like to import this txt file:
Giorno;PM10
2006-01-01 10:10;10.3
2006-02-02 20:22;50.3
2006-03-03 23:33;20.1
.
As it's an irregular time series i use zoo as follow:
require(zoo)
z - read.table(c:\\1.csv, sep=;, na.strings=-999, header=TRUE)
q - zoo(z$PM10,
Hi!
I've been trying to compile R with gcc version 3.4.2 on a box running Solaris
2.6.
I keep on hitting this error during the compile:
gcc -I. -I../../src/include -I../../src/include -I/usr/local/include
-DHAVE_CONFIG_H -g -O2 -std=gnu99 -c dynload.c -o dynload.o
In file included from
Dear R-users,
are there pre-programmed packages that implement integer maximization? A
preliminary google search didnt uncover anything that seemed to match.
I want to solve something like Max{U(x) + beta*U(x+y)} with respect to y,
where y can take only integer values.
thanks.
Alnuddy
On Mon, 8 May 2006, Sul, Young L wrote:
Hi!
I've been trying to compile R with gcc version 3.4.2 on a box running Solaris
2.6.
I keep on hitting this error during the compile:
gcc -I. -I../../src/include -I../../src/include -I/usr/local/include
-DHAVE_CONFIG_H -g -O2 -std=gnu99 -c
On Mon, 8 May 2006 22:09:45 +0200 [EMAIL PROTECTED] wrote:
Hello,
i would like to import this txt file:
Giorno;PM10
2006-01-01 10:10;10.3
2006-02-02 20:22;50.3
2006-03-03 23:33;20.1
.
As it's an irregular time series i use zoo as follow:
require(zoo)
z -
Sul, Young L [EMAIL PROTECTED] writes:
Hi!
I've been trying to compile R with gcc version 3.4.2 on a box running Solaris
2.6.
I keep on hitting this error during the compile:
gcc -I. -I../../src/include -I../../src/include -I/usr/local/include
-DHAVE_CONFIG_H -g -O2 -std=gnu99 -c
Perhaps it is just a glitch in the GUI as it does not happen in the
console version.
Calls to demo() without arguments and help(packages=X) return
nothing in the GUI but work in the console. I d/l'd the new 3114-2.3
GUI app (the version just says 3114).
Hello all,
I am trying to create a matrix of 1s and -1s without any repetitions for a
specified number of columns.
e.g. 1s and -1s for 3 columns can be done uniquely in 2^3 ways.
-1 -1 -1
-1 -1 1
-1 1 -1
-1 1 1
1 -1 -1
1 -1 1
1 1 -1
1 1 1
and for 4 columns in 2^4 ways and so on.
I
Dear all,
I have to handle a large matrix (1000 x 10001) where in the
last column i have a value that all the preceding values in the same row
has to be compared to.
I have made the following code :
# generate a (1000 x 10001) matrix, testm
# generate statistics matrix 1000 x 4:
qnt - c(0.01,
Try expand.grid:
expand.grid( c(-1,1), c(-1,1), c(-1,1) )
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
[EMAIL PROTECTED]
(801) 408-8111
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Nameeta Lobo
Sent: Monday,
On Mon, 2006-05-08 at 16:32 -0500, Nameeta Lobo wrote:
Hello all,
I am trying to create a matrix of 1s and -1s without any repetitions for a
specified number of columns.
e.g. 1s and -1s for 3 columns can be done uniquely in 2^3 ways.
-1 -1 -1
-1 -1 1
-1 1 -1
-1 1 1
1 -1 -1
1 -1
Alternate solution with the combinat package:
m - hcube(c(2, 2, 2))
m[m==1] - -1
m[m==2] - 1
Gabor
On Mon, May 08, 2006 at 03:45:37PM -0600, Greg Snow wrote:
Try expand.grid:
expand.grid( c(-1,1), c(-1,1), c(-1,1) )
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
expand.grid(rep(list(c(-1, 1)), 4)) suffices I believe.
-- Bert Gunter
Genentech Non-Clinical Statistics
South San Francisco, CA
-Original Message-
From: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] On Behalf Of Marc
Schwartz (via MN)
Sent: Monday, May 08, 2006 2:50 PM
To:
Bert,
That will result in a data frame, rather than a matrix:
str(expand.grid(rep(list(c(-1, 1)), 4)))
`data.frame': 16 obs. of 4 variables:
$ Var1: num -1 1 -1 1 -1 1 -1 1 -1 1 ...
$ Var2: num -1 -1 1 1 -1 -1 1 1 -1 -1 ...
$ Var3: num -1 -1 -1 -1 1 1 1 1 -1 -1 ...
$ Var4: num -1 -1
That will result in a data frame, rather than a matrix:
Ah, indeed. I suspect that Nameeta doesn't care about the distinction, but
you're certainly right. I think the reason for the extra step was worth
explicitly mentioning, as many who are new to R are hazy about the
distinction (they're
Seth Falcon [EMAIL PROTECTED] wrote in message
news:[EMAIL PROTECTED]
If you want to do more than one thing, you want to create an anonymous
function. Here's an example:
Thanks for the example. Between what you did and some reading that I did I
have got a vague idea of how to use the {}
Sehr geehter
Mein Name ist Jonas Thabo und ich bin der Leiter des Research Department
Committee in der Standard Bank in S|dafrika. Zur Zeit halte ich mich in den
Niederlanden zu einer Fortbildung auf.
Ich kontaktiere Sie bez|glich des Transfers einer sehr gro_en Summe Geldes
vom Konto eines
Is testm really of class matrix? If its a data.frame then manipulation
of matrices is often faster.
On 5/8/06, Monty B. [EMAIL PROTECTED] wrote:
Dear all,
I have to handle a large matrix (1000 x 10001) where in the
last column i have a value that all the preceding values in the same row
has
G’day,
I am trying to read in a table and am getting an error message stating that R
is unable to open a connection to the file.
➢ Avonvegen- read.table(Y:\Study
Sites\Avon\nonspatial\datafiles\archive\Avon_VegEnh.dat , sep=,,
na.string=-, header=TRUE)
Error in file(file, r) : unable to
See 2.16 of the R Windows FAQ.
On 5/8/06, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote:
G'day,
I am trying to read in a table and am getting an error message stating that R
is unable to open a connection to the file.
➢ Avonvegen- read.table(Y:\Study
I guess you use R under windows,then use \\ instead of \.
or use file.choose() to choose the file directly.
2006/5/9, [EMAIL PROTECTED] [EMAIL PROTECTED]:
G'day,
I am trying to read in a table and am getting an error message stating that R
is unable to open a connection to the file.
➢
Hello,
I've been trying to figure out how to access the individual elements
from an object of class phylo.
I am reading in 201 trees created by paup as below.
read.nexus(A_30knj_200t.txt, tree.names= NULL) - anj30
anj30[1]
$tree1
$edge
[,1] [,2]
[1,] -1 -2
[2,] -2 -3
[3,] -3 1
or use forward slashes, exactly the same as in Unix. This has the advantage
of being portable, especially if you make it relative to a starting directory.
read.table(c:/a/b/c.dat)
read.table(c:\\a\\b\\c.dat)
Original message
Date: Tue, 9 May 2006 09:42:03 +0800
From: ronggui [EMAIL
From posts in Sep 2004 and Feb 2005, glm() was raising Error in x[good, ] * w
: non-conformable arrays. I can reproduce this error using:
df1 = data.frame(u=1:10,
v=rpois(10,10),
z=array(1,10,dimnames=list(1:10)))
glm(v~u+offset(log(z)), data=df1,
Dear R Users,
It works fine. while i read the dataset from the text file. (refered R-data
import -export*.pdf)
read.table(c:\\DataES.txt)
V1 V2
1 ACW261 33
2 ACW311 7
3 ACW321 31
4 ACW342 39
5
6
7
Now my question is how do i /*convert text file Data into *.RData* which i
can use
Hi
I have some code up and runnning in Matlab and need it for some stuff better
handled in R. So R.matlab sounded terrific. However, I am having trouble
getting it to go.
I am using R v 2.2.1 (i know i know i need to upgrade), and Windows XP and
Matlab 6.0.0.88 (R12).
here is what I get.
# in
Hi Alex,
author here - thanks for this. Indeed, there is a typo in
MatlabServer.m that gives the error iff the port number is out of
range. So if you set MATLABSERVER_PORT in range it should work.
However, the Matlab code contains a syntax error in that statement
which shouldn't be there. FYI,
Hi All,
I'm trying to do a GARCH simulation in R 2.3.0 release
in Windows XP. I've seen ?garchsim function but that
is
for garch (1,1) and ?garch gives an example for ARCH
simulation. Can anyone help me how can i extend the
help shown in ?garch to GARCH simulation? Please help
me in this regard.
Hi, R users:
I have built an package, and the package depends on two other packages,
I want to install other two packages at the same time when I install my package,
like the following form,
utils:::menuInstallPkgs()
--- Please select a CRAN mirror for use in this session ---
also
80 matches
Mail list logo