Ana Marija
2020-May-21 22:23 UTC
[R] how to show percentage of individuals for two groups on histogram?
Hello, I have a data frame like this:> head(a)FID IID FLASER PLASER DIABDUR HBA1C ESRD pheno 1 fam1000-03 G1000 1 1 38 10.2 1 control 2 fam1001-03 G1001 1 1 15 7.3 1 control 3 fam1003-03 G1003 1 2 17 7.0 1 case 4 fam1005-03 G1005 1 1 36 7.7 1 control 5 fam1009-03 G1009 1 1 23 7.6 1 control 6 fam1052-03 G1052 1 1 32 7.3 1 control> dim(a)[1] 1698 8 I am doing histogram plot via: ggplot(a, aes(x=HBA1C, fill=pheno)) + geom_histogram(binwidth=.5, position="dodge") there is 848 who have "case" in pheno column and 892 who have "control" in pheno column. I would like to have on y-axis shown percentage of individuals which have either "case" or "control" in pheno instead of count. Please advise, Ana
Ana Marija
2020-May-21 22:52 UTC
[R] how to show percentage of individuals for two groups on histogram?
the result would basically look something like this on in attach or the overlay of those two plots On Thu, May 21, 2020 at 5:23 PM Ana Marija <sokovic.anamarija at gmail.com> wrote:> > Hello, > > I have a data frame like this: > > head(a) > FID IID FLASER PLASER DIABDUR HBA1C ESRD pheno > 1 fam1000-03 G1000 1 1 38 10.2 1 control > 2 fam1001-03 G1001 1 1 15 7.3 1 control > 3 fam1003-03 G1003 1 2 17 7.0 1 case > 4 fam1005-03 G1005 1 1 36 7.7 1 control > 5 fam1009-03 G1009 1 1 23 7.6 1 control > 6 fam1052-03 G1052 1 1 32 7.3 1 control > > > dim(a) > [1] 1698 8 > > I am doing histogram plot via: > ggplot(a, aes(x=HBA1C, fill=pheno)) + geom_histogram(binwidth=.5, > position="dodge") > > there is 848 who have "case" in pheno column and 892 who have > "control" in pheno column. > > I would like to have on y-axis shown percentage of individuals which > have either "case" or "control" in pheno instead of count. > > Please advise, > Ana-------------- next part -------------- A non-text attachment was scrubbed... Name: Screen Shot 2020-05-21 at 5.49.37 PM.png Type: image/png Size: 52888 bytes Desc: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20200521/740baf08/attachment.png>
Jim Lemon
2020-May-22 04:08 UTC
[R] how to show percentage of individuals for two groups on histogram?
Hi Ana,
My apologies for the pedestrian graphics, but it may help.
# a bit of fake data
aafd<-data.frame(FID=paste0("fam",1000:2739),
IID=paste0("G",1000,2739),FLASER=rep(1,1740),
PLASER=c(rep(1,892),rep(2,848)),
DIABDUR=sample(10:50,1740,TRUE),
HBAIC=rnorm(1740,mean=7.45,sd=2),ESRD=rep(1,1740),
pheno=c(rep("control",892),rep("case",848)))
par(mfrow=c(2,1))
casepct<-table(cut(aafd$HBAIC[aafd$pheno=="case"],breaks=0:15))
controlpct<-table(cut(aafd$HBAIC[aafd$pheno=="control"],breaks=0:15))
par(mar=c(0,4,1,2))
barpos=barplot(100*casehist,names.arg=names(casepct),col="orange",
space=0,ylab="Percentage",xaxt="n",ylim=c(0,25))
text(mean(barpos),23,
"Cases: n=848, nulls=26, median=7.3, mean=7.45, sd=1.96")
box()
par(mar=c(3,4,0,2))
barplot(100*controlhist,names.arg=names(controlpct),
space=0,ylab="Percentage",col="orange",ylim=c(0,25))
text(mean(barpos),23,
"Controls: n=892, nulls=7, median=7.3, mean=7.45, sd=1.12")
box()
Jim
On Fri, May 22, 2020 at 9:08 AM Ana Marija <sokovic.anamarija at
gmail.com> wrote:>
> the result would basically look something like this on in attach or
> the overlay of those two plots
>
>
> On Thu, May 21, 2020 at 5:23 PM Ana Marija <sokovic.anamarija at
gmail.com> wrote:
> >
> > Hello,
> >
> > I have a data frame like this:
> > > head(a)
> > FID IID FLASER PLASER DIABDUR HBA1C ESRD pheno
> > 1 fam1000-03 G1000 1 1 38 10.2 1 control
> > 2 fam1001-03 G1001 1 1 15 7.3 1 control
> > 3 fam1003-03 G1003 1 2 17 7.0 1 case
> > 4 fam1005-03 G1005 1 1 36 7.7 1 control
> > 5 fam1009-03 G1009 1 1 23 7.6 1 control
> > 6 fam1052-03 G1052 1 1 32 7.3 1 control
> >
> > > dim(a)
> > [1] 1698 8
> >
> > I am doing histogram plot via:
> > ggplot(a, aes(x=HBA1C, fill=pheno)) + geom_histogram(binwidth=.5,
> > position="dodge")
> >
> > there is 848 who have "case" in pheno column and 892 who
have
> > "control" in pheno column.
> >
> > I would like to have on y-axis shown percentage of individuals which
> > have either "case" or "control" in pheno instead
of count.
> >
> > Please advise,
> > Ana
> ______________________________________________
> R-help at r-project.org mailing list -- To UNSUBSCRIBE and more, see
> 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.