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 ce
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 w
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 t
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
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 et
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 ch
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
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 with
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.
>
Tha
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 1
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
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. T
On 20 Nov 2014, at 23:02 , Duncan Murdoch 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 and could use a face
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 regard
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.
Ma
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 wrote:
> 1 my_min=
> min(outcome_data[which(outcome_data$State==my_state),11],na.rm=TRUE)
> 2 print(my_min)
> 3 jkr=0
>
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 comp
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"
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
http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproduc
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 t
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 in
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 + geo
> 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
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)
T
On Fri, Nov 21, 2014 at 10:36 AM, Henrik Bengtsson
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 = getCurlCon(con), follow
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"
wrote:
> Anyone can help?
>
> On Thu, Nov 20, 2014 at 6:44 PM, Jennifer Sabatier <
> plessthanpointohf...@gmail.com> wrote:
>
> > Hi R-Help,
> >
>
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
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 [mailto:r-help-boun...@r-proj
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 obse
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 rat
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 simpl
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
Scott Rifkin 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 to realize t
On Tue, 18 Nov 2014 22:06:03 -0600
David Doyle 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 <- circular(mydata$JointsRad)
> r
34 matches
Mail list logo