While putting my R code into functions, I've encountered a ddply function
nesting issue and need a bit of advice on the proper way to fix it. I've tried
several approahces, but neither worked and I need to have the ability to
include the "cut", "range", and "fullseq" methods within ddply. (For a bit of
that explanation refer to
http://finzi.psych.upenn.edu/Rhelp08/2009-February/187331.html)
Thus, in order to preserve that functionality, and put my code within
functions, I needed to have an architecture similar to the following
implemented, where you end up running:
function_nesting()
Unfortunately this produced errors within the ddply where it does not appear to
be recognizing or allowing variables or functions to be processed within side
its function.
Thank you for any advice about how to proceed forward.
determine_counts<-function()
{
min_range<-1
max_range<-30
bin_range_size<-5
Me_df<-data.frame(Data = c(1:15), Person = "Me")
You_df<-data.frame(Data = c(10:20), Person = "You")
Them_df<-data.frame(Data = c(15:25), Person = "Them")
Group_df_tmp<-rbind(Me_df,You_df)
Group_df<-rbind(Group_df_tmp,Them_df)
Group_df$Person <- factor(Group_df$Person, levels = c("Them", "You",
"Me"))
#counts <- ddply(Group_df, .(cut(Data, breaks=fullseq(range(Data), 5)),
Person), nrow)
# Approach 1
counts <- ddply(Group_df, .(cut(Data,
breaks=fullseq(range(c(Group_df$Data, min_range, max_range)), bin_range_size)),
Person), nrow)
# Approach 2
range_tmp<-range(c(Group_df$Data, min_range, max_range))
counts <- ddply(Group_df, .(cut(Data, breaks=fullseq(range_tmp,
bin_range_size)), Person), nrow)
names(counts) <- c("Bin", "Person", "Frequency")
qplot(Person, Frequency, data = counts, fill = Person, geom="bar",
stat="identity", width = 0.9, xlab="Person") + facet_grid(. ~ Bin)
}
function_nesting<-function()
{
determine_counts()
}
However, if the code is just run straight through without being nested it works
fine:
min_range<-1
max_range<-30
bin_range_size<-5
Me_df<-data.frame(Data = c(1:15), Person = "Me")
You_df<-data.frame(Data = c(10:20), Person = "You")
Them_df<-data.frame(Data = c(15:25), Person = "Them")
Group_df_tmp<-rbind(Me_df,You_df)
Group_df<-rbind(Group_df_tmp,Them_df)
Group_df$Person <- factor(Group_df$Person, levels = c("Them", "You",
"Me"))
#counts <- ddply(Group_df, .(cut(Data, breaks=fullseq(range(Data), 5)),
Person), nrow)
counts <- ddply(Group_df, .(cut(Data,
breaks=fullseq(range(c(Group_df$Data, min_range, max_range)), bin_range_size)),
Person), nrow)
Unfortunately this is not within a function, so thanks again for any advice on
how to approach this issue.
______________________________________________
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, reproducible code.
