Как я могу получить графики плотности в бобовом заговоре (скрипке), чтобы они соответствовали друг другу?

Question

Я использую пакет beanplot , чтобы создать ряд из 16 beanplot. Я использую его немного необычным образом, потому что у меня нет 2 подгрупп на данные, я просто смотрю на 6 групп данных в матрице со всеми остальными группами. Каждую группу следует сравнивать друг с другом - таким образом, 16 бобовых участков.

При построении графика AxB, BxC, CxA и т. Д. Мои значения остаются в одном и том же универсальном масштабе, скажем, от 1 до 100. Однако сами графики плотности немного изменяются для одной и той же группы данных. Например, AxB (например, набор данных Upland) покажет кривую для группы данных A немного шире, чем на графике для CxA (например, набор данных Lowland).

Как мне сохранить согласованность графиков плотности?

Я пытался использовать maxwidth в пакете beanplot, однако это не решает проблему плотности графика.

Вот пример кода:

Данные по нагорью:

plotnum year    actual  type
prf024  2012    14.3    THW
prf031  2012    33.2    THW
prf033  2012    24.8    CU
prf036  2012    40.6    THW
prf040  2012    28.1    THW
prf053  2012    11.7    CU
prf054  2012    12.7    CU
prf075  2012    24.6    CU
prf076  2012    25.5    CU
prf083  2012    38.3    CU
prf159  2012    38.7    THW
prf161  2012    27.9    THW
prf176  2012    49.1    THW
prf178  2012    55.5    THW
prf179  2012    34.2    THW
prf188  2012    43.4    THW
prf189  2012    53.7    THW
prf192  2012    68.6    THW
prf205  2012    2.9 THW
prf209  2012    28.8    CU
prf211  2012    14.4    CU
prf219  2012    18.6    THW

Данные о низменностях:

plotnum year    actual  type
prf024  2012    14.3    THW
prf031  2012    33.2    THW
prf036  2012    40.6    THW
prf040  2012    28.1    THW
prf052  2012    6   CL
prf069  2012    14  CL
prf092  2012    31  CL
prf159  2012    38.7    THW
prf161  2012    27.9    THW
prf176  2012    49.1    THW
prf178  2012    55.5    THW
prf179  2012    34.2    THW
prf188  2012    43.4    THW
prf189  2012    53.7    THW
prf192  2012    68.6    THW
prf205  2012    2.9 THW
prf206  2012    1.3 CL
prf219  2012    18.6    THW

Вот графические параметры, которые я использовал (у меня больше наборов данных, поэтому mfrow равно 16):

{par(mfrow=c(1,16))
par(mai=c(.6,.10,.6,.10))  
par(omi=c(.05,.05,.05,.05) + 0.5) 
plot(0,ylim=c(-20,100), xlab="", ylab="", type="n", bty="n", yaxt="n", axes=FALSE)
}

beanplot(actual ~ type , data = Upland,
        #main = "               *", 
         side = "both", 
         #xlab="Species Group",
         #ylab=expression("Carbon Stocks (t "~ ha^{-1} ~ ")"),
         col = list("#808080", "grey"),
         overallline="median",
          beanlines="median",
          ylim=c(-40,120),
         xlab= "THW+CU",
         #axes= F,
          #maxwidth=.1,
         what=c(1,1,1,1))
mtext(text="*", side=3, outer=FALSE,  cex=1, las=1, font=4, line=-2, at=1)
mtext(text="l", side=1, outer=FALSE,  cex=1, las=1, font=1, line=2, at=1)

beanplot(actual ~ type , data = Lowland,
        #main = "               *", 
         side = "both", 
         #xlab="Species Group",
         #ylab=expression("Carbon Stocks (t "~ ha^{-1} ~ ")"),
         col = list("#808080", "grey"),
         overallline="median",
          beanlines="median",
          ylim=c(-40,120),
         xlab= "THW+CL",
         axes= F,
         what=c(1,1,1,1))
mtext(text="*", side=3, outer=FALSE,  cex=1, las=1, font=4, line=-2, at=1)
mtext(text="l", side=1, outer=FALSE,  cex=1, las=1, font=1, line=2, at=1)

Итак, здесь вы можете видеть, что левый график плотности для первой группы, Upland, немного больше (но в том же масштабе!), Что и для группы Lowland. Отдельные значения, обозначенные маленькими галочками, выровнены и в том же масштабе. Графики плотности должны быть одинаковой формы и размера.

Помощь !!