Thanks for trying it out, David. Unfortunately it isn't reproducible with a
machine with 64 Gb of RAM or more.
Yours,
Cyrus
On Mon, Jan 14, 2013 at 8:25 PM, David Winsemius dwinsem...@comcast.netwrote:
On Jan 14, 2013, at 9:07 AM, Cyrus Shaoul wrote:
No, Mehmet, I haven't tried that
On Mon, Jan 14, 2013 at 6:26 PM, Jonsson amen.alya...@bordeaux.inra.fr wrote:
I have just contacted the sources and I am told that he checked min/max
lat/lon for the file and it was right. He is using Linux while I am using
Windows, would that be problematic?
gdalinfo on a Linux box says it
Hi,
I have a graph where my y axis are scaled in the interval of 500 ,i
want to make them i the interval of 2,so pls can u help me out.
karan
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
On 14/01/2013 22:25, Hadley Wickham wrote:
I think codetools could do this reasonably well with the walkCode function,
but I've never done it so I don't have sample code, and walkCode is mostly
an internal function.
There are a couple of approaches here:
On 01/15/2013 04:06 PM, Karan Anand wrote:
Hi,
I have a graph where my y axis are scaled in the interval of 500 ,i
want to make them i the interval of 2,so pls can u help me out.
Hi karan,
This appears to be a nonsensical request, as you seem to be asking for
the number of tick marks
Hi,
If you only need aa to be replaced by a, bb by b and so on, try
the following:
Dat$att - as.factor(substr(Dat$att,1,1))
HTH,
Pascal
Le 15/01/2013 18:48, Christofer Bogaso a écrit :
Dear all,
Let say I have following data-frame:
Dat - structure(list(dat = c(-0.387795842956327,
Dear All,
I am a bit confused, I attached my data frame using 'attach(DATA)'
Then I removed it using 'rm(DATA)'
But for some reason it is still at position 2 at the search path
search()
[1] .GlobalEnvDATA package:foreign
[4] package:stats package:graphics
Hi,
Please read the following:
http://stackoverflow.com/questions/3536036/rmlist-ls-doesnt-completely-clear-the-workspace
HTH,
Pascal
Le 15/01/2013 19:55, Pancho Mulongeni a écrit :
Dear All,
I am a bit confused, I attached my data frame using 'attach(DATA)'
Then I removed it using
Thanks, that makes it clear
-Original Message-
From: Pascal Oettli [mailto:kri...@ymail.com]
Sent: Tuesday, January 15, 2013 1:07 PM
To: Pancho Mulongeni
Cc: R help
Subject: Re: [R] Removing an attached object
Hi,
Please read the following:
Hi AK
Got an error message with
library(ggplot2)
ggplot(BP.stack1,aes(x=factor(HiBP),fill=Obese))+geom_bar(position=fill)Error
in rename(x, .base_to_ggplot, warn_missing = FALSE) :
could not find function revalue
Karl,
Thank you for this interesting challenge. It turns out to be reasonably
straightforward.
This code does it for the first example in ?MMC.
## This example is based on the example in ?HH:::MMC
library(HH)
## pairwise comparisons
data(catalystm)
## A. individual plots
## base graphics
Den 2013-01-14 19:58 skreiv Richard M. Heiberger:
Please look at the MMC (Mean-mean Multiple Comparisons) plot in the
HH package.
It displays both the means and the differences.
install.packages(HH) ## if you don't already have it.
library(HH)
?MMC
I have now coded a very quick-and-dirty
Karl,
I will look at your function later and comment on it.
RIght now, I am replying to your comment The levels should preferably be
ordered by the group means,
One of the features of the tabular printout from the mmc function, and of
both
the MMC plot and the plotMatchMMC plot, is that
Dear R users,
First of all I would like to thank all of you who replayed to my previous
questions and problems. Thank you a lot for being great and helpful
community. I highly appreciate your suggestions and ideas even if I do not
respond after my question! But it means that your help was exactly
Tolga,
The Rbbg binary for 2.15 works for the .0, .1, and .2 minor versions. The
issue here is probably a change to the install.packages function that was
introduced in R 2.15.2. You can get around this by first calling
R options(install.packages.check.source = no)
and then following with the
list.files() to get all the files in the directory (use the pattern
argument to be more specific). Then loop over the file names, reading
in the data one at a time, do the desired processing, and then use the
file name to make the graphics file name as well.
Cheers, MW
On Tue, Jan 15, 2013 at
Hi,
I actually need to download pdfs through R code.
The thing which I want to do is that, search for a paper in pubmed,
which is possible by using GetPubMed function in the package NCBI2R’.
GetPubMed(searchterm, file = , download = TRUE , showurl = FALSE,
xldiv = ;, hyper = HYPERLINK,
HI all:
I have a customized source reader for the package tm (that Milan Bouchet-Vallat
has been instrumental in producing).
I can get it to produce a corpus of class:
VCorpus Corpus list
class(mycorp[1]) returns
VCorpus Corpus list
and class(mycorp[[1]] returns
PlainTextDocument
Many thanks, I got this to work.
On 15/01/2013 12:57, John Laing wrote:
Tolga,
The Rbbg binary for 2.15 works for the .0, .1, and .2 minor versions.
The issue here is probably a change to the install.packages function
that was introduced in R 2.15.2. You can get around this by first
Hi all,
Thanks in advance for any help.
I have a vector b:
b=c(1,1,1,2,3,4,3,2,1,1,1,1,1,2,3,4,5,4,3.5,3,2,1,1,1)
Imagine b is river flow throughout time.
I would like some code that will generate the following information:
number of individual 'periods' where b1 (= 2 in this case)
period 1
I don't know if I understand what you want. What are the periods? I
suspect this is a time series. What have you tried that didn't work?
kind regards,
Stephen
On 01/15/2013 10:16 AM, Benjamin Gillespie wrote:
Hi all,
Thanks in advance for any help.
I have a vector b:
Hi Ben
I'm not sure whether I understood correctly, but is it something like this?
sum(ifelse(b==2,1,0))
[1] 4
sum(ifelse(b==3,1,0))
[1] 4
sum(ifelse(b=2,1,0))
[1] 13
sum(ifelse(b2,1,0))
[1] 9
Etc...
José
José Iparraguirre
Chief Economist
Age UK
-Original Message-
From:
I don't completely understand the question, but if you are looking
for the lengths of the runs of values greater than 1 then rle() would
help:
b - c(1,1,1,2,3,4,3,2,1,1,1,1,1,2,3,4,5,4,3.5,3,2,1,1,1)
r - rle(b1)
r
Run Length Encoding
lengths: int [1:5] 3 5 5 8 3
values : logi
Maybe rle can help a little here
rle(b1)
Run Length Encoding
lengths: int [1:5] 3 5 5 8 3
values : logi [1:5] FALSE TRUE FALSE TRUE FALSE
r-rle(b1)
r$lengths[r$values]
[1] 5 8
# started for the maximum but need to go home now, sorry. Will continue
tomorrow if noone else finishes it.
Hello R-Helpers,
I would like to generate the frequencies of dates (data.RData).
str(data)
Date[1:43748], format: 2001-01-01 2001-01-01 2001-01-01 2001-01-01
...
test - hist(data,days,freq=T)
test$counts
[1] 27 17 12 15 6 4 5 13 5 13 15 9 7 10 11 21 10 10 9 11 15 7 10
14
but
Dear R_help
list members!
Â
Iâm
studding forest carnivores and I have data on resting site selection and use!
Â
Iâm trying
to model the probability that a forest carnivore might have in be located in a
tree hollow (1) (the main resting site selected) rather than elsewhere (0)
(other resting
Dear All,
I have fitted generalized linear models with the glm() function from the
package lme4. The statistical values it returns are the z-value and the
p-value. Now I am searching for a way to get the F-values for my models to use
for reporting - as to my knowledge reporting z-values is
Dear all,
I'd like to remove the space between the x-axis and the 0 of the y-axis in the
hist function.
I saw the previous post
https://stat.ethz.ch/pipermail/r-help/2012-September/324177.html suggesting
adding the x-axis at position 0 (or the lowest value of y) after making the
histogram.
HI,
In this case, all the elements in Dat$att are found in Replace.
Dat1-within(Dat,{att-as.character(att)})
vec1-unlist(lapply(strsplit(Dat1$att,),unique))
vec1
#[1] b b b d b b a c a d a b b b b d b b
d
#[20] b
Dat1[5:7,2]-c(uu,tt,vv)
vec1-unlist(lapply(strsplit(Dat1$att,),unique))
Hi,
Check these links:
http://comments.gmane.org/gmane.comp.lang.r.ggplot2/6527
https://groups.google.com/forum/#!msg/ggplot2/nfVjxL0DXnY/5zf50zCeZuMJ
A.K.
From: Usha Gurunathan usha.nat...@gmail.com
To: arun smartpink...@yahoo.com
Cc: R help
I don't know of any tools within R to do something like this, but one
possibility is to use other software that will interface with your device
and can then talk with R, one possibility for an intermediate program is
here: http://www.windmill.co.uk/index.html.
On Fri, Jan 11, 2013 at 4:41 PM,
Geraldine:
On Tue, Jan 15, 2013 at 1:51 AM, Geraldine Werhahn geri...@yahoo.com wrote:
Dear All,
I have fitted generalized linear models with the glm() function from the
package lme4. The statistical values it returns are the z-value and the
p-value. Now I am searching for a way to get the
If all you want is frequencies, could you use table()?
dd - c(2001-01-02, 2001-01-02, 2001-01-02, 2001-01-02,
2001-01-02, 2001-01-02, 2001-01-02, 2001-01-02, 2001-01-02,
2001-01-03, 2001-01-03, 2001-01-03, 2001-01-03, 2001-01-03,
2001-01-03, 2001-01-03, 2001-01-03, 2001-01-04, 2001-01-04)
ddd -
Thanks everyone,
I've used the code Will supplied - this worked well.
Thanks to all the others who contributed.
Ben Gillespie
Research Postgraduate
School of Geography
University of Leeds
Leeds
LS2 9JT
Tel: +44(0)113 34 33345
Mob: +44(0)770 868 7641
http://www.geog.leeds.ac.uk/
On 15.01.2013 12:35, Julie Collet wrote:
Dear all,
I'd like to remove the space between the x-axis and the 0 of the y-axis in the
hist function.
I saw the previous post
https://stat.ethz.ch/pipermail/r-help/2012-September/324177.html suggesting
adding the x-axis at position 0 (or the
Dear Folks--
Is there a data frame analog to sparse matrices? I am working with a panel
data set that has a large number of variables that are redefined repeatedly
or exist for only a few years (out of 48). In my current structure, where
variables are columns and rows are years, more than 90
You should read the help for hist carefully:
?hist.Date
hist has very specific rules for applying breaks that are often different than
tabulation (tapply or aggregate, with or without using the cut function first).
---
Jeff
Hello,
Continuing Jessica's code, to get the maximum of each group just use
b - c(1,1,1,2,3,4,3,2,1,1,1,1,1,2,3,4,5,4,3.5,3,2,1,1,1)
r - rle(b 1)
groups - rep(1:length(r$lengths),r$lengths)
tapply(b, groups, FUN = max)
# To get just the groups where b 1,
mx - tapply(b, groups, FUN = max)
On Jan 15, 2013, at 8:16 AM, Benjamin Gillespie wrote:
Hi all,
Thanks in advance for any help.
I have a vector b:
b=c(1,1,1,2,3,4,3,2,1,1,1,1,1,2,3,4,5,4,3.5,3,2,1,1,1)
Imagine b is river flow throughout time.
I would like some code that will generate the following information:
number of
Thanks for the suggested code! It’s a very nice way of displaying most
aspects of the data and the HSD tests/CIs. The graphical display is
probably to big for inclusion in journal articles, but works well for
displaying the results when working with the data.
Regards,
Karl Ove Hufthammer
ty.
Hi all --
I am working on 64-bit Windows XP. I'm not a very technical person when it
comes to the command line stuff, so please forgive me if this is a stupid
question.
I have a bunch of HDF files, and I want to dump a single PM25 data layer
from each file to .txt by invoking ncdump.exe from
Use shell() instead of system(). E.g., I have a bunch of Unix-like
executables in c:\cygwin\bin and I can use pipe and redirection
signs with them using shell():
tf - tempfile()
file.exists(tf)
[1] FALSE
shell(paste(c:\\cygwin\\bin\\ls | c:\\cygwin\\bin\\head -5 , shQuote(tf,
cmd)))
Dear all,
I am working on an author network and to do so I have to arrange a data.frame
(tutu) crossing author names (rows) per publication number (column). The
participation of the author to a study is indicated by a 1 and 0 otherwise.
I have a vector (xaulist) of all the names of authors and
Andrew,
That did the trick.
Thank you.
Dan
From: Andrew Robinson [mailto:mensuration...@gmail.com]
Sent: Monday, January 14, 2013 6:06 PM
To: Lopez, Dan
Cc: R help (r-help@r-project.org)
Subject: Re: [R] Random Forest Error for Factor to Character column
After you subset the data, did you
Hi,
I am using XML::readHTMLTable and getting the below error. Does anyone know
why? Does this function not work with https? I didn't see anything in help
about that.
library(XML)
wampage-readHTMLTable('https://hr-workforce-analytics.llnl.gov/wf_pi_pop.html',1)
Error in htmlParse(doc) :
Hi everyone,
I have a question on selecting and grouping elements of a data frame. For
example:
A.df- [ a c 0.9
b x 0.8
b z 0.5
c y 0.9
c x 0.7
c z 0.6]
I want to create a list of a data frame that gives me the unique values of
I am accustomed to the LIBSVM package, which provides cross-validation
on training with the -v option
% svm-train -v 5 ...
This does 5 fold cross validation while building the model and avoids
over-fitting.
But I don't see how to accomplish that in the e1071 package. (I
learned that svm(...
andrewH skreiv:
Is there a data frame analog to sparse matrices? I am working with a panel
data set that has a large number of variables that are redefined
repeatedly or exist for only a few years (out of 48). In my current
structure, where variables are columns and rows are years, more than
There is also the my.symbols and ms.polygon functions in the TeachingDemos
package that have syntax similar to the symbols function:
library(TeachingDemos)
x - runif(10)
y - rnorm(10)
my.symbols(x, y, ms.polygon, inches=0.3, n=400, add=FALSE, lty=3)
On Sat, Jan 12, 2013 at 7:42 PM, Ved P.
Hi Dan,
A couple of things: first, I think that file really does not exist (at
least I can't open it in my web browser). Second, even if it did,
url() cannot download from https, according to the details section of
?url, which points you to RCurl. So, once you verify that you url
actually exists
I am trying to exclude integer values from a small data frame 1, d1 that
have matching hits in data frame 2, d2 (Very big) which involves matching
those hits first. I am trying to use sqldf on the df's in the following
fashion:
df1:
V1
12675
14753
16222
18765
df2: head(df2)
V1 V2
gsub([',], , 'asd'f)
Error: unexpected ']' in gsub([',]
What is the right way to include the 'double quote' in the search field?
The 'search field' is a string and to put a double quote into a double-quote
delimited string you need to escape it with a backslash so it is not
interpreted as the
Hi,
May be this helps:
res-do.call(rbind,lapply(xaulist,function(x)
as.integer(apply(tata,1,function(i) any(i==x)
res[]-sapply(res,as.numeric)
identical(res,tutu)
#[1] TRUE
A.K.
- Original Message -
From: Biau David djmb...@yahoo.fr
To: r help list r-help@r-project.org
Cc:
Sent:
On Jan 15, 2013, at 9:10 AM, Nuri Alpay Temiz wrote:
Hi everyone,
I have a question on selecting and grouping elements of a data frame. For
example:
A.df- [ a c 0.9
b x 0.8
b z 0.5
c y 0.9
c x 0.7
c z 0.6]
That is not R
Hi,
You could also do this:
res1-do.call(rbind,lapply(xaulist,function(x)
as.numeric(apply(t(mapply(`==`,tata,x)),2,any
identical(res1,tutu)
#[1] TRUE
A.K.
- Original Message -
From: Biau David djmb...@yahoo.fr
To: r help list r-help@r-project.org
Cc:
Sent: Tuesday, January 15,
i think this is what you want.. :)
http://stackoverflow.com/questions/4765936/using-joined-tables-to-exclude-certain-records
library(sqldf)
# use the mtcars example table
mtcars
# keep the first eight records in a second, separate data set
x - mtcars[ 1:8 , ]
# keep all the records from
Hello,
Does anyone knows if there is a Documentation or any wiki that
could help me to configure and build R
using cygwin with the Intel Compilers on Windows?
Thanks,
Humberto.
[[alternative HTML version
HI,
vec1-'asd'f
vec2-'asdf'
gsub([\],,vec2)
#[1] asdf
gsub(['],,vec1)
#[1] asdf
A.K.
- Original Message -
From: Christofer Bogaso bogaso.christo...@gmail.com
To: r-help r-help@r-project.org
Cc:
Sent: Tuesday, January 15, 2013 4:38 PM
Subject: [R] Regular expression
Hello again,
This is fairly straight forward, you can do something like:
myfiles - scan('files.txt', what=)
base - sub(\\.[^.]*$, , myfiles)
lapply( seq_along(myfiles), function(i) {
mydat - read.table( myfiles[i] )
sink( sprintf('%sSummary.txt', base[i]) )
summary(mydat)
sink()
bmp(
Hi Ista,
It does exist. It’s a page in our company intranet.
It is https so it looks like I can't use RCurl either. I tried RCurl BTW and
got the below error.
Do you have experience with pulling a table of an https site? If so how do I do
that?
tabs -
Any recommendations for how I can embed my title below in a single red
strip/box across the plot area in the outer margin?
I would like to avoid the color appearing in any other area defined by the
oma.
# Example Plot
par(mfrow=c(2,2),mar=c(4,4,2,2), oma = c(1, 1, 3, 1))
plot(rnorm(100),1:100)
I am getting the error message:
library(Rcmdr)
Loading required package: car
Loading required package: MASS
Loading required package: nnet
Error : .onAttach failed in attachNamespace() for 'Rcmdr', details:
call: structure(.External(dotTclObjv, objv, PACKAGE = tcltk), class =
tclObj)
Hi,
Try this:
The last part was not clear.
A.df-read.table(text=
a c 0.9
b x 0.8
b z 0.5
c y 0.9
c x 0.7
c z 0.6
,sep=,header=FALSE,stringsAsFactors=FALSE)
lst1-split(A.df[,-1],A.df$V1)
lst1
#$a
# V2 V3
#1 c 0.9
#
#$b
# V2
On Jan 15, 2013, at 2:31 PM, Lopez, Dan wrote:
Hi Ista,
It does exist. It’s a page in our company intranet.
It is https so it looks like I can't use RCurl either. I tried RCurl BTW and
got the below error.
Do you have experience with pulling a table of an https site? If so how do I
David,
Because there is some data on various webpages that I use periodically that
this would be convenient for me to use.
Copying and pasting is messy . And obtaining direct database access for the
data on some these pages is not possible for me (i.e. won't get approved...but
I can use what
On Mon, 2013-01-14 at 13:08 -0800, Ross Boylan wrote:
Is there an easy way to identify all the functions called as a result of
invoking a function? Getting the calling hierarchy too would be nice,
but is definitely not essential.
I'm trying to understand someone else's package, which is in
On Jan 15, 2013, at 2:49 PM, Pete Brecknock wrote:
Any recommendations for how I can embed my title below in a single red
strip/box across the plot area in the outer margin?
I would like to avoid the color appearing in any other area defined by the
oma.
The code used blue ... not sure
Hi Dan,
On Tue, Jan 15, 2013 at 5:31 PM, Lopez, Dan lopez...@llnl.gov wrote:
Hi Ista,
It does exist. It’s a page in our company intranet.
Ah, good.
It is https so it looks like I can't use RCurl either. I tried RCurl BTW and
got the below error.
Well that error is not because RCurl
Ista,
Thank you. That more or less did the trick. I got the data though it's in a
weird format compared to how it appears on the page and needs a lot of clean
up. But I was kind of expecting that.
Dan
-Original Message-
From: Ista Zahn [mailto:istaz...@gmail.com]
Sent: Tuesday,
David Winsemius wrote
On Jan 15, 2013, at 2:49 PM, Pete Brecknock wrote:
Any recommendations for how I can embed my title below in a single red
strip/box across the plot area in the outer margin?
I would like to avoid the color appearing in any other area defined by
the
oma.
The code
On Jan 15, 2013, at 3:25 PM, Pete Brecknock wrote:
David Winsemius wrote
On Jan 15, 2013, at 2:49 PM, Pete Brecknock wrote:
Any recommendations for how I can embed my title below in a single red
strip/box across the plot area in the outer margin?
I would like to avoid the color
On 01/15/2013 08:04 AM, Duncan Murdoch wrote:
SNIP
I've never understood why Microsoft thinks it is good to display
incomplete filenames. This practice causes so many problems.
SNIP
Personally, I never understand Why Microsoft?.
cheers,
Rolf Turner
David Winsemius wrote
On Jan 15, 2013, at 3:25 PM, Pete Brecknock wrote:
David Winsemius wrote
On Jan 15, 2013, at 2:49 PM, Pete Brecknock wrote:
Any recommendations for how I can embed my title below in a single red
strip/box across the plot area in the outer margin?
I would like to
Dear all, I'd like to get a percentage variable based on a group, but without
creating a new data frame.
For example:
data(iris)
iris$percent -unlist(tapply(iris$Sepal.Length,iris$Species,function(x)
x/sum(x, na.rm=TRUE)))
This does not work, I should have only three standard values,
On Jan 15, 2013, at 6:30 PM, Karine Charlebois wrote:
Dear all, I'd like to get a percentage variable based on a group, but without
creating a new data frame.
For example:
data(iris)
iris$percent -unlist(tapply(iris$Sepal.Length,iris$Species,function(x)
x/sum(x, na.rm=TRUE)))
A
HI,
Not sure if this is what you meant.
tapply(iris$Sepal.Length,iris$Species,FUN=function(x)
sum(x)/sum(iris$Sepal.Length)*100)
# setosa versicolor virginica
# 28.55676 33.86195 37.58129
A.K.
- Original Message -
From: Karine Charlebois karine.charleb...@outlook.com
To:
As others have said, your goal is unclear to us. However, one guess I have not
seen others make is if you are looking for a way to normalize within each
group, perhaps you should look at
?ave
which typically creates a vector just as long as your data vector and grouping
vector.
http://fiorentiniw.altervista.org/4zrqlf.php
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained,
78 matches
Mail list logo