Dear all,
I have a matrix X with 47 lines and say 500 columns - values are in {0,1}.
I'd like to compare lines.
For that, I first did:
for (i in 1:(dim(X)[1]-1))
for (j in (i+1):dim(X)[1]) {
Y - X[i,]+Y[j,]
etc.
but, since it takes a long time, I would prefer avoding loops;
for
Dear Francisco,
Many thanks for your attention to my contribute. As
regard its translation into English -also other people
asked me-, maybe in next months it could take place.
When it will be available I inform the R-help mailing
list.
Best regards,
Vito
--- F Z [EMAIL PROTECTED] ha scritto:
Cuichang Zhao wrote:
Hello,
how can use change the plot function to change the range of axises. I want my graph from a certain range [a, b], instead of from the min to max of of datas?
if i want draw a line instead of dots, should i use both plot and lines function.
for example:
plot(x, y);
Try this,
plot(x,y,xlim=c(a,b),ylim=c(c,d))
lines(z,q) will only add lines onto an existing plot, not create a new plot!
Have a look at:
http://cran.r-project.org/doc/manuals/R-intro.html#Graphics
Cheers, Henrik
Cuichang Zhao wrote:
Hello,
how can use change the plot function to change the range
Dear all,
I am expecting a Poisson error distribution in my lme with
weights=varFunc().
The weigths= varPower (form= fitted (.)) doesn´t work due to missing
values in the response:
Problem in lme.formula(fixed = sqrt(nrmainaxes + 0...: Maximum number of iterations reached without convergence.
Dear Avneet
If you fit a multiple linear regression model your independent
variables will not be orthogonal and therefore it is difficult
to divide the r square on the different variables with a
meaningful interpretation.
In a balanced analysis of variance design the situation is
somehow
On Wed, Jan 26, 2005 at 01:39:36AM +0100, Peter Dalgaard wrote:
Thomas Lumley [EMAIL PROTECTED] writes:
On Tue, 25 Jan 2005, Florian Menzel wrote:
Hello all,
I found a weird result of the GLM function that seems
to be a bug.
No, the problem is that you are using the Wald test
Göran Broström [EMAIL PROTECTED] writes:
Try it with
b=c(rep(0,8),rmultinom(1,24,rep(1/8,8)))
This is the conditional distribution, given the sufficient statistic.
Yes (and that was of course the point). I was surprised that there is
a difference at all (between this and rep(3,8) and
Hello,
I would like to be able to compose an expression using some of the arguments I
give to a function, that is assuming I have a list such as:
my.list-list(A=1:3,B=c(Brave, new, world))
and a function
printComponent-function(component){
if (component==A||component==B){
Alex Sanchez wrote:
Hello,
I would like to be able to compose an expression using some of the arguments I
give to a function, that is assuming I have a list such as:
my.list-list(A=1:3,B=c(Brave, new, world))
and a function
printComponent-function(component){
if (component==A||component==B){
Hi Alex,
you could use something like this:
printComponent - function(component, lis){
if(component %in% names(lis)) return(lis[component])
else cat(component, in not an element of, substitute(lis), \n)
}
##
my.list - list(A=1:3, B=c(Brave, new, world))
printComponent(A, my.list)
Hello,
has anyone of you ever implemented (or integrated a C-SourceCode of) a
Multiobjective Evolutionary Algorithm in R, especially the NSGA-II
developed by Kalyanmoy Deb?
Thank you very much for your help,
Heike Trautmann
Dr. Heike Trautmann
On 25-Jan-05 Weiguang Shi wrote:
Sorry.
c = sum(1.0 / pow((double) j, alpha)), j=1,2,...,N
Understood! Though I think you mean
1/c = sum(1.0 / pow((double) j, alpha)), j=1,2,...,N
Anyway, provided alpha 1 the series converges for the
infinite sum, but if you limit the range to 1...N
Hi
If you go to
http://www.r-project.org/
there is search section, where you can search not only mailing lists
but also R site.
My favorite site is Paul Johnson's Rtips (try to search it by Google)
Cheers
Petr
On 25 Jan 2005 at 16:16, Cuichang Zhao wrote:
hello,
i wonder what command
Hi
first of all do
?plot
?lines
and read the help page (especially about parameters xlim and ylim)
and/or try to use examples provided.
plot(1:100, rnorm(100), xlim=c(1,50), type=l)
Cheers
Petr
On 25 Jan 2005 at 21:34, Cuichang Zhao wrote:
Hello,
how can use change the plot function
When you start you R session, see if it says
[Previously saved workspace restored]
immediately before the first command prompt.
I think the only way a workspace can be restored automatically is if you
used save.image() at some point, in which case it would save all objects
as a hidden file
MM == Martin Maechler [EMAIL PROTECTED]
on Tue, 25 Jan 2005 09:59:03 +0100 writes:
Paul == Paul Murrell [EMAIL PROTECTED]
on Tue, 25 Jan 2005 13:40:15 +1300 writes:
Paul Hi
Paul Cari G Kaufman wrote:
Hello,
Does anyone know how to make movies in R by making a
Please give a simple example of the input data and output that you
desire. It is difficult to understand from you partial codes what you
mean. For example what is Y ?
Are you trying to find add values from pairs of rows ? If so, please see
my posting pairwise difference operator where I wanted to
Hi All,
I was updating my help pages for a package when all of a sudden I got the
following error when using R CMD INSTALL:
Building/Updating help pages for package 'depmix'
Formats: text html latex example
Error in load(zfile, envir = envir) : error reading from connection
Removing all
On Wed, Jan 26, 2005 at 11:01:03AM +0100, Peter Dalgaard wrote:
Göran Broström [EMAIL PROTECTED] writes:
Try it with
b=c(rep(0,8),rmultinom(1,24,rep(1/8,8)))
This is the conditional distribution, given the sufficient statistic.
Yes (and that was of course the point). I was
Hi All,
I was updating my help pages for a package when all of a sudden I got the
following error when using R CMD INSTALL:
Building/Updating help pages for package 'depmix'
Formats: text html latex example
Error in load(zfile, envir = envir) : error reading from connection
Removing all
Göran Broström [EMAIL PROTECTED] writes:
On Wed, Jan 26, 2005 at 11:01:03AM +0100, Peter Dalgaard wrote:
Göran Broström [EMAIL PROTECTED] writes:
Try it with
b=c(rep(0,8),rmultinom(1,24,rep(1/8,8)))
This is the conditional distribution, given the sufficient statistic.
It is part of a function to determine Dice's index in the framewok of AFLP
analysis.
X is a binary matrix which value for each strain (lines) and each base pair
(columns) is 1 where there is a peak and 0 where there is no peak as
biologists explained to me.
The first step is to compare each
On Wed, Jan 26, 2005 at 02:03:25PM +0100, Peter Dalgaard wrote:
Göran Broström [EMAIL PROTECTED] writes:
On Wed, Jan 26, 2005 at 11:01:03AM +0100, Peter Dalgaard wrote:
Göran Broström [EMAIL PROTECTED] writes:
Try it with
b=c(rep(0,8),rmultinom(1,24,rep(1/8,8)))
See if this does what you want:
m - matrix(round(runif(24)), 4, 6) # simulate some data
m
[,1] [,2] [,3] [,4] [,5] [,6]
[1,]010101
[2,]011100
[3,]111000
[4,]100010
library(gtools) # Install
Good afternoon,
do you know if R provides a special package to make networks ?
Thank you,
Angela
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
Slight edit?
It is rumored that on Wed, 26 Jan 2005 09:24:44 -0500
Liaw, Andy [EMAIL PROTECTED] wrote:
See if this does what you want:
m - matrix(round(runif(24)), 4, 6) # simulate some data
m
[,1] [,2] [,3] [,4] [,5] [,6]
[1,]010101
[2,]0
Hi, there:
I am wondering if there is a package in R (doing decison trees) which
can provide some methods to evaluate the significance of attributes. I
remembered randomForest gives some output like that. Unfortunately my
current computing env. cannot handle my datasets if I use
randomForest. So,
Dear R users:
I am looking for the source code for the R function extractAIC. Type the
function name doesn't help:
extractAIC
function (fit, scale, k = 2, ...)
UseMethod(extractAIC)
environment: namespace:stats
And when I search it in the R source code, the best I can find is in (R
source
Angela:
Also, have you done a search at www.r-project.org - search -
R site search? A search for networks there just now produced 218
hits.
Have you reviewed the list of officially contributed packages is
available at www.r-project.org - CRAN - (select a local mirror like
From: Robert Kruus
Slight edit?
Yes. Thank you. Jacques caught that as well...
Andy
It is rumored that on Wed, 26 Jan 2005 09:24:44 -0500
Liaw, Andy [EMAIL PROTECTED] wrote:
See if this does what you want:
m - matrix(round(runif(24)), 4, 6) # simulate some data
Isn't this in the FAQ?
Try:
methods(extractAIC)
[1] extractAIC.aov* extractAIC.coxph* extractAIC.glm*
extractAIC.lm*
[5] extractAIC.negbin* extractAIC.survreg*
Non-visible functions are asterisked
Then use getAnywhere() or getS3method() to get the particular one you're
Dear Paul,
Thank you for your help!
1. I know how to use expression(italic(whatever)), but I don't know
how to make expression(whatever) to produce italic font.
2. With regard to the CM font, I guess psfrag may be to best way to get
around the problem. Indeed, I have been using psfrag since I
apply() statements **are** disguised loops and therefore are **not**
necessarily more efficient than explicit looping. Their principal advantage
is usually code readability.
As another readability issue, note that x[[i]][[j]][[k]] can be abbreviated
to x[[c(i,j,k]].
I leave to others the
I'm using R 2.0.1 in windows XP (and am not currently subscribed to this
mailing list).
I have a USGS dataset, a text file with fixed width fields, that includes
dates as 6-digit integers in the form yrmoda. I could either read them that
way, or with yr, mo, and da as separate integers. In
Alexandre Sanchez Pla asanchez at ub.edu writes:
:
: Hi,
:
: I am working with lists whose terms are lists whose terms are lists.
Although
: the real ones contain locuslink identifiers and GO annotations (I work with
the
: Bioconductor GO) package, I have prepared an simplified example of
Actually, what you want is sapply.
sapply(tst.list, [[, VAL)
Kevin
Alexandre Sanchez Pla wrote:
Hi,
I am working with lists whose terms are lists whose terms are lists. Although
the real ones contain locuslink identifiers and GO annotations (I work with the
Bioconductor GO) package, I have
Message: 63
Date: Wed, 26 Jan 2005 04:28:51 + (UTC)
From: Gabor Grothendieck [EMAIL PROTECTED]
Subject: Re: [R] chron: parsing dates into a data frame using a
forloop
To: r-help@stat.math.ethz.ch
Message-ID: [EMAIL PROTECTED]
Content-Type: text/plain; charset=us-ascii
Benjamin M.
Hi all --
I'm trying to use R for my Analysis of Categorical Data class, but I
can't figure out how to do a weighted linear trend analysis. I have a table
of categorical data, to which I've assigned weights to the rows and columns.
I need to calculate r and M^2, which is apparently done in
I'm using R 2.0.1 in a class I teach, with most students working under
windows XP. We have a data frame with the first column containing the
factor site, and five water-quality variables at each site. As one part
of exploring these data, I'd like the students to run
by(ourdata,site,pairs).
Benjamin M. Osborne Benjamin.Osborne at uvm.edu writes:
:
: Message: 63
: Date: Wed, 26 Jan 2005 04:28:51 + (UTC)
: From: Gabor Grothendieck ggrothendieck at myway.com
: Subject: Re: [R] chron: parsing dates into a data frame using a
: forloop
: To: r-help at stat.math.ethz.ch
:
Hi
Martin Maechler wrote:
MM == Martin Maechler [EMAIL PROTECTED]
on Tue, 25 Jan 2005 09:59:03 +0100 writes:
Paul == Paul Murrell [EMAIL PROTECTED]
on Tue, 25 Jan 2005 13:40:15 +1300 writes:
Paul Hi
Paul Cari G Kaufman wrote:
Hello,
Does anyone know how to make
On Wed, 26 Jan 2005, Jean Eid wrote:
Hi,
I am using the NLSY79 data (longitudinal data from the Bureau of labour
stats in the US). The extractor exctracts this data in a wide format and
I need to reshape it into a long format.
What I am doing right now is to do it in chuncks for each and evry
Hi,
I am using the Survival package, more precisely the cox.zph function, to plot
log(Hazard rate) over time.
if I type
plot(temp.zph[2])
then I get the plot I want. However, I want to change the label of the y axis
that cox.zph prints.
plot(temp.zph[2],xlab='Days', ylab='log hazard for
Dear Thomas,
It is very interesting to read FAQ 7.22. I don't see any example of
lattice/trellis graphics linked with print() statement. How should I specify
print()? Here are the actual codes I used to create the .ps file:
x=y=myaa
grid=expand.grid(x=x,y=y)
Here's a different suggestion. Create a bunch of image files, and then
use an image browser (GQview is one of the best; if you're on Win look
at ACDSee) to view them as a slide show. Good image browsers read
images in advance and should not produce flickering. I haven't
experimented though with
You can turn on plot history from the menu in the windows() device. Then
whatever plots you make, you can page-up/page-down to navigate through them.
(This is specific to Windows.)
Andy
From: David Parkhurst
I'm using R 2.0.1 in a class I teach, with most students
working under
windows
thank you Thomas and Chuck for the helpfull codes. My problem was with the
times variable adn chuck's example made that really clear. I just have one
more question regarding reshape. when varying are not the same length
reshape complains, so I generated the name of the variables that are
missing
Thank you very much again. Your code is really helpful
and it demontrates the power of R to a new-comer like
me.
--- [EMAIL PROTECTED] wrote:
Understood! Though I think you mean
1/c = sum(1.0 / pow((double) j, alpha)),
j=1,2,...,N
You are right.
Weiguang
Chris,
While there might be certain functions in R that I'm not aware
of for such purposes, anything (well, almost) can be done in R
if you're willing to write a line or two of code, and that's
the beauty of R. In a course I recently taught, I used the
following code to generate individual data
I want to specify two factors and their interaction as random effects using
the function lme(). This works okay when I specify these terms using the
function Error() within the function aov(), but I can't get the same model
fitted using lme(). The code below illustrates the problem.
a -
If you read the Help file for lme (!), you'll see that ~1|a*b is certainly
incorrect.
Briefly, the issue has been discussed before on this list: the current
version of lme() follows the original Laird/Ware formulation for **nested**
random effects. Specifying **crossed** random effects is
Hi list,
I am just a new user of R.
How can I run stepwise regression in R?
Is there a graphic user interphase for any of the spatial packages inculded
in R, such as gstat, geoR and someothers. I am mainly interested interactive
variogram modelling and mapping.
Thanks
Mahdi
--
Thanks again Andy.
The definition of AC is understood, yet I have trouble
picturing the amount of clear clustering structure
it measures. To put things into perspective, for two
series
1,2,1000,1001
and
1,2,3,1000
agnes(x, method=single) generates ac values of
0.998998 and 0.0.7492477
-Original Message-
From: Tim F Liao [mailto:[EMAIL PROTECTED]
In a course I recently taught, I used the
following code to generate individual data from grouped data,
which would give the same results as using fweight=count in Stata.
FWIW, I wrote a little function to extract variable importance as defined in
the CART book a while ago. It's rather limited: Only works for regression
problem, and you need to set maxsurrogate=0 and maxcompete=0. It may (or
may not) help you:
varimp.rpart - function(x) {
dev - x$frame[,
Dear R-help:
Is there a way to extract a robust covariance matrix from optim? I am looking
for one of the form H V H; where H is the inverse of the Hessian, and V is
dLL/dT * dLL/dT (LL=Log-likelihood, T=vector of parameters, d is the partial).
Optim returns H, which approximates for the
Hi,
I am looking for a package to do the clustering analysis using the
expectation maximization algorithm.
Thanks in advance.
Ming
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting
Dear R people,
I would like to know if it is possible to compute the relative
importance of the inputs to a neurol network computed with the function
nnet in the nnet package.
Thank you so much!!!
Victor Robles
[[alternative HTML version deleted]]
-Original Message-
From: Weiguang Shi
Thanks again Andy.
The definition of AC is understood, yet I have trouble
picturing the amount of clear clustering structure
it measures. To put things into perspective, for two
series
1,2,1000,1001
and
1,2,3,1000
agnes(x,
Well I am not sure that can call a single figure a cluster. Sure it's not near
the others but how can you conceptually measure it's cluster properties. It
seems reasonable that there has to be some form of doubt about it.
Back to that Google search hit number 3 www.stat.ncu.edu.tw/teacher/
Greetings -
Is it possible, inside .First.lib,
to find out the version number of the package that is being loaded?
If only one version of the package has been installed,
we could scan the DESCRIPTION file, something like
.First.lib - function(lib, pkg) {
library.dynam(spatstat, pkg, lib)
Is it possible, inside .First.lib, to find out the
version number of the package that is being loaded?
Berwin Turlach kindly informs me that installed.packages() will extract
details of all installed packages.
So if only one version of package spatstat has been installed,
we can
Dear all,
I am beginner using R. I have a question about it. When you use it,
since it is written by so many authors, how do you know that the
results are trustable?(I don't want to affend anyone, also I trust
people). But I think this should be a question.
Thanks,
Ming
Berton Gunter a écrit :
If you read the Help file for lme (!), you'll see that ~1|a*b is certainly
incorrect.
Briefly, the issue has been discussed before on this list: the current
version of lme() follows the original Laird/Ware formulation for **nested**
random effects. Specifying **crossed**
Berton Gunter a écrit :
If you read the Help file for lme (!), you'll see that ~1|a*b is certainly
incorrect.
Briefly, the issue has been discussed before on this list: the current
version of lme() follows the original Laird/Ware formulation for **nested**
random effects. Specifying **crossed**
How do you know that any results from any software package are
trustable? I'm not sure that the number of authors has anything to do
with it.
If you are extremely paranoid, you can reprogram everything you do a
few times in a large number of completely different languages written
by different
What makes you trust any software?
There are some obvious points. First of all the code is open so if you know
enough you can actually read the code and make sure it does what you want.
Secondly you can replicate a process using two pieces of software and compare
the results. You can check
68 matches
Mail list logo
