beta.int.R 里写到:
area<-function(r){
s<-3.14*r^2
return(s)
}
退到R的界面：先source这个函数文件:
>source(beta.int.R)
>r<-2
>area(r) 不source是无法用area这个函数的
[1] 12.56
1 普通堆叠图
2 拆成柱形图
3 添加冲积图
安装依赖：
绘制冲积图：
4 添加facet_wrap分面
5 添加geom_segment分组标记
数据准备：准备geom_segment需要的x、x_end值
绘图：
6 翻转90度
7 绘制堆叠面积图
数据准备：给每个样品按数字编号
绘图：
这配色似乎还可以
一组好看的堆叠图参数：
参考：
R语言ggplot2绘制分组箱型图和分组柱状图
Make Grouped Boxplots with ggplot2
“面积”的意思是积分吗？还是说x轴之下的部分也算正的面积？x必须是正的吗？
我就当x必须是正的了哈。
只考虑积分的话
F<-function(x) integrate(sin,0,x)$value-1.5uniroot(F,c(0,pi))$root
不然就是
F<-function(x) integrate(function(t) abs(sin(t)),0,x)$value-1.5uniroot(F,c(0,20))$root #20应该够大了
（虽然结果都一样o(╯□╰)o）
