[R] Help for x axis

[David L Carlson]

Try these:

aggData <- aggregate(age~smoke+gender, sampleData, function(x) c(mean=mean(x), 
	stderr=sd(x)/sqrt(length(x))))
aggData
plotInfo <- barplot(aggData$age[,1], ylim=c(0, max(rowSums(aggData$age))))
axis(1, c(0, plotInfo), c("Gender", "-", "-", "+", "+"), line=.75, 
	lwd=0, cex.axis=1.25, xpd=TRUE)
axis(1, c(0, plotInfo), c("Smoke", "-", "+", "-", "+"), line=2, 
	lwd=0, cex.axis=1.25, xpd=TRUE)
top <- aggData$age[,1]+aggData$age[,2]
bottom <- aggData$age[,1]-aggData$age[,2]
arrows(plotInfo, bottom, plotInfo, top, length=.15, angle=90, code=3)


boxplot(age~smoke+gender, sampleData, xaxt="n")
axis(1, 0:4, c("Gender", "-", "-", "+", "+"), line=.75, 
	lwd=0, cex.axis=1.25, xpd=TRUE)
axis(1, 0:4, c("Smoke", "-", "+", "-", "+"), line=2, 
	lwd=0, cex.axis=1.25, xpd=TRUE)


David

-----Original Message-----
From: Olivier [mailto:olivier.lerouzic at ymail.com] 
Sent: Monday, November 17, 2014 4:39 PM
To: David L Carlson
Subject: Re: [R] Help for x axis

Thank you very much, it is all I want to do. Is it possible with showing 
the error-bars or in a boxplot?
Best regards,

Olivier Le Rouzic

On 2014-11-17, 4:07 PM, David L Carlson wrote:
> This should get you started:
>
>> aggData <- aggregate(age~smoke+gender, sampleData, mean)
>> aggData
>    smoke gender      age
> 1     0      0 39.47231
> 2     1      0 40.09020
> 3     0      1 39.59814
> 4     1      1 42.04092
>> plotInfo <- barplot(aggData$age)
>> axis(1, c(0, plotInfo), c("Gender", "-", "-", "+", "+"), line=.75,
> + lwd=0, cex.axis=1.25, xpd=TRUE)
>> axis(1, c(0, plotInfo), c("Smoke", "-", "+", "-", "+"), line=2,
> + lwd=0, cex.axis=1.25, xpd=TRUE)
>
> To adapt it you will need to read the manual pages for barplot() and axis() and the page on graphical parameters par(). In particular, you will have to allocate more space at the bottom of the plot if you want to add more lines.
>
> -------------------------------------
> David L Carlson
> Department of Anthropology
> Texas A&M University
> College Station, TX 77840-4352
>
>
>
> -----Original Message-----
> From: r-help-bounces at r-project.org [mailto:r-help-bounces at r-project.org] On Behalf Of Olivier
> Sent: Monday, November 17, 2014 9:24 AM
> To: r-help at r-project.org
> Subject: [R] Help for x axis
>
> Hi,
> I want to customize x axis to scientific data. I do experiments with
> different triggers. As others publications, I want that there is one
> line for each trigger with the sign "-" or "+" to show if the trigger is
> used or no. You will find attached an exemple.
> Please find below a data.frame you could use to explain me.
> Thank you for your response,
> Olivier
>
>
>
> set.seed(3)
> sampleData <- data.frame(id = 1:100,gender = sample(c("0", "1"), 100,
> replace = TRUE), smoke = sample (c("0","1"), 100, replace=TRUE), age =
> rnorm(100, 40, 10))
> summary(sampleData)
>
> -> I want to give results with histograms or box.plot (age according to
> sex and smoking status)
> -> x axis may be like something like this :
>
> Gender        -          -         +        +
> Smoke         -          +        -         +
>
>