On 04/02/2015 10:16 AM, michael.laviole...@dhhs.state.nh.us wrote:
I'm preparing some reports for substate regions from BRFSS survey data. I
can get estimates easily enough, but am having problems putting the results
in convenient form. Here's some code using the New Hampshire portion of the
Shanley:
I am trying to run STAR logistic/binomial model using R2BayesX . A
warning message came up every time when i included a random effect,
bs='re', (unstructured spatial effect) into the model. The warning
message is : Warning message: running command
Hello,
Aren't the levels of your example wrong? If the levels are
levels=c('a','b','c'), not c('b', 'a', 'c'), then the following will do
the job.
unname(unlist(tapply(dat$D, dat$S, order)))
Hope this helps,
Rui Barradas
Em 04-02-2015 19:34, Tom Wright escreveu:
Given a dataframe:
Hi,
To complete the last part:
dist - function(x)
+ {
+
return(c(rep(x[1],15),rep(x[2],15),rep(x[3],25),rep(x[4],20),rep(x[5],25)))
+ }
x - dist(c(2,3,2,5.3,7.3))
x
[1] 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 2.0 3.0 3.0
3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0 3.0
[29]
Dear R users,
I am having trouble finding a package and function to remove DNA sequences
from a fasta file that are 99% similar and/or create an output of the
remaining unique sequences. I found the uniquefasta function in phytools
but R can't find this function and also doesn't allow me to set
Given a dataframe:
dat-data.frame(S=factor(c(rep('a',2),rep('b',1),rep('c',3)),levels=c('b','a','c')),
D=c(5,1,3,2,3,4))
where S is a subject identifier and D a visit (actually a date in my
real dataset). I would like to generate another column giving the visit
number
No problem with disguise, I'm looking for pretty.
On Wed, 2015-02-04 at 12:06 -0800, Bert Gunter wrote:
tapply() (of which by() is essentially a wrapper) **is** a (disguised)
loop (at the R level, of course).
Cheers,
Bert
Bert Gunter
Genentech Nonclinical Biostatistics
(650)
Is there a way for me to use the Swirl course Data_Visualization? This
is my main reason for learning R, for the plotting ability but can't use
the course. I've already been through the following (Programming,
Programming Alt, Getting_and_Cleaning_Data).
I get the following
Thanks, I was not aware of order().
I did deliberately mess up the order of S. The following example breaks
your solution
dat_2-data.frame(S=factor(c('a','c','a','b','c','c')),
D=c(5,3,1,3,2,4))
which should give the answer c(2,2,1,1,2,3)
Your solution does indicate that
For the introduction to R I strongly suggest you look at the materials
published by software-carpentry www.software-carpentry.org. The lessons
are all open-source, hosted on github and are under active development.
On Tue, 2015-02-03 at 12:08 +0100, Michael Haenlein wrote:
Dear all,
I am
tapply() (of which by() is essentially a wrapper) **is** a (disguised)
loop (at the R level, of course).
Cheers,
Bert
Bert Gunter
Genentech Nonclinical Biostatistics
(650) 467-7374
Data is not information. Information is not knowledge. And knowledge
is certainly not wisdom.
Clifford Stoll
How about?
ave(dat$D, dat$S, FUN=order)
[1] 2 1 1 1 2 3
ave(dat_2$D, dat_2$S, FUN=order)
[1] 2 2 1 1 1 3
Note, your answer for the second example is incorrect since row 2 (c, 3) and
row 5 (c, 2) are both assigned 2.
-
David L Carlson
Department of
Of course you are correct the second answer should be
c(2,2,1,1,1,3)
Thanks everyone.
On Wed, 2015-02-04 at 20:53 +, David L Carlson wrote:
How about?
ave(dat$D, dat$S, FUN=order)
[1] 2 1 1 1 2 3
ave(dat_2$D, dat_2$S, FUN=order)
[1] 2 2 1 1 1 3
Note, your answer for the second
dat-data.frame(S=factor(c(rep('a',2),rep('b',1),rep('c',3)),levels=c('b','a','c')),
+ D=c(5,1,3,2,3,4))
dat
S D
1 a 5
2 a 1
3 b 3
4 c 2
5 c 3
6 c 4
dat$visit - ave(seq(nrow(dat)), dat$S, FUN = seq_along)
dat
S D visit
1 a 5 1
2 a 1 2
3 b 3 1
4 c 2 1
5 c 3
Dear Z and Nikolaus,
Thank you for your help.
Here is my R codes:
data(MunichBnd)
N - length(MunichBnd); n - N*5
dat - data.frame(x1 = runif(n, -3, 3),id = as.factor(rep(names(MunichBnd),
length.out = n)))
dat$sp - with(dat, sort(runif(N, -2, 2), decreasing = TRUE)[id])
dat$re_poi - with(dat,
It's hard to say from your description what the situation is. The error
simply means the plot area is too small for the figure margins to fit. Try
closing the graphics (plot) window before you run the section that causes
the error, or you can try maximizing the plotting window.
You can also
I have some materials at http://tutorials.iq.harvard.edu , you're welcome
to use or adapt.
On Feb 3, 2015 6:12 AM, Michael Haenlein haenl...@escpeurope.eu wrote:
Dear all,
I am Professor at a business school and I would like to develop a course
about quantitative research using R.
My
Sorry Jim,
That messes up on S=='a'. Should be 2,1 not 1,2
Neat answer though and looks like it should be pretty quick after I
apply some sorting.
On Wed, 2015-02-04 at 15:37 -0500, jim holtman wrote:
dat-data.frame(S=factor(c(rep('a',2),rep('b',1),rep('c',3)),levels=c('b','a','c')),
+
A useful technique when it is easy to compute a vector from an ordered
data.frame but you need to do it for an unordered one is to compute the
order
vector 'ord', compute the vector from df[ord,], and use df[ord,...] -
vector
to reorder the vector. In your case you could do:
A potential problem with
ave(dat_2$D, dat_2$S, FUN=order)
is that it will silently give the wrong answer
or give an error if dat_2$D is not numeric.
E.g., if D is a Date vector we get
dat_3 - dat_2[,1:2]
dat_3$D - as.Date(paste0(2015-02-, dat_2$D))
with(dat_3, ave(D, S, FUN=order))
Steve (and any others still paying attention to this thread),
Larry Wall (author of Perl) said something along the lines of:
things that are similar should look similar, things that are different
should look different.
Ironically one of the first places I saw that quote was in a Perl vs.
Python
On Feb 4, 2015, at 11:42 AM, Nick Jeffery wrote:
Dear R users,
I am having trouble finding a package and function to remove DNA sequences
from a fasta file that are 99% similar and/or create an output of the
remaining unique sequences. I found the uniquefasta function in phytools
but R
Hi All,
I need to loop through and download the past 10 years of met data to a
temporary directory. I then need to unzip it and place it into another
directory.
year = (2005:2015)
for (i in year)
tmpdir = tempdir()
file[i] = file.path(tmpdir, sprintf('724927-23285-%4i.gz', i))
url =
On Sat, Jan 31, 2015 at 1:36 PM, Fisher Dennis fis...@plessthan.com wrote:
followed by a series of errors:
gcc -m32 -shared -s -static-libgcc -o ../compiled/Windows32.so tmp.def
ConvertSAS.o CKHashTable.o readstat_convert.o readstat_bits.o readstat_io.o
readstat_sas.o
Dear R users,
I have three plots, so I tried, for exmple,
par(mfrow=c(2,2))
y1 - rnorm(100)
y2 - rnorm(100)
y3- rnorm(100)
plot(y1);plot(y2);plot(y3)
Here, I'd like to put a legend on the bottom right hand side (empty space).
is it possible?
Thanks for helping,
Kathryn Lord
Buenas.
Abajo pongo la salida de un modelo de cox , dónde he estratificado por
una variable de país (Countryb) y por otra (Q6). Además hay interacción
entre la variable mobilityPDurG2 (es una variable 0,1, y 0 es la
categoría de referencia) país.
La categoría de referencia para país es
¿Es posible que lo que busques sea la función get()?
Un saludo
Isidro Hidalgo Arellano
Observatorio Regional de Empleo
Consejería de Empleo y Economía
http://www.jccm.es
-Mensaje original-
De: R-help-es [mailto:r-help-es-boun...@r-project.org] En nombre de Isa
García Barón
Enviado el:
Hola Xavier.
Me alegro mucho de que el BioStsFLOSS os haya servido para desfacer el
entuerto.
Uno de los objetivos del proyecto era conseguir dotar de funcionalidad a R en
entornos hostiles (aulas de inform�tica, redes corporativas, ...)
Ahora mismo estamos con la fase Beta de EpiLinux 4 y, te
Hello R-Help folks. I'm fairly new to R and particularly the graphing
functions. The figures I have made look great though and so I want to use it
more.
I am using bargraph.CI and have made several charts successfully using raw
data. I have a figure, however that doesn't lend itself to the
Hola, espero explicar bien el problema que tengo.
Estoy intentando hacer loops para glm's. El problema vieneal nombrar cada
glm de una manera y realizar la seleccion por AIC mediante la función step.
Cuando realizo los glm utilizo las funcion assign y paste0, para nombrar a
cada uno distinto:
I'm preparing some reports for substate regions from BRFSS survey data. I
can get estimates easily enough, but am having problems putting the results
in convenient form. Here's some code using the New Hampshire portion of the
public BRFSS SMART data:
library(foreign)
library(survey)
# download
Hi everyone,
I have a data.table with 200 columns and few million rows and am trying to
calculate the .1 and .9 quantiles for each row across all 200 columns.
I have found different ways to do this, all with different performances. The
examples I used are below. I wonder whether there is a
32 matches
Mail list logo