On Tue, 18 Nov 2014 22:06:03 -0600
David Doyle kydaviddo...@gmail.com wrote:
Thank you to David and David for their help. The code below
generated what I needed.
library(circular)
mydata - read.table(http://doylesdartden.com/R/Joints.csv;,
header=TRUE, sep=,,)
x -
Scott Rifkin sarif...@ucsd.edu writes:
Ista,
On the one hand I'd like it to be as flexible as possible so the
students could really come up with whatever they like. On the other
hand, restricting their choices probably would make it easier to do
the backend. The goal would be to get them
Dear all,
I use R 3.1.1 for Windows.
I performed two different meta-analysis assessing the prognostic value of
two different tests in patients with coronary artery disease. The study
included in the two analysis are different.
The variable of interest in dichotomous (normal/abnormal result) for
On 21/11/2014 08:51, Mario Petretta wrote:
Dear all,
I use R 3.1.1 for Windows.
I performed two different meta-analysis assessing the prognostic value of
two different tests in patients with coronary artery disease. The study
included in the two analysis are different.
That makes life
Those hazard ratios and CIs seem a bit strange. On the log-scale, they should
be symmetric, but they are not. Could be due to heavy rounding though. At any
rate, it comes down to this:
hr- c(3.12, 1.15)
ci.lb - c(2.2, 1.03)
ci.ub - c(4.1, 2.6)
meta - c(1,2)
### log-transform hazard ratios
I am aware of the fact that bootstrapping produces different CIs with every
run. I still believe that there is a difference between both types of
procedures. My understanding is that setting w in the boot() function
influences the importance of observations or how the bootstrap selects
the
No. Just use the circular() function to specify that your data are in degrees
and clockwise and the graph will be labeled that way.
David C (I was beginning to think that this thread was only for Davids).
-Original Message-
From: r-help-boun...@r-project.org
Anyone can help?
On Thu, Nov 20, 2014 at 6:44 PM, Jennifer Sabatier
plessthanpointohf...@gmail.com wrote:
Hi R-Help,
So, I will try to provide a reproducible example...I basically made a
dummy spreadsheet that contains the same number of tabs as the spreadsheet
I am really interested in.
Add a print(ts2$i) to see if that really gives what you think it does.
Then take it from there.
Henrik
On Nov 21, 2014 7:19 AM, Jennifer Sabatier plessthanpointohf...@gmail.com
wrote:
Anyone can help?
On Thu, Nov 20, 2014 at 6:44 PM, Jennifer Sabatier
plessthanpointohf...@gmail.com wrote:
On Fri, Nov 21, 2014 at 10:36 AM, Henrik Bengtsson h...@biostat.ucsf.edu
wrote:
print(ts2$i)
You're right. It's not doing what I think it's supposed to be doing. The
print(ts2$i) gives NULL:
[1] Cash Flows211
NULL
[1] Cash Flows211 - 2014-11-21.csv
7: getURL(sheet@cellsfeed, curl =
I got it!
for (i in sheetz$names) {
print(i)
file.name - paste(i, - , Sys.Date(), .csv, sep=)
print(file.name)
#traceback()
i - sheetAsMatrix(ts[[i]], header = TRUE, as.data.frame = TRUE, trim =
TRUE)
write.csv(i, file.name)
i
}
I have to use ts[[i]] rather than ts$i (or ts2, either)
Thanks!
Most field mapping data I have collected has been
either in quadrant form (rarely) or more commonly in azimuthal form
(0-360 degrees order clockwise from the top).
You can specify that when making the 'circular' object, by using the
zero and roation arguments. Compare the plots made by the
Hi Dimitri,
I am not sure I completely understand the issue but here is what I think you
want (using my variable names :) )
library(ggplot2)
dat1 - data.frame(aa = sample(1:20, 100, replace = TRUE), bb = 1:100 )
p - ggplot(dat1, aes(aa, bb)) + geom_point()
p - p +
I got an error when trying to extract a 1-column subset of a data frame (called
my.output) created by dplyr/summarize. The ncol() function says that
my.output has 4 columns, but my.output[4] fails. Note that converting
my.output using as.data.frame() makes for a happy ending.
Is this the
Many many thanks to Michael Dewey and to Viechtbauer Wolfgang for the kindly
and useful replay !!
I only ask to Wolfgang if I should log-transform hazard ratios and compute
standard error only for the summary hazard ratio estimates or, as suggested by
Michael, for each single study, combining
Your code in creating 'frm' is not working for me and it is complicated enough
that I don't want to work it out. See ?dput for a better way to supply data.
Also see:
https://github.com/hadley/devtools/wiki/Reproducibility
On Nov 21, 2014, at 6:52 AM, ivan wrote:
I am aware of the fact that bootstrapping produces different CIs with every
run. I still believe that there is a difference between both types of
procedures. My understanding is that setting w in the boot() function
influences the importance of
1 my_min= min(outcome_data[which(outcome_data$State==my_state),11],na.rm=TRUE)
2 print(my_min)
3 jkr=0
4 if (jkr= 4706) {jkr=jkr+1
5 if (identical(outcome_data[jkr,11],my_min)
identical(outcome_data[jkr,7],my_state)) {
6 print((outcome_data[jkr,2]))
7 break
8 }
Dear Experts,
My
Seems to me you probably wanted a while in line 4.
N.b. There's also a missing } and totally messed up formatting :-(
B.
On Nov 21, 2014, at 1:19 PM, Aditya Singh aps...@yahoo.com wrote:
1 my_min=
min(outcome_data[which(outcome_data$State==my_state),11],na.rm=TRUE)
2 print(my_min)
3
Dear guRus
How can I round of time in R to the nearest 30th minute please?
For example suppose if
Sys.time()
[1] 2014-11-21 22:49:05.59042 GMT
then I would like a function that outputs 22:30:00.
if Sys.time is 12:13:22 then I would like to get 12:00:00 etc.
Any help would be appreciated.
Many
Hi all,
I need some assistance regarding the use of C project (set of programming
files) as R functions in Windows OS.
By now I really would like to avoid package-building.
What are the steps to undergo or where can I check to perform that
successfully?
Thanks in advance,
My best
On 20 Nov 2014, at 23:02 , Duncan Murdoch murdoch.dun...@gmail.com wrote:
On 20/11/2014, 4:35 PM, Ista Zahn wrote:
This prompted me to do something I've been meaning to do for a long
time, which is to try to get the R website updated in terms of visual
style. It really is showing it's age
Hi everyone.
I have a non-linear model specified as this (6 variables, a0, S, K, p0, p1, p2):
fit - nlsLM(formula = dat ~ a0 * exp(-S*(x - 250)) + K + ( C ) )
where C = (p0*exp(-0.5*((x-p1)/p2)^2))
The problem is that I do not know in advance how many component (C) I will have
in the model.
Sorry I forgot to mention it clearly. I like to round it to the
nearest 30th minute that is past. So 12:28:59 will be again 12:00:00
and
12:59:59 will be 12:30:00 etc. Apologies for the lack of clarity in
the beginning.
Many thanks
Raghu
On Fri, Nov 21, 2014 at 10:52 PM, Raghuraman Ramachandran
On Nov 21, 2014, at 2:52 PM, Raghuraman Ramachandran wrote:
Dear guRus
How can I round of time in R to the nearest 30th minute please?
For example suppose if
Sys.time()
[1] 2014-11-21 22:49:05.59042 GMT
then I would like a function that outputs 22:30:00.
if Sys.time is 12:13:22 then
On Nov 21, 2014, at 2:55 PM, Raghuraman Ramachandran wrote:
Sorry I forgot to mention it clearly. I like to round it to the
nearest 30th minute that is past. So 12:28:59 will be again 12:00:00
and
12:59:59 will be 12:30:00 etc. Apologies for the lack of clarity in
the beginning.
That's
Is there a way to shrink the size of RandomForest-class (an S4
object), so that it requires less memory during run-time and less disk
space for serialization?
On my system the data slot is about 2GB, which is causing problems,
and I'd like to see whether predict() works without it.
# example
Perhaps this ...
roundToHalf - function(t) {
t - as.POSIXlt(t)
if (t$min 15) {
t$min - 0
}
else if (t$min 45) {
t$min - 30
}
else {
t$min - 0
t$hour - t$hour + 1
On 21/11/2014, 3:02 PM, Diogo André Alagador wrote:
Hi all,
I need some assistance regarding the use of C project (set of programming
files) as R functions in Windows OS.
By now I really would like to avoid package-building.
What are the steps to undergo or where can I check to
On Nov 21, 2014, at 3:19 PM, David Winsemius wrote:
On Nov 21, 2014, at 2:55 PM, Raghuraman Ramachandran wrote:
Sorry I forgot to mention it clearly. I like to round it to the
nearest 30th minute that is past. So 12:28:59 will be again 12:00:00
and
12:59:59 will be 12:30:00 etc.
Hi all,
I think this is probably a very simple question related to ANOVA or
quantile regression.
Let's say if I'm doing ANOVA on 2 factors (f_a and f_b) and trying to
figure out whether they have an impact on the output. The result shows me
that f_a has a significant impact while f_b does not as
Diogo,
For a gentle introduction to package creation see Hadley Wickham's book in
progress on the subject at http://r-pkgs.had.co.nz. This site is particularly
accessible if you use RStudio, but is not dependent on its use.
I use very simple packages that are not designed to be used by others
Dear Boris and R-Experts,
I have a variable my_state which is a 2-letter character string telling which
American state the user inputs. This I am do a if(identical(database
entry,my_state)) to check for occurrences in the database.
The problem is that the database entry[i,j] shows as a factor
This is off topic here.Post to a statistics forum like
stats.stackexchange.com instead -- or talk to your professor or TA (if
you're a student).
Cheers,
Bert
Bert Gunter
Genentech Nonclinical Biostatistics
(650) 467-7374
Data is not information. Information is not knowledge. And knowledge
is
Nueva publicación en blog: Ofertas de trabajo, Buscamos profesores para Curso
de R+Hadoop (BigData Analytics), por luispenya en 21/11/14 13:04h
Ver el blog en:
http://r-es.org/tiki-view_blog_post.php?blogId=2postId=92
Si no desea recibir estas notificaciones siga este enlace:
Muchas gracias a todos, con la dodificacion UTF-8 todo queda solucionado.
Una pena que en pleo año 2014 aun existan estas incompatibilidades y más
para gente como yo que acabamos utilizando los tres principales OS.
Muchas gracias de nuevo.
El 18 de noviembre de 2014, 17:05, Carlos J. Gil
36 matches
Mail list logo