By

R语言之前看数据分析的时候想学来着，但是没有操作，最后选择了python。近日帮同学弄一个meta分析。
其实本质上就是利用meta库，求解
（1）加载meta库，镜像源选择“中国上海”
（2）看一下成功了没
（3）数据加载到rate变量
（4）看一下rate成功没
（5）进行换
（6）写入几个变换到csv文件里
（7）进行测试
（8）利用metaprop函数进行求解得到meterate
(9)后处理得到森林图，“Egger’s”图，和Gallbraith图
代码如下：
install.packages(“meta”)
library(meta)
rate<-read.csv(‘rate.csv’,sep=’,’,header=TRUE)
rate
transform(rate,p=event/n,log=log(event/n),logit=log((event/n)/(1-event/n)),arcsin=asin(sqrt(event/(n+1))),darcsin=0.5*(asin(sqrt(event/(n+1)))+asin(sqrt((event+1)/(n+1)))))->rate
write.csv(rate,file=’rate.csv’)
rate
shapiro.test (rate$p)
metarate<-metaprop(event,n,study,data=rate,sm=’PLOGIT’,incr=0.5,allincr=TRUE,addincr=FALSE,tile=”lanjie”)
funnel(metarate)
forest(metarate,digits=2)
metabias(metarate,method=”linreg”,plotit=TRUE,k.min=10)
metainf(metarate)|
表格的整理格式如下

table1

Study Event Num
M. Muzza et al 2014 13 30
F. Wang et al 2013 10 67
HyeYoung Jin et al 2014 20 43
Feng Sun et al 2018 66 110
C. Fu et al 2015 13 45
Kyoung-Jin Park 2015 19 112
Chunyun Fu 2016 13 45
Matsuo et al 2016 11 48
Hong Jiang 2016 10 12
M. Tan et al 2016 60 96
Long et al 2018 49 106
X. Fan et al 2017 21 66

tabel2

表头 表头 表头
内容 内容
内容		 内容
内容 内容 内容

By

前两天看师弟在求一个线性方程组的解，只不过这个方程组的系数是待定的，因此可以算做是利用matlab求函数的解析解。利用同样的方法也可以求数值解。

1.求线性方程组的解析解

syms W1 W2 g1 g2 q1 q2 P0;
syms a A1 A2 b1 b2 c1 c2 d1 d2 e1 e2 G X1 X2 P1 P2 P3;
[WH1, WH2, gH1, gH2, qH1, qH2, PH0]= solve((2-2P1b1/A1)W1==P1e1+2X1-2G+2g1-a+q2-P3,(2-2P2b2/A2)W2==P2e2+2X2-2G+2g2-a+q1-P3, 2c1g1==2P1b1W1/A1+P1e1-P3,2c2g2==2P2b2W2/A2+P2e2-P3,q1==X1-W1-G+g1,q2==X2-W2-G+g2,P0==a-q1-q2,[W1, W2, g1, g2, q1, q2, P0])

已知，a A1 A2 b1 b2 c1 c2 d1 d2 e1 e2 G X1 X2 P1 P2 P3; 并且知道W1 W2 g1 g2 q1 q2 P0和已知量的关系，求解W1 W2 g1 g2 q1 q2 P0;

代码：
syms W1 W2 g1 g2 q1 q2 P0;
syms a A1 A2 b1 b2 c1 c2 d1 d2 e1 e2 G X1 X2 P1 P2 P3;
[WH1, WH2, gH1, gH2, qH1, qH2, PH0]= solve((2-2P1b1/A1)W1==P1e1+2X1-2G+2g1-a+q2-P3,(2-2P2b2/A2)W2==P2e2+2X2-2G+2g2-a+q1-P3, 2c1g1==2P1b1W1/A1+P1e1-P3,2c2g2==2P2b2W2/A2+P2e2-P3,q1==X1-W1-G+g1,q2==X2-W2-G+g2,P0==a-q1-q2,[W1, W2, g1, g2, q1, q2, P0])

主要是利用matlab中的函数slove;

[aa,bb]=slove(关系式子1，关系式子2，[要求的量1，要求的量2])
aa和bb是要求的量1和2的结果

2.利用1中的结果，给出的是方程组，那么[aa,bb]就是方程组的解。具体点说，当方程组是两条线，那么前面给出的就是交点坐标
代码1：
syms x,y;
[aa,bb]=solve(xx+yy==4,x+y=5,[x,y])

代码2：
syms x,y,c,d;
[aa,bb]=solve(xx+yy==c,x+y=d,[x,y])

3类似的还有求微分方程（dslove),具体可以再百度，或者看matlab的帮助文档。

By

写在病情快要好的时间

从手术到现在已经6个月过去了。今天早上去中南检查，范医生说他觉得基本上已经好了，但是校医院的护士说里面还有一个线性小口子，今天看照片，还是有一块，皮都没有长好。不知道这个算是好还是不好，起码这个痊愈肯定是算不上了。

病太烦了，现在还是不能坐着来工作。

要努力克服“一坨二懒三不读书”的问题。

By

今天是出高考成绩的日子，顺便也回忆一下当时我高考成绩的那天，我记得是一个午饭后，爷爷在沙发上躺着，而我在一边又一边地在电脑上刷新页面，查我的高考成绩。最后差不多也是在12点40的时候刷新出来了，考了637分，全省2731或者2173现在已经记不清楚了。当我刷新出的时候，爷爷急得从沙发上坐起来过来看我的成绩。

没有考到自己想要的成绩，从我内心来说，就是这样。

但是想想看，自己也不想复习，所以最后选择了“武汉大学”，然后被调剂到了水利水电学院，最后现在选了农业水利工程，读了博。

然而，比较幸运的是，自己在读研的时候，遇到了井雨，

最近搬东西，又翻到了同学录。看了看，十年了，变化还是真大，但是想一想，也挺快的，朋友们现在都是天南海北，而自己孤悬武汉。平常联系地也少，但是感觉，情谊还是依旧，很抱歉，但是也是客观存在的，一些同学现在几乎很少联系。还真有点“无可奈何花落去的感觉”。

By

今天是妹妹高考成绩出来的日子，学了三年，现在一朝出了成绩。感觉有些失落，毕竟考得不太好。
现在决定要复习了，鲸鱼提了一下要求：

（1）是要自己选择，自己愿意
（2）自己要有学的意愿，不要太大的压力
我现在刚刚又和她讲了一下，要把学习的方法思路，以及基本功，还有基础知识。

小时候，觉得妹妹很辛苦，每天都跟着妈妈在县里。得到的关爱比较少。现在长大了，希望高考考好一点，有个好的归宿。

By

今天下午和去小米授权店换电池。说是主板漏电，平时充电也有问题。最后觉得对方说的可能不对，然后就业没换主板。只换了电池。

完了后和鲸鱼去吃巴厘龙虾，很是不错。