
笔记
仿生bug
临时记录,过程笔记。【主要是自己忘记了,来这里找找】
展开
-
R语言 方差分析 线性混合模型(简单初探)Linear Mixde Model
使用本地数据 鸢尾花(yuān wěi huā)初探线性混合模型Linear Mixde Model施工ing 交互作用的#载入示例数据data(iris)#载入包library(tidyverse) library(lme4) #生成个体idiris$ID <- rep(1:30,each=5) #假设一朵花测五次iris$ID <- as.factor(iris$ID)#-------------【ANOVA】------------iris_aov1<原创 2022-04-16 10:00:52 · 1615 阅读 · 0 评论 -
R语言 t分布的推导 初级统计学 学生t分布理论
t分布的推导那我们来写写代码,实践这个过程。我设定一个总体均数为0,标准差=1,样本量为1000的人群(图A是这个总体的概率分布)。图B,也是一次抽取三个人,抽了200个,图C一次抽取6个人。都分别计算t值。我们看到图C,一次抽取6个人,分布更多集中在总体均数为0。 还有一个,问题,图B中画绿框的这一部分,怎么解释呢?课后可以想想。(1.抽样误差、2.开端的,无闭合)继续写代码实践,。如一次,,,他们分布会有什么特点。从图片可以得出什么规律:发现一次抽取3次和一次抽取6次比。一次抽取3次的曲线更低,向原创 2022-04-14 16:50:02 · 1454 阅读 · 0 评论 -
R语言 lavaan包进行结构方程模型SEM
学习资料/开发者Github-yrosseel/lavaan原创 2022-03-30 15:51:26 · 2102 阅读 · 0 评论 -
ggplot2 画图参数快速查找,交互式画图,导出为PPT编辑。
随机创建个数据data1<-data.frame(n = rnorm(6),mu = rnorm(6)+1)data1加载包library(ggplot2)第一种 cheatsheet 速查表 图片在下方p=ggplot(data1,aes(x=n, y=mu))+geom_point()p第二种 交互式可以复制控制台下面的 网址到网站打开 library(ggThemeAssist) ggThemeAssistGadget(p)#第三种library(es原创 2021-09-05 11:24:54 · 886 阅读 · 0 评论 -
二代测序16S扩增子测序 双端测序拼接 如何切取碱基质量切割连接 DADA2 简易概念
==(一些思考记录) DADA2切取示例个人控制连续三个低于25%的碱基处,的第一个作为分割点首先要了解自己的测序是可变区的哪几个区,长度多少切割的长度应该要大于可变区的长度,具体大多少没有官方说明最好DADA2进行多次不同位置切割,然后比较table的最小feature前面要去掉引物 引物去掉也去掉Barcode?qiime cutadapt trim-paired \--i-demultiplexed-sequences 01demux/demux.qza \--o-t原创 2021-09-04 17:18:32 · 2038 阅读 · 0 评论 -
R语言markdown的一个 暗黑版本 htm R转网页markdown 结果报告
可以隐藏代码查看结果---title: "c34d的笔记"output: html_document: code_folding: show theme: bg: "#202123" fg: "#B8BCC2" primary: "#EA80FC" base_font: google: Prompt heading_font: google: Proza Libre ver原创 2021-09-04 16:56:25 · 214 阅读 · 0 评论 -
linux服务器 qsub .pbs 提交作业 无特殊报错,但是提交不上的 bug
需要在每次输完代码空一行原创 2021-09-03 18:26:55 · 1906 阅读 · 0 评论 -
R语言 统计描述 分组统计(各种小方案)
#加载包library(tidyverse)#创建一系列数据data=data.frame( row.names = c('QWE','TYE','WUN','EJN','JDB', 'wuy','okf','Aqe','Bdf','wty') , Chinese = c(50, 60, 64, 55, 59, 54, 61, 63, 64, 66 ), Math = c(72, 68, NA, 29, NA, #这里故意设原创 2021-08-23 15:01:57 · 4117 阅读 · 0 评论 -
R语言 trim=0.1 tirm=0.2 什么意思呢?
RT:表示去掉的极端值的比例 下面实例理解首先我编写一组数据,改动第一列中的两位数 99和130> 表示整列数据算平均值> apply(mydata, 2, mean, ) col1 col2 27.3 15.5 > #表示去掉第一列 10*0.1=1个异常值,算平均值,观察第二列没有变化> apply(mydata, 2, mean,trim=0.1) col1 c原创 2021-08-20 18:50:34 · 2135 阅读 · 0 评论 -
R-Studio cheatsheets 速查表 小抄 Github下载
## 下面这个网站直接下载,然后解压找到里面的PDF链接点击下载原创 2021-08-16 11:33:54 · 417 阅读 · 2 评论