Hello,
I have a problem with read.affy.mixed function. I want to read in
together a set of CEL files from chip types Affymettrix HGU133A_2 and
HGU133_Plus_2. I have my files to be read in in one directory together
with a white space delimited file describing them (covdesc). In this
Dear R-users,
I was wondering if it possible to produce clustergrams in R (further
information about clustergrams in the references below).
Schonlau, M. Visualizing Hierarchical and Non-Hierarchical Cluster
Analyses with Clustergrams. Computational Statistics: 2004; 19(1):
95-111.
Schonlau M.
Dear R users,
I want more information about the package Brew.
The example which is given here (
http://learnr.wordpress.com/2009/09/09/brew-creating-repetitive-reports/) is
really complex to understand.
Are there any simple examples through which I can learn about the Brew
package?
Thanks in
Hi,
Did you read the documentation that comes with the package.
require( brew )
?brew
Romain
On 12/08/2009 10:40 AM, Shreyasee wrote:
Dear R users,
I want more information about the package Brew.
The example which is given here (
Hi Romain,
I read that but it seems to be too vague.
Is there any detailed explanation available with simple examples?
Thanks,
Shreyasee
On Tue, Dec 8, 2009 at 5:55 PM, Romain Francois
romain.franc...@dbmail.comwrote:
Hi,
Did you read the documentation that comes with the package.
Hi everybody!
I'm using the mailing list since a few days and I'm starting to wonder
how you guys can know so much about the different functions and packages.
Is there some kind of search tool for functions/packages or did you all
learn through the online help and the mailing list?
Thanks
Hi,
Maybe this would suit you (from the help page)
file.show( system.file(catprint.brew,package=brew) )
file.show( system.file(brew-test-1.brew,package=brew) )
file.show( system.file(brew-test-2.brew,package=brew) )
brew syntax is quite simple, once you read the points 1:5 from the
Hi,
I'm trying to recreate a sensitivity-specificity graph using the
histbackback function. The only problem is that these graphs are typically
drawn with vertical rather than horizontal bar plots (and the histbackback
function only seems to work with horiz=TRUE argument, using horiz=FALSE
Hello Ivan,
there are pretty good manuals for all packages:
http://cran.r-project.org/web/packages/
Johannes
On Tue, Dec 8, 2009 at 11:04 AM, Ivan Calandra ivan.calan...@uni-hamburg.de
wrote:
Hi everybody!
I'm using the mailing list since a few days and I'm starting to wonder how
you guys
http://www.rseek.org/
It is particularly useful to search the mailing list archives of
r-help with rSeek. No matter what kind of problem you have, somebody
has had it before and asked on r-help.
Titus
__
R-help@r-project.org mailing list
Hi,
I am looking for a R package which is capable to process and analysis pictures
of tissues (stained) in an automatic way.
I had a look on biops and EBImage (Bioconductor) but they are not automatic...
Did you already use/know a such package ?
Thanks,
- Martial
Hi,
Depending on what you want to do, the RImageJ package might be for you.
It is a very simple wrapper around ImageJ: http://rsbweb.nih.gov/ij/
Romain
On 12/08/2009 11:40 AM, Martial Sankar wrote:
Hi,
I am looking for a R package which is capable to process and analysis pictures
of
R-help,
I have a list whose elements are data frames.
I want to change the colnames attribute in each element of this list but an
error message
comes up:
lapply(LD_strataNew,function(x) dimnames(x)[[2]][-1]) - as.roman(1:9)[-6]
Error in lapply(LD_strataNew, function(x) dimnames(x)[[2]][-1]) -
You can try this instead:
lapply(LD_strataNew, 'colnames-', as.character(as.roman(1:9)[-6]))
On Tue, Dec 8, 2009 at 9:38 AM, Luis Ridao Cruz lu...@hav.fo wrote:
R-help,
I have a list whose elements are data frames.
I want to change the colnames attribute in each element of this list but an
Using Rscript under Windows, I would like
to use help from the script. But since this
should appear on a window, nothing is done.
I tried the following:
help1 = function (x)
{windows(); options(locatorBell=F)
plot(0,0,axes=F,frame.plot=F)
print(help(x)) # The script may be called by source.
try this:
myDat - read.table(textConnection(group id
1 101
1 201
1 301
2 401
2 501
2 601
3 701
3 801
3 901),header=TRUE)
closeAllConnections()
corr_mat -as.matrix(read.table(textConnection(1 1 .5 0 0 0 0
0 0 0
2 .5 1 0 0 0 0 0 0 0
3 00 1.0 0 0 0 0 0 0
Rolf Turner wrote:
...
To paraphrase Eric Blair: ``SAS-thinkers unbellyfeel R-speak.''
cheers,
Rolf Turner
This *has* to be a fortune.
Pat
##
Attention:\ This e-mail message is privileged and
Hi List.
It may be a very simple question, but I couldnt find an answer on the
internet. Which function (maybe in a specific package) would I use to
perform a Cochran C Test for Homogeneity of Variances?
Stats package have the mantelhaen.test, but Im almost sure that Its not
what I want,
Hello Martial,
Cell segmentation is not easy in many cases. Cells may overlap and
contain complex structures such as organelles or cytoskeletal fibers
that hinder the determination of cell membranes.
There is no one-size-fits-all algorithm that can process and analyse
cellular tissues in an
Hi, I have a data frame and want to merge adjacent rows if some condition is
met. There's an obvious solution using a loop but it is prohibitively slow
because my data frame is large. Is there an efficient canonical solution for
that?
head(d)
rt dur tid mood roi x
55 5523 200 4 subj
On Mon, Dec 7, 2009 at 9:53 AM, Christophe Genolini
cgeno...@u-paris10.frwrote:
Hi all,
In the package rtlu, I use the function savePlot. It is convenient since it
let the user decide in which graphic format he wants his graph to be export.
But when I run R CMD check, I get the following
I have tried the method proposed by Dave, and I must say it works very well.
Not to yield starting estimates for an nls-fit, but as an independent method
for calculating E (which is, by the way, the only paramater that I am
actually interested in).
The calculated values for E (Esp[3]) are on
Code works fine. Thanks!
Muhammad Rahiz | Doctoral Student in Regional Climate Modeling
Climate Research Laboratory, School of Geography the Environment
Oxford University Centre for the Environment, University of Oxford
South Parks Road, Oxford, OX1 3QY, United Kingdom
Tel: +44
Thank you for your reply.
Unfortunately, I don't have any experience with this kind of mathematics. I
don't understand what the wiki page tries to tell me.
Please don't misinterpret his, but I simply don't have the time to learn a)
how to interpret the math description on wikipedia and b) how
Hi Antoon,
now that you mention trying out different methods, maybe you should consider
fitting a
sigmoidal curve to the entire dataset and not only the exponential part (which
constitutes
a very small dataset) as seems to have been the endeavour that initiated the
posting to
R-help.
One
Hello
I have a Linux machine (Ubuntu 8.04 hardy, Gcc version 4.2.4
(i486-linux-gnu) currently running R 2.6.2. I'd like to upgrade to 2.10.
First Question): What is the appropriate way to remove the old version of
R?
Part 2.
After downloading r-base_2.10.0.orig.tar.gz and opening the
Hi Zuzanna --
Zuzanna wrote:
Hello,
I have a problem with read.affy.mixed function. I want to read in
together a set of CEL files from chip types Affymettrix HGU133A_2 and
HGU133_Plus_2. I have my files to be read in in one directory together
with a white space delimited file describing them
On Dec 8, 2009, at 8:29 AM, Geert Janssens wrote:
Thank you for your reply.
Unfortunately, I don't have any experience with this kind of
mathematics. I
don't understand what the wiki page tries to tell me.
Please don't misinterpret his, but I simply don't have the time to
learn a)
how
Dear R users,
I need some advices on the Cox proportional hazard model fitting, as I
don#39;t fully understand the mecanism behind.
The dataset I#39;m working with have individualsbwho can have very long
censored or event time (can be multiple).
For the problem at hand, I#39;m interested
Hi all,
Is there function on R for calculating Modula generators? For example for
primes above 100, e.g 157, i want to know which number generates the group
under multiplication mod 157. i.e i want to find an element whose order is
156. The problem I occur is that modular arithmetic becomes
R-help,
I have a data frame:
mydata
strata nh NhWhfh
1 I 10 26 0.048 0.385
2 II 32 84 0.154 0.381
3 III 16 42 0.077 0.381
4 IV 4 11 0.020 0.364
5 V 10 26 0.048 0.385
7 VII 64 168 0.309 0.381
8VIII 49 129 0.237 0.380
9 IX 22 58 0.107
Try this:
mydata[order(as.roman(as.character(mydata$strata))),]
On Tue, Dec 8, 2009 at 11:59 AM, Luis Ridao Cruz lu...@hav.fo wrote:
R-help,
I have a data frame:
mydata
strata nh Nh Wh fh
1 I 10 26 0.048 0.385
2 II 32 84 0.154 0.381
3 III 16 42 0.077 0.381
4
On Tuesday 8 December 2009, David Winsemius wrote:
On Dec 8, 2009, at 8:29 AM, Geert Janssens wrote:
Thank you for your reply.
Unfortunately, I don't have any experience with this kind of
mathematics. I
don't understand what the wiki page tries to tell me.
Please don't misinterpret
Hi Steve
Have you tried:
apt-cache search gfortran
in a terminal window.
Then
sudo apt-get install theRelevantPackage
I think you also need the Universe repos enabled.
HTH
Iain
--- On Tue, 8/12/09, steve_fried...@nps.gov steve_fried...@nps.gov wrote:
From: steve_fried...@nps.gov
On Dec 8, 2009, at 9:14 AM, Geert Janssens wrote:
On Tuesday 8 December 2009, David Winsemius wrote:
On Dec 8, 2009, at 8:29 AM, Geert Janssens wrote:
Thank you for your reply.
Unfortunately, I don't have any experience with this kind of
mathematics. I don't understand what the wiki page
You need to give your criteria for preferable. For normal-linear
models, REML estimates of variances are unbiased, whereas ML estimates are
downwardly biased. My intuition is that the ML-induced bias would be
worse in small samples. I don't know about other distributions. Likewise I
don't
Hi Luis,
You can't assign to lapply.
--Gray
On Tue, Dec 8, 2009 at 5:38 AM, Luis Ridao Cruz lu...@hav.fo wrote:
R-help,
I have a list whose elements are data frames.
I want to change the colnames attribute in each element of this list but an
error message
comes up:
Dear List,
I need to print out each of 'k' levels of a factor 'n' times each, where
'n' is the number of elements belonging to each factor.
I know that this can normally be done using the gl() command,
but in my case, each level 'k' has an unequal number of elements.
Example with code is as
Sam K upperhalfplane at yahoo.co.uk writes:
Hi all,
Is there function on R for calculating Modula generators? For example for
primes above 100, e.g 157, i want
to know which number generates the group under multiplication mod 157. i.e i
want to find an element whose
order is 156. The
I think I forgot to send the original to the mailing list, so I'm
forwarding it (see below). Sorry about that (and sorry if I did
remember and this is a duplicate). After a few more minutes of
thought, I realized that you should probably make sure that rt, tid,
and mood are also the same in
Dear Sayan,
there is a vcovHC method for panel models doing the White-Arellano covariance
matrix, which is robust vs. heteroskedasticity *and* serial correlation,
although in a different way from that of vcovHAC. You can supply it to coeftest
as well, just as you did. The point is in
I have a set of parameter estimates for a multivariable Cox model predicting
survival duration and a data-frame of new measurements for the variables in
the model, as well as the actual survival duration.
Is there a function to estimate the error the model has on predicting
survival on this new
Hi
A quick question. Standard errors reported by gee/yags differs from the ones in
geeglm (geepack).
require(gee)
require(geepack)
require(yags)
mm - gee(breaks ~ tension, id=wool, data=warpbreaks,
corstr=exchangeable)
mm2 - geeglm(breaks ~ tension, id=wool, data=warpbreaks,
Your question is well taken. I did not give any criteria because I realized
there might be different answers based upon different criteria. Certainly one
fundamental criteria would be that the estimates are BLUE, but this is not the
only criteria one might be used.
John
-Original
jlu...@ria.buffalo.edu wrote:
You need to give your criteria for preferable. For normal-linear
models, REML estimates of variances are unbiased, whereas ML estimates are
downwardly biased.
I suspect that you can't actually say anything general about the
direction of the bias, except for
How about creating an index using multiple columns.
a - with(d, aggregate(dur, list(rt=rt,tid=tid,mood=mood,roi=roi),
sum))
b - with(d, aggregate(x, list(rt=rt,tid=tid,mood=mood,roi=roi),
mean))
c - merge(a, b, by=c(rt,tid,mood, roi))
I suppose one could save some time by not running
I have a set of parameter estimates for a multivariable Cox model predicting
survival duration and a data-frame of new measurements for the variables in
the model, as well as the actual survival duration.
Is there a function to estimate the error the model has on predicting
survival on this new
I think there might be a problem with this approach if roi, tid, rt,
and mood are the same for nonconsecutive rows.
--Gray
On Tue, Dec 8, 2009 at 9:29 AM, Nikhil Kaza nikhil.l...@gmail.com wrote:
How about creating an index using multiple columns.
a - with(d, aggregate(dur,
# BEGIN CODE ##
#!/usr/bin/perl
##
#
# --start, -s = The date you would like to start generating regressors
#--end, -e = When to stop generating holiday regressros
# --scope, -c = D, W for Daily or Weekly respectively (e.g. Does this week
have a particular holiday)
# --file, -f =
On Tue, Dec 8, 2009 at 4:50 PM, Gray Calhoun gray.calh...@gmail.com wrote:
I think there might be a problem with this approach if roi, tid, rt,
and mood are the same for nonconsecutive rows.
True, but I can use the index of my reshape solution. Aggregate was
the crucial ingredient. Thanks
Here is the correct version. The old version is the redirect only version of
the script.
### BEGIN SCRIPT
#!/usr/bin/perl
##
# --start, -s = The date you would like to start generating regressors
#--end, -e = When to stop generating holiday regressros
# --scope, -c = D, W for Daily or
A contrarian point of view:
If you have so little data (relative to the number of parameters to be
estimated, especially NONLINEAR parameters like covariance estimates)that
the ml vs reml bias could be large, then there's so little information
anyway that such bias is the least of your problems
I am having a problem with the censor tick marks with plot and prodlim
I am using the example code from ?prodlim but I have added mark.time=T to
the plot command to get tick marks at censor times.
The problem is the tick marks occur at exactly the same point in each of the
arms.
This seems wrong
Hi Tal!
First let me say that I deeply appreciate the work that you're putting into
this. You're doing good things for our community, and that's great!
I put the planet-R stuff together rather hastily a few years ago, as a way
of seeing whether it was of enough use for the community for
it to
On Dec 8, 2009, at 11:49 AM, bnorth wrote:
I am having a problem with the censor tick marks with plot and prodlim
I am using the example code from ?prodlim but I have added
mark.time=T to
the plot command to get tick marks at censor times.
The problem is the tick marks occur at exactly the
Statistical Programmer - USDA Center for Veterinary Biologics
The USDA Center for Veterinary Biologics is seeking a statistical
programmer with good R skills. The position is a two-year term position in
the CVB Statistics Section with the possibility of becoming permanent.
Applicants must be
I am not sure if its still the case but one of the problems with the
Planet R feed was that it had material in it not related to R or
statistics or any technical subject at all so if you harvest it be
sure to exclude such sources.
On Tue, Dec 8, 2009 at 12:12 PM, Elijah Wright
On Dec 8, 2009, at 12:17 PM, David Winsemius wrote:
On Dec 8, 2009, at 11:49 AM, bnorth wrote:
I am having a problem with the censor tick marks with plot and
prodlim
I am using the example code from ?prodlim but I have added
mark.time=T to
the plot command to get tick marks at censor
Hi list,
This was asked a couple of years ago but I can't find a resolution. Is
there any way to get the coefficients from one of the local polynomial fits
in locfit. I realize that locfit only constructs polynomials at a handful
of intelligently selected points and uses interpolation to
On 8 December 2009 at 12:24, Gabor Grothendieck wrote:
| I am not sure if its still the case but one of the problems with the
| Planet R feed was that it had material in it not related to R or
| statistics or any technical subject at all so if you harvest it be
| sure to exclude such sources.
On Dec 8, 2009, at 12:33 PM, David Grabiner wrote:
Hi list,
This was asked a couple of years ago but I can't find a resolution.
Is
there any way to get the coefficients from one of the local
polynomial fits
in locfit. I realize that locfit only constructs polynomials at a
handful
of
worrying about df (ml vs reml) is just a silly obsession of statisticians (of
which I'm one)
I too have often wondered about the importance of such tertiary issues. My
half-baked understanding is that the main practical difference between ML vs
REML is with regards to ease of computing the
I believe the prediction is done some some sort of grid, then
interpolated to fill in the rest. This is, however, purely for
computational reason, and not for any threoretical reasons. The formal
definition of local polynomials is to do a weighted fit of polynomial at
each point.
Andy
On Tue, Dec 8, 2009 at 2:25 PM, Iain Gallagher
iaingallag...@btopenworld.com wrote:
Hi Steve
Have you tried:
apt-cache search gfortran
in a terminal window.
Then
sudo apt-get install theRelevantPackage
I think you also need the Universe repos enabled.
HTH
Iain
--- On Tue,
Ok, that seems to work. Thanks for the help.
Steve
Steve Friedman Ph. D.
Spatial Statistical Analyst
Everglades and Dry Tortugas National Park
950 N Krome Ave (3rd Floor)
Homestead, Florida 33034
steve_fried...@nps.gov
Office (305) 224 - 4282
Fax (305) 224 - 4147
On Tue, Dec 8, 2009 at 6:24 PM, stephen's mailinglist account
stephen4mailingli...@googlemail.com wrote:
On Tue, Dec 8, 2009 at 2:25 PM, Iain Gallagher
iaingallag...@btopenworld.com wrote:
Hi Steve
Have you tried:
apt-cache search gfortran
in a terminal window.
Then
sudo apt-get
Am doing some vector plots with the arrows() function. Works well. But
what I need to do is supply an arrow for scaling for the reader. I need
to plot an arrow of some known magnitude somewhere on the page
(preferably outside the bounds of the plot, so that it can be seen
clearly) with some
On Tue, 2009-12-08 at 13:42 -0500, Cable, Samuel B Civ USAF AFMC
AFRL/RVBXI wrote:
Am doing some vector plots with the arrows() function. Works well. But
what I need to do is supply an arrow for scaling for the reader. I need
to plot an arrow of some known magnitude somewhere on the page
Hi all,
I used rpart to fit a model, where the covariates are categorical variables.
Then I plotted the tree (mytree) and used the command text to add labels
to the tree.
In the nodes of the tree, the values of the covariates are represented with
a, b or c (tree attached).
Is there a way to show
Thank you for your help!
On Mon, Dec 7, 2009 at 9:53 AM, Viechtbauer Wolfgang (STAT)
wolfgang.viechtba...@stat.unimaas.nl wrote:
If you just want a forest plot, then the forest() function.
If you have the betas and corresponding variances, then you can create a
forest plot with:
I'm doing some data manipulations. I thought originally that I should
use data.frame, as the elements in a data.frame can have multiple
types but the elements in a matrix has to be the same, which all the
data have to convert to strings if there is a single column that is
string.
However, when I
Is it possible?
I was hoping to find something like:
lwd
for the different bars in the barplot but couldn't find it.
Does it exist ?
Thanks,
Tal
Contact
Details:---
Contact me: tal.gal...@gmail.com | 972-52-7275845
Read me:
Hi,
I know there are older threads discussing the quadratcount function in
spatstat. Unfortunately, I could not find a solution to my problem there.
I'm analyzing a point pattern in an irregular polygonal window. Both the
window (an entire country) and the points are projected using WGS84.
On Dec 8, 2009, at 3:26 PM, Tal Galili wrote:
Is it possible?
I was hoping to find something like:
lwd
for the different bars in the barplot but couldn't find it.
Does it exist ?
?box
Thanks,
Tal
--
David Winsemius, MD
Heritage Laboratories
West Hartford, CT
On Dec 8, 2009, at 2:26 PM, Tal Galili wrote:
Is it possible?
I was hoping to find something like:
lwd
for the different bars in the barplot but couldn't find it.
Does it exist ?
Thanks,
Tal
barplot() calls rect() via a wrapper function internally to draw the
rectangles. rect() uses
On Tue, Dec 8, 2009 at 8:52 PM, stephen's mailinglist account
stephen4mailingli...@googlemail.com wrote:
On Tue, Dec 8, 2009 at 7:14 PM, steve_fried...@nps.gov wrote:
Hi Stephen,
After running the script
sudo apt-get update
sudo apt-get install r-base
I launch R and find the it still
Hi all,
I'm a beginner user of R. I am stuck at what I thought was a very obvious
problem, but surprisingly, I havent found any solution on the forum or
online till now.
My problem is simple. I have a file which has entries like the following:
#ID Value1List_of_values
ID1
Two questions:
What is your code?
What do you get with:
options()$dec
decimal_point
.
--
David
On Dec 8, 2009, at 3:55 PM, Gaurav Moghe wrote:
Hi all,
I'm a beginner user of R. I am stuck at what I thought was a very
obvious
problem, but surprisingly, I havent found any
Hi Marc,
You answered my question in depth, leaving me to go with another solution.
Thank you for the detailed answer!
Best,
Tal
Contact
Details:---
Contact me: tal.gal...@gmail.com | 972-52-7275845
Read me: www.talgalili.com
Gaurav -
Here's one way:
x = textConnection('ID1 0.342 0.01,1.2,0,0.323,0.67
+ ID2 0.010 0.987,0.056,1.3,1.5,0.4
+ ID3 0.146 0.1173,0.1494,0.211,0.1257
+
+ ')
y = read.table(x,stringsAsFactors=FALSE)
res =
Sebastian: Can you send me (off-list) your point pattern (ppp_cameroon)
so that I can experiment with it and try to figure out what's
going wrong? (I am one of the maintainers of spatstat.)
Save the point pattern using dput(), e.g.
dput(ppp_cameroon,cameroon.dput)
and then attach
Hi David,
1) My code is as follows:
x=read.table(file,sep=\t)
prate-x[,2]
lrates-(x[,3])
When I do:
print (typeof(lrates)): I get integer
When I do:
for (line1 in lrates) {
lsp-unlist(strsplit(line1,\\,))
}
I get some intermediate
Here are a couple of others to try (using the lattice package):
dotplot(Factor1 ~ Value | Factor2, data=foo, groups=Factor3, auto.key=T)
dotplot(Factor1 ~ jitter(Value) | Factor2, data=foo, groups=Factor3, auto.key=T)
dotplot(Factor3 ~ Value | Factor2*Factor1, data=foo )
dotplot(Factor1:Factor3 ~
-Original Message-
From: r-help-boun...@r-project.org
[mailto:r-help-boun...@r-project.org] On Behalf Of Gaurav Moghe
Sent: Tuesday, December 08, 2009 12:56 PM
To: r-help@r-project.org
Subject: [R] Split comma separated list
Hi all,
I'm a beginner user of R. I am stuck at what
Hi R Users,
I'm trying to re-order the site names (Waseca, Morris, ...). I'm using
following code:
libarry(lattice)
barchart(yield ~ variety | site, data = barley,
groups = year, layout = c(6,1), aspect=.7,
ylab = Barley Yield (bushels/acre),
scales =
On Dec 8, 2009, at 4:12 PM, Gaurav Moghe wrote:
Hi David,
1) My code is as follows:
x=read.table(file,sep=\t)
prate-x[,2]
lrates-(x[,3])
When I do:
print (typeof(lrates)): I get integer
You've already received several solutions from people more R-savvy
than I, so I will try instead
Depending on your sample size, you might be able to just label the nodes
by drawing a random sample from the variable names :-)
Frank
kaida ning wrote:
Hi all,
I used rpart to fit a model, where the covariates are categorical variables.
Then I plotted the tree (mytree) and used the command
On Sun, Dec 06, 2009 at 04:34:18PM -0800, Brock Tibert wrote:
I have successfully created and analyzed my network data. I am
new to R, and Network Analysis too, but I want to color my
vertex based on some of the centrality measures calculated.
Can someone point me in the right direction?
Thanks a lot to everyone who replied. I used some of Stephan Kolassa's
suggestions (thanks, Stephan). Here's my final code which worked (with
comments to explain novices like me who come looking for the same answer)
x=read.table(file,sep=\t,colClasses=character)
prate-x[,2]
lrate-x[,3] #this
On Dec 8, 2009, at 4:42 PM, Gary Miller wrote:
Hi R Users,
I'm trying to re-order the site names (Waseca, Morris, ...).
I'm using
following code:
libarry(lattice)
# slip this code (or one with your preferred ordering) in before the
plot call:
barley$site - factor(barley$site,
I'm having trouble using split on a very large data-set with ~1400 levels of
the factor to be split. Unfortunately, I can't reproduce it with the simple
self-contained example below. As you can see, splitting the artificial
dataframe of size ~13MB results in a split dataframe of ~ 144MB, with an
On Tue, Dec 8, 2009 at 5:19 PM, Nikhil Kaza nikhil.l...@gmail.com wrote:
I suppose that is true, but the example data seem to suggest that it is
sorted by rt.
I was not very clear on that. Sorry.
d$count - 1
a - with(d, aggregate(subset(d, select=c(dur, x, count),
I have xrange which is a range of values from 1 to a max of 162.
I have a yrange of values which really could be any values, but there's a
min and a max. I'd like to create N number of steps between the min and the
max so the length matches the xrange, so that I can plot them together.
Any tips?
On Fri, Dec 4, 2009 at 11:17 PM, Greg Snow greg.s...@imail.org wrote:
I am sure that you mentioned before that your are using 2.7.1, and possibly
even why, but with the number of posts to this list each day and the number
of different posters, I cannot keep track of what version everyone is
On 12/09/2009 05:42 AM, Cable, Samuel B Civ USAF AFMC AFRL/RVBXI wrote:
Am doing some vector plots with the arrows() function. Works well. But
what I need to do is supply an arrow for scaling for the reader. I need
to plot an arrow of some known magnitude somewhere on the page
(preferably
On Fri, Dec 4, 2009 at 6:12 PM, Gavin Simpson gavin.simp...@ucl.ac.uk wrote:
On Fri, 2009-12-04 at 15:18 -0600, Peng Yu wrote:
On Fri, Dec 4, 2009 at 3:06 PM, Peng Yu pengyu...@gmail.com wrote:
On Fri, Dec 4, 2009 at 2:35 PM, Greg Snow greg.s...@imail.org wrote:
The invert argument seems a
I've sent last message only to Titus. Sorry :)
Below my proposition:
Instead of aggregate, try summaryBy from doBy package. It's much
faster. And this package made my life easier :)
try:
summaryBy(dur+x~index, data=d, FUN=c(sum, mean)),
but index should be in data.frame as I remember.
I haven't
All the solutions presented work, but, gosh,
why not fix the source file in the first place?
One way: Open the text file in Excel or a clone thereof, and tell the
importer dialog box to use both tabs and commas as delimiters. Then
save the result, which will have your 'third column' nicely
Tena koe Oliver
?seq
should help if I understand your question correctly
Peter Alspach
-Original Message-
From: r-help-boun...@r-project.org
[mailto:r-help-boun...@r-project.org] On Behalf Of Wells Oliver
Sent: Wednesday, 9 December 2009 11:19 a.m.
To: r-help@r-project.org
Or perhaps just:
plot(xrange, yrange, type =s)
--
On Dec 8, 2009, at 6:00 PM, Peter Alspach wrote:
Tena koe Oliver
?seq
should help if I understand your question correctly
Peter Alspach
-Original Message-
From: r-help-boun...@r-project.org
1 - 100 of 135 matches
Mail list logo