hello,
subset(swiss, Agriculture 60 Examination != c(14,16), select =
c(Agriculture,Examination,Catholic))
Agriculture Examination Catholic
Broye 70.2 16 3.30
Glane 67.8 14 4.20
Aigle 62.0 21 5.16
elyakhlifi mustapha wrote:
hello,
subset(swiss, Agriculture 60 Examination != c(14,16), select =
c(Agriculture,Examination,Catholic))
Try %in% :
subset(swiss, Agriculture 60 Examination %in% c(14,16), select =
c(Agriculture,Examination,Catholic))
Agriculture
Daniel,
Thanks for answering.
your AIC is monotonic increasing
That was the reason I suspected something is wrong, but what?
Either the way I calculate the number of estimiated parameters?
Or the (scale of the) Likelihood itself?
Have you tried a dimensional reduction technique or to
Sorry,
of course it should read
subset(swiss, Agriculture 60 !(Examination %in% c(14,16)), select
= c(Agriculture,Examination,Catholic))
Agriculture Examination Catholic
Aigle 62.0 21 8.52
Avenches60.7 19 4.43
Cossonay
Tony == Tony Chiang [EMAIL PROTECTED]
on Fri, 4 May 2007 00:07:04 +0100 writes:
Tony Hi all,
Tony I am wondering if this is a bug in the Matrix package
Tony or if it something that I am just getting wrong...
Tony here is an example:
[..]
It's a bug.
A
Dear R users,
I tried to fit a cox proportional hazard model to get estimation of stratified
survival probability. my R code is as follows:
cph(Surv(time.sur, status.sur)~ strat(colon[,13])+colon[,18]
+colon[,20]+colon[,9], surv=TRUE)
Error in if (!length(fname) || !any(fname == zname)) { :
For me, the simplest way to find, what is wrong, would be tracing the R code:
library(debug)
mtrace(cph)
cph(Surv(time.sur, status.sur)~ strat(colon[,13])+colon[,18]
+colon[,20]+colon[,9], surv=TRUE)
... then find the place of the error and analyze how to adjust the function
call arguments to
Here is some information on this regression in R
http://tolstoy.newcastle.edu.au/R/help/06/02/21153.html
Abdus Sattar wrote:
I am want to use tobit regression for left censored panel/longitudinal
data. Could you please provide me the name of library and/or package
that will give me option
I am trying to use the survfit() function with the newdata argument to
produce predicted survivor curves for a particular covariate profile.
The main purpose of the plot will be to visualise the effect of snp1,
coded 0 and 1. In my Cox model I have stratified by one variable, edu, and
so I
Vladimir Eremeev wrote:
For me, the simplest way to find, what is wrong, would be tracing the R
code:
library(debug)
mtrace(cph)
cph(Surv(time.sur, status.sur)~ strat(colon[,13])+colon[,18]
+colon[,20]+colon[,9], surv=TRUE)
... then find the place of the error and analyze how to
Hi, how can i use data to forecast next time period value, if data has
been influenced by a change in legislation?
thank you.
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
I suggest reading the chapter on debugging in 'Writing R Extensions' and
using the tools described there (such as options(error=recover)) would be
at least as effective. Even traceback() would have helped the readers.
Note that there is no function 'cph' in R: presumably we are supposed to
T. Kounouni wrote:
Hi, how can i use data to forecast next time period value, if data
has been influenced by a change in legislation? thank you.
Well, you could use chicken entrails.
cheers,
Rolf Turner
Hello,
I have got two matrices with different length. The matrices have 3
columuns. The first two are coordinates. The third is a measurement.
Now I want to get a subtraction between every single value of the
third column (between matrix1 and matrix2), but only if the two
first coordinates in
The site is up again. This was a wrong manipulation. Sorry.
..°}))
) ) ) ) )
( ( ( ( (Prof. Philippe Grosjean
) ) ) ) )
( ( ( ( (Numerical Ecology of Aquatic Systems
) ) ) ) ) Mons-Hainaut University, Belgium
( ( ( ( (
Hello,
it might be a very simple question but I cannot find the solution (I tried a ||
b, any(a,b)... but none works). My problem is:
I have two vectors,
a - c(TRUE,FALSE,FALSE)
b - c(TRUE,FALSE,TRUE)
and I would like to obtain a vector that indicates if it is TRUE in any of the
two vectors.
I think I found the solution (just after sending the email).
The following apparentely works:
a == TRUE | b ==TRUE
Regards,
Federico
Federico Abascal [EMAIL PROTECTED] escribió: Hello,
it might be a very simple question but I cannot find the solution (I tried a ||
b, any(a,b)... but none
try this:
ind1 - do.call(paste, c(as.data.frame(mat1[, 1:2]), sep = \r))
ind2 - do.call(paste, c(as.data.frame(mat2[, 1:2]), sep = \r))
mat1[ind1 %in% ind2, 3] - mat2[ind2 %in% ind1, 3]
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of
Or you can use the operator:
a - c(TRUE,FALSE,FALSE)
b - c(TRUE,FALSE,TRUE)
a b
[1] TRUE FALSE FALSE
On 5/4/07, Federico Abascal [EMAIL PROTECTED] wrote:
Hello,
it might be a very simple question but I cannot find the solution (I tried a
|| b, any(a,b)... but none works). My problem
Meant the | operator
a | b
[1] TRUE FALSE TRUE
On 5/4/07, jim holtman [EMAIL PROTECTED] wrote:
Or you can use the operator:
a - c(TRUE,FALSE,FALSE)
b - c(TRUE,FALSE,TRUE)
a b
[1] TRUE FALSE FALSE
On 5/4/07, Federico Abascal [EMAIL PROTECTED] wrote:
Hello,
it might be a
you need:
a | b
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuven, Belgium
Tel: +32/(0)16/336899
Fax: +32/(0)16/337015
Web: http://med.kuleuven.be/biostat/
Hello,
The answer is in one page of the R Wiki (just created to address such a
question, by the way):
http://wiki.r-project.org/rwiki/doku.php?id=tips:data-manip:calc_on_two_tables
..°}))
) ) ) ) )
( ( ( ( (Prof. Philippe Grosjean
) )
Dear Sattar,
You can use the survreg function in the survival package, which is part of
the standard R distribution. [BTW, help.search(Tobit) would have led you
to that.] There are other possibilities as well -- e.g., the tobit function
in the VGAM package.
I hope this helps,
John
hello,
how can I do to drop decimal after the comma please for example for tthis line
print(P)
[1] 62.00 1.00 7.661290 5.20 17.10 2.318801
how canI do to keep only 62 1 7.66 5.2 17.12.32
thanks
__
ble contre les
I have just learned how to play with map, but something weird
(or not) is happening.
Suppose I want to draw a map of two countries (that have disconnected
components), like Argentina and Brazil.
If I command:
library(maps)
library(mapdata)
map(worldHires, c(Argentina, Brazil))
It works fine.
Hallo,
just look ?round. It should help.
Corinna
-Ursprüngliche Nachricht-
Von: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED] Im Auftrag von elyakhlifi
mustapha
Gesendet: Freitag, 4. Mai 2007 14:59
An: R-help@stat.math.ethz.ch
Betreff: [R] decimal values
hello,
how can I do to drop
OK, so I just want to go on record as noting that following the
instructions exactly (exactly, exactly, exactly, emphatically) on the
web given by Duncan, Gabor, and Tony at the page below, and following
the instructions in Writing R Extensions results in a successful Windows
build of a package.
Thanks very much for your help.
I remain puzzled by the original behavior of reshape though...
Does anyone have an explanation to this ?
Regards,
David
-Message d'origine-
De : Gabor Grothendieck [mailto:[EMAIL PROTECTED]
Envoyé : jeudi 3 mai 2007 19:27
À : GOUACHE David
Cc :
hi,
hi, good morning everyone. I have a time series with binary outputs like :
000100100.etc. Now I want to forecast the future values of
that. Can anyone please tell me whether there is any tools exist in literature
for dealing with this kind of binary observation? If
Can anybody explain me how do i get Correlation Coefficient R² in a
non-linear regresion analisys performed with nls()?. Thanks in advance.
--
Adrián J. Montero Calvo
[EMAIL PROTECTED]
__
R-help@stat.math.ethz.ch mailing list
Might there be an (semi-)automated procedure to create a minimal,
personal package, for my eyes only, that I can load with a
libray(MyStuff) command? This would be preferable to having to
source() the files. Is there already such a procedure?
Joe
-Original Message-
From: [EMAIL
Regarding your comments on changing paths,
1. Rcmd.bat in the batchfiles distribution
http://code.google.com/p/batchfiles/
can be used to avoid having to change the path (other than path
changes required for perl and tex which you would probably want
anyways to be able to use those). You
hongyuan cao wrote:
Dear R users,
I tried to fit a cox proportional hazard model to get estimation of
stratified survival probability. my R code is as follows:
cph(Surv(time.sur, status.sur)~ strat(colon[,13])+colon[,18]
+colon[,20]+colon[,9], surv=TRUE)
Error in if (!length(fname) ||
Hi Gabor:
I tried the link below, but it seems to be broken.
-Original Message-
From: Gabor Grothendieck [mailto:[EMAIL PROTECTED]
Sent: Friday, May 04, 2007 10:05 AM
To: Doran, Harold
Cc: Duncan Murdoch; r-help@stat.math.ethz.ch
Subject: [SPAM] - Re: [R] R package development
After cleaning up your workspace to keep only the objects that you will need
to have in your package, you can save it as MyStuff.RData
save(list=ls(), file=file.path(my.package.dir, MyStuff.RData))
To use this package in future sessions all you need is
attach(file.path(my.package.dir,
I have a function definition such as
f - function (pattern = .*\\.txt) {}
in the manpage this has to be documented as
f - function (pattern = .*.txt)
in order to get the correct display (with double backslash) in the R console
when issuing `?f', but this causes complains from `R CMD
I was able to click the link in your reply and it worked so
there is some problem with your email viewer. You can google for:
code google batchfiles
and use the first hit or type the URL in by hand.
On 5/4/07, Doran, Harold [EMAIL PROTECTED] wrote:
Hi Gabor:
I tried the link below, but
On Fri, 4 May 2007, Alberto Monteiro wrote:
I have just learned how to play with map, but something weird
(or not) is happening.
Suppose I want to draw a map of two countries (that have disconnected
components), like Argentina and Brazil.
If I command:
library(maps)
library(mapdata)
hi! All:
I can do up or down arrow on keyboard to browse through command history
on R console in windows. However, I can't do that on a linux xterm or
console... I wonder how to make this feature work on linux... it will
make working a lot more efficient...
Thanks
Hao
I guess it depends on what you want to be able to do with such a private
package; e.g., does it not need to have any documentation (i.e., the Rd
files)? If all you want is to be able to access the objects, you can
just save() all those objects (mostly functions, I presume) in a .rda
file, and
On 5/4/07, Adrian J. Montero Calvo [EMAIL PROTECTED] wrote:
Can anybody explain me how do i get Correlation Coefficient R² in a
non-linear regresion analisys performed with nls()?. Thanks in advance.
It may seem obvious how to define the multiple correlation
coefficient R^2 for a non-linear
On 5/4/2007 9:30 AM, Doran, Harold wrote:
OK, so I just want to go on record as noting that following the
instructions exactly (exactly, exactly, exactly, emphatically) on the
web given by Duncan, Gabor, and Tony at the page below, and following
the instructions in Writing R Extensions results
On Fri, 2007-05-04 at 10:48 -0400, Hao Liu wrote:
hi! All:
I can do up or down arrow on keyboard to browse through command history
on R console in windows. However, I can't do that on a linux xterm or
console... I wonder how to make this feature work on linux... it will
make working a
Hi Megh,
One good reference for your problem could be
Harvey(1989), Forecasting, Structural Time Series Models and the Kalman
Filter, Cambridge University Press.
In section 6.6.2, he includes the treatment of binomial and multinomial
observations using state space models.
Hope it helps,
I have unzipped the R-2.5.0.tar.gz
gzip -dc R-x.y.z.tar.gz | tar xvf -
2. then #./configure
3. ./configure
checking build system type... x86_64-unknown-linux-gnu
checking host system type... x86_64-unknown-linux-gnu
loading site script './config.site'
loading build specific script './config.site'
We are announcing a short course, which is scheduled to take place in
June 2007 at The University of Reading, UK. Summary information is given
below.
For more detailed information and registration forms please see
http://www.ssc.rdg.ac.uk providing your address and/or fax number, or
email
Every things that moves, moves because something have moved, Thomas d'Aquin
said. Explain cleary in detail you problem with more details.
Justin BEM
Elève Ingénieur Statisticien Economiste
BP 294 Yaoundé.
Tél (00237)9597295.
- Message d'origine
De : [EMAIL PROTECTED] [EMAIL PROTECTED]
On 5/3/07, ivo welch [EMAIL PROTECTED] wrote:
dear R experts:
sorry, I have to ask this again. I know that the answer is in section
7.2 of S Programming, but I don't have the book (and I plan to buy
the next edition---which I hope will be titled S/R programming ;-) ).
I believe the
?format
elyakhlifi mustapha wrote:
hello,
how can I do to drop decimal after the comma please for example for tthis line
print(P)
[1] 62.00 1.00 7.661290 5.20 17.10 2.318801
how canI do to keep only 62 1 7.66 5.2 17.12.32
thanks
Megh == Megh Dal [EMAIL PROTECTED]
on Fri, 4 May 2007 00:12:25 -0700 (PDT) writes:
Megh hi, good morning everyone. I have a time series with
Megh binary outputs like :
Megh 000100100.etc. Now I want to
Megh forecast the future values of that. Can anyone
On 5/4/07, Duncan Murdoch [EMAIL PROTECTED] wrote:
On 5/4/2007 9:30 AM, Doran, Harold wrote:
OK, so I just want to go on record as noting that following the
instructions exactly (exactly, exactly, exactly, emphatically) on the
web given by Duncan, Gabor, and Tony at the page below, and
On Fri, 2007-05-04 at 10:06 -0500, Pramod Anugu wrote:
I have unzipped the R-2.5.0.tar.gz
gzip -dc R-x.y.z.tar.gz | tar xvf -
2. then #./configure
3. ./configure
checking build system type... x86_64-unknown-linux-gnu
checking host system type... x86_64-unknown-linux-gnu
loading site script
Given the following, one of the things I am trying to see is what % of
draws are below a certain number:
lambda - 3
rate - 5
n - 5
set.seed(123)
v - replicate(n, rexp(rpois(1,lambda), rate))
vv - unlist(v)
cat(% of draws below 0.1:, round(length(subset(vv, vv
0.1))/length(vv)*100,0), %\n)
In
On 5/4/2007 11:31 AM, Gabor Grothendieck wrote:
On 5/4/07, Duncan Murdoch [EMAIL PROTECTED] wrote:
On 5/4/2007 9:30 AM, Doran, Harold wrote:
OK, so I just want to go on record as noting that following the
instructions exactly (exactly, exactly, exactly, emphatically) on the
web given by
We can also say that in recent versions of R you have to choose
deliberately to disable readline when building R. (This is made quite
hard: I am currently building R on a brand new Solaris 10 system, and I
had to work to get the right 64-bit libreadline linked in. Just having
I had a question about Random Forests. I have a text file with 10
dichotomous variables and a bivariate response vector. I read this file
into R as a data frame, and then used the command randomForest(Response ~.,
dataset, etc.. where Response is the column header of the response
variable and
On 5/4/07, Duncan Murdoch [EMAIL PROTECTED] wrote:
On 5/4/2007 11:31 AM, Gabor Grothendieck wrote:
On 5/4/07, Duncan Murdoch [EMAIL PROTECTED] wrote:
On 5/4/2007 9:30 AM, Doran, Harold wrote:
OK, so I just want to go on record as noting that following the
instructions exactly (exactly,
First, your example works happily in a 1Gb Windows machine under R 2.5.0:
your R (2.2.1) is well overdue for an update.
What you are generating is a random number of rexp(rate=10) random
variables. So all you need is
n - 100
N - sum (rpois(n, 26))
vv - rexp(N, 10)
mean(vv[vv 0.1])
which
You can try a probit on panel data or see it as survival models if you want to
analyse the time spent before the appearance of the event (code 1). a good
reference is Econometric Analysis of Cross Section and Panel Data, 2nd Edition
Jeffrey M. Wooldridge
Justin BEM
Elève Ingénieur
The best, of course, would be to get rid of Perl altogether.
In Python, it is possible to make standalone executables. Is it possible
to also do this in Perl, then one could eliminate a perl install. Or, is
it possible to use Python to accomplish what perl is currently doing? I
may be getting
Hi,
I like to know the simple way to list the R package names in our linux
system.
Thanks
Ying
[[alternative HTML version deleted]]
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
Jacques Wagnor said the following on 5/4/2007 8:53 AM:
Given the following, one of the things I am trying to see is what % of
draws are below a certain number:
lambda - 3
rate - 5
n - 5
set.seed(123)
v - replicate(n, rexp(rpois(1,lambda), rate))
vv - unlist(v)
cat(% of draws below
On Fri, 2007-05-04 at 13:06 -0400, Hu, Ying (NIH/NCI) [E] wrote:
Hi,
I like to know the simple way to list the R package names in our linux
system.
Thanks
Ying
See ?installed.packages
HTH,
Marc Schwartz
__
R-help@stat.math.ethz.ch
Just googling I found this:
http://www.perlmonks.org/?node_id=186402
On 5/4/07, Doran, Harold [EMAIL PROTECTED] wrote:
The best, of course, would be to get rid of Perl altogether.
In Python, it is possible to make standalone executables. Is it possible
to also do this in Perl, then one
On Fri, 2007-05-04 at 12:05 -0500, Bill Vorias wrote:
I had a question about Random Forests. I have a text file with 10
dichotomous variables and a bivariate response vector. I read this file
into R as a data frame, and then used the command randomForest(Response ~.,
dataset, etc.. where
Bill,
A couple more points:
1. Please use an informative subject line. I'd deleted the original
post w/o reading if I didn't catch Marc's reply.
2. Are you sure you have bivariate response? To me bivariate means
two variables, and randomForest surely does not handle that (at least
I'm trying to install rgl 0.71 on a redhat enterprise 4, x86_64. I have
tried using R 2.2.1, 2.3.1, and 2.5.0. I have successfully installed
this version of rgl, using R 2.2.1 on an rhel4 i386 host. On the x86_64
host, I receive the following configuration error:
checking GL/gl.h
See:
packs - .packages(all=T)
length(packs)
--
Henrique Dallazuanna
Curitiba-Paraná-Brasil
25° 25' 40 S 49° 16' 22
Ohttp://maps.google.com/maps?f=qhl=enq=Curitiba,+Brazillayer=ie=UTF8z=18ll=-25.448315,-49.276916spn=0.002054,0.005407t=kom=1
On 04/05/07, Hu, Ying (NIH/NCI) [E] [EMAIL PROTECTED]
Hi,
I'd like to partition a 2d probability density function into regions of
equal probability. It is straightforward in the 1d case, like
qnorm(seq(0,1,length=5)) but for 2d I'd need more constraints.
Any suggestions for how to approach this? Is seems like a spatial
sampling problem but I'm
One can also use:
library()
Henrique Dallazuanna wrote:
See:
packs - .packages(all=T)
length(packs)
__
R-help@stat.math.ethz.ch mailing list
hi,
is there a way to save a R object into workspace instead of into a
file during a running of function?
thanks,
--
Weiwei Shi, Ph.D
Research Scientist
GeneGO, Inc.
Did you always know?
No, I did not. But I believed...
---Matrix III
__
hi doug: yikes. could I have done better? Oh dear. I tried to make
my example clearer half-way through, but made it worse. I meant
set.seed(1);
fe = as.factor( as.integer( runif(100)*10 ) ); y=rnorm(100); x=rnorm(100);
print(summary(lm( y ~ x + fe)))
deleted
Coefficients:
I have used the pp/par combination for Perl before. It is pretty straight
forward to convert an existing perl script into a stand alone windows
executable.
Both the Activestate licence and the Perl Artistic licence allow for embedding
a script and perl interpreter together and distributing
On Fri, May 04, 2007 at 03:45:10PM -0400, Weiwei Shi wrote:
hi,
is there a way to save a R object into workspace instead of into a
file during a running of function?
if I understand the question correctly you want the 'super-assignment'
operator `-' as in
-cut---
R
The problem with the size could be handled by making each
perl script a separate function and combing them into one
large script so that the overhead only gets incurred once.
On 5/4/07, Greg Snow [EMAIL PROTECTED] wrote:
I have used the pp/par combination for Perl before. It is pretty straight
sorry for my English, staying in US over 7 years makes me think I were
ok, now :(
anyway, here is an example and why I need that:
suppose I have a function like the following:
f1 - function(){
line1 - f2() # assume this line takes very long time, like more than 30 min
# then I need to save
On 04/05/2007 4:25 PM, Greg Snow wrote:
I have used the pp/par combination for Perl before. It is pretty straight
forward to convert an existing perl script into a stand alone windows
executable.
Both the Activestate licence and the Perl Artistic licence allow for
embedding a script
On 04/05/2007 3:43 PM, Smith, Troy (NIH/NCI) [C] wrote:
I'm trying to install rgl 0.71 on a redhat enterprise 4, x86_64. I have
tried using R 2.2.1, 2.3.1, and 2.5.0. I have successfully installed
this version of rgl, using R 2.2.1 on an rhel4 i386 host. On the x86_64
host, I receive the
Hello,
I need to import a data set.
I have never imported data files with R. I have always worked on simulated data.
I have looked at R Data Import/Export manual.
It is a bit peculiar because my data base is already an R object called japan.
I guess it is not yet a data set, and I don't know
[for those that worry about these things, this _is_ a homework
assignment. However, it's not an R homework, it's a Geography
and History homework... and I want to use R to create a pretty
map]
Roger Bivand wrote:
Is there any way to associate one color to each country?
Try:
map_poly_obj -
It certainly would be excellent if installing perl could be eliminated.
One additional thing that I really dislike about the R installation is that
one needs find on one's path and that conflicts with find on Windows
so other applications unrelated to R that use scripts can suddenly break
because
Hi David,
This is an interesting question! I was wondering is the density function like?
What sort of partitions are you looking for?
Sorry I am not of much help, but this question perhaps better-defined?
Best,
Ranjan
On Fri, 4 May 2007 14:48:36 -0500 (CDT) David Forrest [EMAIL PROTECTED]
Dear R Developers:
What is the proper log(likelihood) for 'factanal'? I believe it
should be something like the following:
(-n/2)*(k*(log(2*pi)+1)+log(det(S)))
or without k*(log(2*pi)-1):
(-n/2)*log(det(S)),
where n = the number of (multivariate) observations.
On 04/05/2007 9:00 PM, Alberto Vieira Ferreira Monteiro wrote:
[for those that worry about these things, this _is_ a homework
assignment. However, it's not an R homework, it's a Geography
and History homework... and I want to use R to create a pretty
map]
Roger Bivand wrote:
Is there any
--- [EMAIL PROTECTED] wrote:
Hello,
I need to import a data set.
I have never imported data files with R. I have
always worked on simulated data.
I have looked at R Data Import/Export manual.
It is a bit peculiar because my data base is already
an R object called japan.
Then you
Does anyone know if there is a way to specify regressors dynamically
rather than explicitly?
More specifically, I have a data set in long format that details a
number of individuals and their responses to a question (which can be
positive, negative, or no answer). Each individual answers as many
86 matches
Mail list logo