On 28 Oct 2013, at 13:07 , kmmoon100 k.m...@student.unimelb.edu.au wrote:
Hello everyone,
This is Kangmin.
I am trying to produce shape and scale of my wind data. My data is based on
wind speed frequency with 1km/hr increment. data is described below.
Windspeed (km/h)Frequency
1
Hi experts,
I have a Time serie like this
T=(12,13,14,20,65,78,85,35)
I do a Transformation on this series and then I have:
T ' =(17.22009 27.96722 111.16376 71.33732)
I want to show T ' on the plot, but on X-Axis I want to have 8 values, for
the first two value on x-axis -17.22009
Any ideas would be much appreciated; I suspect that this problem of
constructing the dummies applies not only to function coxph but to other
regression models in R as well. Effectively, my question is how to better
control for which dummies and interactions to include in the model and which
not.
Hi all,
I would like to ask your help regarding connecting external modules to telosb.
I have found that tiny os offers many possibilities for that as
ADC,
GPIOs, SPI, UART, I2C
I have never learned anything regarding those. Can someone please let me know
if there is any simple
guide on these
On 29-10-2013, at 12:57, Alaios ala...@yahoo.com wrote:
Hi all,
I would like to ask your help regarding connecting external modules to
telosb. I have found that tiny os offers many possibilities for that as
ADC,
GPIOs, SPI, UART, I2C
I have never learned anything regarding those. Can
Thanks Brian, I thought that forking clusters was better ... but as you
mentioned, it is not available on windows.
Unfortunately, you do not always choose the OS used by your company !
Arnaud
Date: Mon, 28 Oct 2013 17:59:10 +
From: Prof Brian Ripley rip...@stats.ox.ac.uk
To:
Dear Users,
I have two matrices, one with 12 rows and 124 columns(A) and the other with 1
row and 124 column(B). i want to replace the maximum value in all columns of A
with each (single) column value of B.
How can i do it??
Thanks indeed in advance,
Eliza
Hi,
Try:
sapply(seq_len(ncol(A)),function(i) {indx - which(A[,i]%in% max(A[,i]));
A[,i][indx] - B[,i]; A[,i]})
A.K.
On Tuesday, October 29, 2013 10:16 AM, eliza botto eliza_bo...@hotmail.com
wrote:
Dear Users,
I have two matrices, one with 12 rows and 124 columns(A) and the other with 1
I tried to change fonts for lattice graphics:
trellis.par.set(list(axis.text=list(fontfamily=Monaco),
par.strip.text=list(fontfamily=Monaco)))
It works for axis.text but not for par.strip.text.
trellis.par.set() does not set par.strip.text neither strip.text while I
can change font for strip
Hi there,
I just updated my R to 3.0.2 and ran
R CMD check --as-cran on the just produced new version of fpc.
I got an error
Error: package MASS was built before R 3.0.0: please re-install it
- but I actually *did* re-install MASS without error just before that and
within R library(MASS)
Perhaps check your R_LIBS* variables?
http://stat.ethz.ch/R-manual/R-devel/library/base/html/libPaths.html
---
Jeff NewmillerThe . . Go Live...
DCN:jdnew...@dcn.davis.ca.us
To Eliza: What if the max in a column is not unique?
Given the small size of A, the solution given by Arun seems completely
adequate. However, I was wondering if it could be done without the
R-level loop in sapply by taking advantage of pmax() . Of course it
can. Here's code to illustrate how:
Thanks bert!!!it worked out perfectly well.thankyou onceagain,
Eliza
Date: Tue, 29 Oct 2013 08:30:13 -0700
Subject: Re: [R] maximum value replacement
From: gunter.ber...@gene.com
To: smartpink...@yahoo.com; eliza_bo...@hotmail.com
CC: r-help@r-project.org
To Eliza: What if the max in a
Dear R People:
I'm on a Centos 5 Red Hat system and I'm trying to install such packages as
Cairo, Rserve, etc.
However, I keep getting an error: sh:/bin/sh bad interpreter.
I'm logged in as root, so it shouldn't be a permissions error.
Has anyone else run into this, please?
Any help would be
... and while I'm being OCD, note that the which() call in my code can
and should be omitted. It's completely superfluous. :-(
-- Bert
On Tue, Oct 29, 2013 at 8:46 AM, eliza botto eliza_bo...@hotmail.com wrote:
Thanks bert!!!
it worked out perfectly well.
thankyou onceagain,
Eliza
Date:
On Tue, 29 Oct 2013, Erin Hodgess wrote:
I'm on a Centos 5 Red Hat system and I'm trying to install such packages as
Cairo, Rserve, etc.
However, I keep getting an error: sh:/bin/sh bad interpreter.
Erin,
A Web search shows several possible causes, including incorrect options in
Hi All,
if memory serves me well I recall some paper comparing the relative success in
getting mainstream acceptance (as mainstream as statistics can be) of both R
and Octave. I remember vaguely that the fact the development strategies (core
team vs one main developer) played a major role in
On Oct 29, 2013, at 4:38 AM, rm wrote:
Any ideas would be much appreciated; I suspect that this problem of
constructing the dummies applies not only to function coxph but to other
regression models in R as well. Effectively, my question is how to better
control for which dummies and
Hi,
So I just took an intro to R programming class and one of the lectures was on
Regular Expressions. I've been playing around with various R functions that use
Regular Expressions.
But this has me stumped. This was part of a quiz and I got it right through
understanding the syntax. But when
Sorry, this relates to ?mgcv-FAQ number 5, unfortunately. 'ti' terms
were introduced as a much better and cleaner way of allowing smooth main
effects and interactions where the interactions are based on 'te' terms:
this required some re-engineering of the tensor product smooth objects
(the
On Tue, 29 Oct 2013, Erin Hodgess wrote:
I'm on a Centos 5 Red Hat system and I'm trying to install such packages
as
Cairo, Rserve, etc.
However, I keep getting an error: sh:/bin/sh bad interpreter.
Erin, just to add to what Rich wrote, this may be a disk related result as
well. I get
(1) Am I correct in understanding that Heteroscedasticity is a problem for
Generalized Additive Models as it is for standard linear models? I am
asking particularly about the GAMs as implemented in the mgcv package.
Based upon my online search it seems that some forms of penalized splines
can
Thank you Simon that's quite helpful! I'll compare that with the GLMSS
models.
Best,
Collin.
(1) Am I correct in understanding that Heteroscedasticity is a problem
for
Generalized Additive Models as it is for standard linear models? I am
asking particularly about the GAMs
This covers the topic you mention, but from the perspective of
the role of the R Core team. The point about Octave is a single
sentence/footnote:
Fox, John. 2009. Aspects of the Social Organization and
Trajectory of the R Project. The R Journal 1/2: 5-13.
On Tue, Oct 29, 2013 at 1:13 PM, Lopez, Dan lopez...@llnl.gov wrote:
grep(^([a-z]+) +\1 +[a-z]+ [0-9],lines)
Your expression has a typo:
R grep(^([a-z]+) +\\1 +[a-z]+ [0-9],lines)
[1] 2
--
Sarah Goslee
http://www.functionaldiversity.org
__
Please read and follow the Posting Guide, in particular re plain text email.
You need to keep in mind that the characters in literal strings in R source
have to make it into RAM before the regex code can parse it. Since regex needs
a single backslash to escape normal parsing and interpret 1 as
Thanks for the good suggestions!
This is how I solved it (with lots of internet help):
Create a tempdir. Do a chmod 777 on it.
Within R
Sys.setenv(TMPDIR=/home/erin/tempdir)
install.packages(Cairo,depen=TRUE)
and all was well.
Thanks,
Erin
On Tue, Oct 29, 2013 at 12:24 PM,
This is from the other perspective
http://www.r-project.org/conferences/DSC-2001/Proceedings/Eaton.pdf
I can’t spot any direct comparison (and there is no mention of R in the
references), but I recall the ideas contrasting the two projects being bandied
about at the time. That discussion is
From ?regex
(do remember that backslashes need to be doubled when entering
R character strings, e.g. from the keyboard).
lines[grep(^([a-z]+) +\\1 +[a-z]+ [0-9],lines)]
[1] night night at 8
-
David L Carlson
Department of Anthropology
Texas AM University
Hi Jeff,
I was reviewing my old lecture notes and see that the professor did use \1 so I
think he was talking about regex in a non-platform specific context. But
obviously \\1 is the way to do it in R.
The examples you gave me to study really helped.
I was also going to ask how to identify
Christian Hennig ucakche at ucl.ac.uk writes:
I just updated my R to 3.0.2 and ran
R CMD check --as-cran on the just produced new version of fpc.
I got an error
Error: package MASS was built before R 3.0.0: please re-install it
- but I actually *did* re-install MASS without error just
Hello,
I would like to draw a map of Europe. Each country should be colored
depending on how it scores in an index called GPIndex.
Say a dark red for real bad countries a light red for those which are
not so bad, light blue for the fairly good ones and so on up to the
really good ones in a
Hi,
I'm having the following loop:
result - vector(list,100)
for (i in 1:max(dat$simNumber))
{
result[[i]]-survfit(Surv(dat[dat$simNumber==i,]$TAFD,dat[dat$simNumber==i,]$DV)~1)
}
In a next step, I would like to calculate the mean, 5% and 95% PI of the
Kaplan-Meier estimates of the 100
Dear Jeff,
thanks. Somehow R didn't install MASS where it later looked for it. I still
haven't understood properly what caused the problem but I managed to fix it now
(by specifying lib when installing it).
Best wishes,
Christian
*** --- ***
Christian Hennig
University College London,
Thanks everyone for the replies to my question. The issue turns out to be
that I'm on a Rocks cluster head node, and the Rocks distribution disables
alternate repos by defaut so it's using the Rocks-6.1 repo which has the
old R.
In the meantime I've built R 3.0.2 from source, which seems a better
Hello!
I have a very large matrix of results: 5x10. I saved it as RDS, but I
would also need to save it as txt or csv. Is there a way to do it? Now, with
write.table I am receiving an error:
Error in .External2(C_writetable, x, file, nrow(x), p, rnames, sep, eol, :
long vectors not
Dear expeRts,
If I specify group = as.factor(rep(1:2, each=n)) in the below
definition of dat, I get the expected behavior I am looking for. I
wonder why I
don't get it if group is *not* a factor... My guess was that,
internally, factors are treated as natural numbers (and this indeed
seems to be
Hi,
I am looking for some R function which will tell me, whether a particular sheet
in an Excel file (.xlsx/.xls) exists or not. I just need to get some TRUE/FALSE
type of answer.
Can somebody give me any pointer if such function exists or not?
Thanks and regards,
Check out this link for some examples
http://www.r-bloggers.com/maps-in-r-choropleth-maps/
Jean
On Tue, Oct 29, 2013 at 12:02 PM, palad...@trustindata.de wrote:
Hello,
I would like to draw a map of Europe. Each country should be colored
depending on how it scores in an index called
Think about it. How can one define a smooth term with a factor???
Further discussion is probably offtopic. Post on
stats.stackexchange.com if it still isn't obvious.
Cheers,
Bert
On Tue, Oct 29, 2013 at 1:16 PM, Marius Hofert
marius.hof...@math.ethz.ch wrote:
Dear expeRts,
If I specify group
Have you tried write.csv() or write.matrix()? I really don't know, but
they may be more efficient than write.table() with large matrices.
Jean
On Tue, Oct 29, 2013 at 2:27 PM, Petar Milin pmi...@gmail.com wrote:
Hello!
I have a very large matrix of results: 5x10. I saved it as RDS,
On 10/30/2013 04:02 AM, palad...@trustindata.de wrote:
Hello,
I would like to draw a map of Europe. Each country should be colored
depending on how it scores in an index called GPIndex.
Say a dark red for real bad countries a light red for those which are
not so bad, light blue for the fairly
On 10/29/13 19:44, peter dalgaard wrote:
SNIP
There really is no substitute for knowledge and understanding! Did it not occur
to you that the Windspeed column needs to enter into your analysis?
SNIP
Fortune!
cheers,
Rolf Turner
Hi, I'm running a model in ASRemlR and the conditional wald statistics table
is producing groupings for the marginality of some variables (A's B's and
C's in the tables below) that we are struggling to understand. From our
understanding of the model we are fitting and the reference manual, we were
Hi here is the solutions using XLConnect package:
library(XLConnect)
wb - loadWorkbook(path to your Excel file)
c(particular sheet name)%in%getSheets(wb)
Andrija
On Tue, Oct 29, 2013 at 9:26 PM, Ron Michael ron_michae...@yahoo.comwrote:
Hi,
I am looking for some R function which will
You could use the XLConnect package to do this. For example, something
like this might do the trick ...
library(XLConnect)
mysheet - Sheet4
wb - loadWorkbook(C:/temp/MyData.xlsx)
wbsheets - getSheets(wb)
mysheet %in% wbsheets
Jean
On Tue, Oct 29, 2013 at 3:26 PM, Ron Michael
Hello,
You can use the argument to write.csv or write.table append = TRUE to
write the matrix in chunks. Something like the following.
bigwrite - function(x, file, rows = 1000L, ...){
passes - NROW(x) %/% rows
remaining - NROW(x) %% rows
k - 1L
Hi Thorn,
it is not entirely clear (at least for me) what you want to accomplish.
an easy and fail safe way of extracting used terms in a (g)lm-object is
names(model.frame(l))
if you want to extract terms to finally select a model, have a look at
drop1 and/or MASS::dropterm
Hth
Am 28.10.2013
Hi,
Try either:
res1 - apply(mydata[,1:2],2,mean)
res2 - colMeans(mydata[,1:2])
identical(res1,res2)
#[1] TRUE
# Also if you need to find means for each group (Ungrazed vs. Grazed)
by(mydata[,-3],mydata[,3],colMeans)
#or if column names are V1, V2, V3
aggregate(.~V3,mydata,mean)
#or
On 29/10/2013 20:42, Rui Barradas wrote:
Hello,
You can use the argument to write.csv or write.table append = TRUE to
write the matrix in chunks. Something like the following.
That was going to be my suggestion. But the reason long vectors have not
been implemented is that is rather
Hello,
On Oct 29, 2013, at 10:16 PM, Prof Brian Ripley rip...@stats.ox.ac.uk wrote:
On 29/10/2013 20:42, Rui Barradas wrote:
Hello,
You can use the argument to write.csv or write.table append = TRUE to
write the matrix in chunks. Something like the following.
That was going to be my
Hi,
I need to read an Excel file which can be available in following link:
http://www45.zippyshare.com/v/43626889/file.html
Now I wanted to read the 1st sheet of this Excel file. Below are my code so far
(I saved that file in 'F:' drive):
library(XLConnect)
Loading required package: rJava
Hi Petar,
If you're going to share this matrix across R sessions, save()/load() is
probably one of your best options.
Otherwise, you could try the rhdf5 package from Bioconductor:
1. Install the package with:
source(http://bioconductor.org/biocLite.R;)
biocLite(rhdf5)
2. Then:
Dear list,
I need to use the function deparse in a specific situation.
But, it always replace any occurence of by \.
For example:
arg - deparse(args(cov), width.cutoff = 100L)[1]
arg
[1] function (x, y = NULL, use = \everything\, method =
c(\pearson\, \kendall\, \spearman\))
Some
On 13-10-29 7:42 PM, Jose Claudio Faria wrote:
Dear list,
I need to use the function deparse in a specific situation.
But, it always replace any occurence of by \.
No it doesn't. That's just how print() displays quotes. Use cat() and
you can see what's really there.
Duncan Murdoch
Dear All,
is it known that source works much faster in R 2.15.2 than in R 3.0.2 ?
In the example below I observe e.g. for a data.frame with 10^7 rows the
following timings:
R version 2.15.2 Patched (2012-11-29 r61184)
length: 1e+07
user system elapsed
62.040.22 62.26
R version
I am not sure that I fully understand your question, but the mvrnorm
function in the MASS package may do what you want. It will generate
multivariate normal data with a specified correlation(covariance). If
that is not what you want then try to explain a bit more about what
you want your final
Hi all,
I am trying to add a color legend to my plot. As an example I am giving you a
bit of code that you can run.
I am sharing for everyone a small data snipset that you can load
https://www.dropbox.com/s/fh8jhwujgunmtrb/DataToPlotAsImage.Rdata
load(DataToPlotAsImage.Rdata)
require(plotrix)
Hi,
May be this helps.
library(plyr)
res - join_all(lapply(my.list,function(x)
as.data.frame(t(unlist(x,type=full)
res
# AICc Intercept Burned StandAge TreeDensity RoadDensity Intercept.SE
#1 108.2303 -1.3358063 1.351866 0.05606852 -0.1886327 -0.03904008 0.8392739
#2
Dear R-users,
I am trying to fit my data using one or more horizontal lines. If my data
is in y, I understand that lm(y~1) will fit a single horizontal line at
mean(y). However, I want to try and fit the data with multiple horizontal
lines if that reduces the error while still keeping the number
Hello again!
I'm using python with a module rpy2 to call functions from R.
It works fine on built in R functions like rnorm.
However, I would like to access user-defined functions as well. For those
of you who use this, I have:
import rpy2.robjects as R
x = R.r.buzz(3)
R object as no
61 matches
Mail list logo
