在教师节收到学生提问,刷我B站74小时视频的时候看到我演示了RNA-seq差异分析只用了一行代码就完成了3大R包的全部分析,并且输出了对应的图表结果,觉得很神奇,但是B站视频并没有配套讲义和代码还有测试数据...## 表达矩阵来自于R包: airway if(F){ library(airway) data(airway) exprSet=assay(airway) group_list=colData...下面的图表是如何自动出来的呢? ? 因为这个 run_DEG_RNAseq 函数的代码非常长,这里我就不贴在公众号了哈,大家可以在我的GitHub的GEO项目找到它!...这个时候是没有标准答案的,因为每个R包都非常热门,引用量都是好几千,你选择哪个都符合市场规律,不过,我这里有一个代码,对3个结果根据阈值筛选交集。...当然是啊,都会写代码了,还有什么是不能为所欲为的呢? 同样的,代码也是在GitHub,需要你仔细理解,不过我有一个小小的要求,请不要把我的代码雪藏,或者刻意隐瞒。
p=14854 一种类型的平滑称为样条平滑。柔性金属(通常是铅),可以用作绘制平滑曲线的参考。将选择一组点(称为结),然后将样条线压在特定的x,y点,然后弯曲以通过下一个点,依此类推。...由于金属的柔韧性,此过程将生成通过这些点的平滑曲线。 在数学上,可以通过选择结点并使用(通常是三次)回归来估计结之间的点,并使用演算来确保每条单独的回归线连接在一起时曲线都平滑,从而重现该过程。...平滑的程度由参数控制,通常在0和1之间的范围内。 为了说明,我们考虑由来自1910至2004年的小麦产量数据集 。 生成数据图,并叠加样条曲线平滑度。...在R中,可通过 supsmu 函数获得 超级平滑器。 为了说明这一点,考虑汽车数据。以下几行产生了 重量 与 MPG的关系图,并叠加了一条超平滑线。...6.使用SAS,Stata,HLM,R,SPSS和Mplus的分层线性模型HLM 7.R语言中的岭回归、套索回归、主成分回归:线性模型选择和正则化 8.R语言用线性回归模型预测空气质量臭氧数据 9.R
前言 又是一年考研季,今天有朋友问我去年的学校录取情况。作为一只“菜狗”,这那是我该操心的事,但我还是去官网找了下去年的结果。...具体来说,主要包括如下技术要点: Stata16 调用 Python 下载 PDF,转化为 Excel 表格; Stata 爬取学校官网的学院名单,整理和生成学院代码; Stata 中的 putdocx...在开始之前,要注意使用的是 Stata 16,而且需要先安装 Python 的 pdfplumber 库;Stata 的 tbl2putdocx.ado 文件,这是来自 statalist[1] 解决...在公众号后台回复【zuledocx】可以获取我整理好的 ado 文件,之后放入 Stata 所在路径的 ./ado/plus/t 文件夹中即可。...在制作学院代码过程中,一是需要剔除不招硕士的学院,二是要和名单中的专业进行比照,确保拟录取专业和学院相对应。 copy "http://www.zuel.edu.cn/schools/" ".
Stata较好地实现了使用简便和功能强大两者的结合。尽管其简单易学,它在数据管理和许多前沿统计方法中的功能还是非常强大的。...如果在一个程序中出现一个错误,找到并改正这个错误将是困难的。在所有的统计软件中,SAS有最强大的绘图工具,由SAS/Graph模块提供。...R语言:科研首选 R是用于统计分析、绘图的语言和操作环境,属于GUN系统的一个自由、免费、源代码开放的软件,它是一个用于统计计算和统计制图的优秀工具,多用于论文,科研领域。...因此R有很多最新的模型和检验方法,但是非常难自学,对英语的要求很高。R与SAS的区别在于,R是开放免费的,处理更灵活,同时对编程要求较高。 界面展示: ?...如果你是SEM新手,建议使用Amos,因为有防呆装置,一般不容易犯错,确定是知其然却不知其所以然。如果你对数学很有兴趣尤其对矩阵概念好的人,建议使用Lisrel,上手虽麻烦,但功能强大易操作。
所以,我的理解是:有其他语言的基础,学习新语言时,应该迁移和对比相结合。具体而言:对比着学,尽快找到感觉;注重工具之间的特性;问题导向,根据需求选择合适的工具。...本文以 Stata 自带 auto.dta (1978年美国汽车数据) 数据为例,对照着 Stata 的完成多元线性回归的过程,展示在 Python 中如何跑回归。...Stata 中运行回归 在 Stata 中,完成整个实证的过程大致如下: cd .....接下来,将在 Python 中按照此流程重现。...本文演示的还仅是最简单的多元线性回归,一些复杂和前沿的计量模型, Python 中可能还没有现成的包,需要自己编写代码。
本文的学习方法是在追求实用性和全面性的不断平衡中建立起来的。这里有许多优质的 R 语言免费资源,但可惜并非全部。...它源于学界,但今天你会在越来越多的商业环境中看到 R 语言的身影,它现在成了商业软件公司如 SAS,STATA 和 SPSS 的贡献者。...步骤 3:R 语言的核心 - >packages 每个 R 语言包只是一组用于特定目的代码,旨在被其他开发人员重复使用。除了主代码库之外,程序包通常还包括数据、文档和测试。.../) haven package 可以将 SAS,STATA 和 SPSS 数据文件导入 R 语言中。...有很多软件包资源可供学习,比如,互动代码教程,是 Hadley Wickham 做的一个备忘清单和即将出版的新书。
pwd=ny0u 提取码:ny0u –来自百度网盘超级会员V4的分享 正式安装 1.下载Stata 17安装包 2.双击Stata17.exe开启安装 3.接受同意条款,然后继续安装 4.选择想要安装的版本...6.等待软件安装完成,先不要打开软件;安装过程中不要进行其他操作,耐心等待几分钟~ 7.复制软件包中的另一个程序:StataMP-64.exe 8.切换到你的安装路径.安装路径就是你存放所有Stata安装程序和文件的地方...结尾 在本篇博客中,我们详细介绍了如何安装最新版本的Stata17。无论您是新手还是有经验的用户,我们提供的步骤和注意事项都将帮助您顺利完成安装过程。...Stata17不仅带来了许多新功能和改进,还提升了用户体验,使数据分析变得更加高效、准确和愉快。现在,您可以开始探索Stata17的各种功能,将其应用于您的研究、工作或学习中。...祝您使用Stata17愉快,愿本文对您有所帮助!如果您有任何疑问或反馈,请随时联系我。
下图是 Stata 界面中工具栏的“图形”所包含的内容(图1)。 Stata Window_Graphics .png 2....要说明的是,由于绘图命令十分“庞大”,在学习和应用中,不断积累各方资料中的图形代码很有必要;同时在绘图中也要善用 Graph Editor 对图形进行局部细节的优化,毕竟我们不可能记得所有绘图命令的选项...graph_structure.jpg Stata 的绘图代码主要包括四个部分:(1)命令(Graph Commands);(2)选项(Options);(3)风格(Styles);...下图是 twoway 家族的所有成员(图3),图2只展示了部分常用的图形类型。" [ ] " 表示代码中可以省略的部分。虽然可以省略,但这部分却是掌握绘图命令的核心。...下面,我们用一组容易混淆的示例进行Stata绘图的展示,以帮助我们了解 Stata 的功能与绘图类型。 3.
该方法被我用于19年发表的文章(见参考文献)中,不过当时我是基于metafor包用几行代码试了一下,然后用邮件咨询文章作者,他使用的是Stata,开发了一个叫ipdmetan的包(比我厉害多了~),最后也是使用...stata画的图放在文章里。...因为Stata的结果和我使用几行代码的结果一致,所以我在去年7,8月左右想实现一个R版本,在Github创建了相应的仓库。...然后我跟他说我这个包根本就没写,已经弃疗一年了,叫他去找Stata包,然后附送了之前使用的R代码。因为他是一个R用户,所以回信跟我说有个R包就好了。...于是我就简单粗糙地造了一个,就有了真正可以用的metawho package,并把我文章的数据放在里面重复了一下结果。
刚有朋友问我怎么调整 boxplot 中 x 轴标签,用上图重现了他的问题。换句话说,问题是如何解决 graph 中 x 轴重叠的问题。...简单解释代码中的 {char 0xa0},这是 SMCL(Stata Markup and Control Language) 语言,是 Stata 中的一种标记语言, help 文档大都是用它编写的。...char 表示字符,0xa0 是 ASCII 编码的空格。所以这里实际是“偷梁换柱”,将 x 轴刻度每隔 5 个单位换为空格,这样绘制出来的图就实现了肉眼不可见的空白。...我刚开始也顺着这个思路考虑是否能通过 SMCL 语句更改 x 轴刻度的倾斜角度,但 SMCL 似乎没有并不能实现文本倾斜。对绘图中可用的 SMCL 语句,可自行 help text 查看。...附录:演示数据生成代码 clear set obs 126 egen year = seq(), block(6) from(1990) to(2010) set seed 2803 gen y = exp
STATA与EVIEWS都是我偏好的计量软件。...关于实验 在随机实验中,样本被随机分成两组,一组经历处理条件(进入干预组),另一组接受控制条件(进入对照组),然后比较两组样本的效果指标均值是否有差 异。...关于拟合优度、变量选择原则及估计值绝对大小的意义 在人人的“数据分析”小站中,某同学提出这样一个问题:“多元回归分析中,怎么选择自变量和因变量,可以使R方达到80%以上?”...,根本没必要追求R方到80%的程度,一般来说,有个20%、30%就非常大了。...如果一定要增大R方,那么最应该做的的确是对纳入模型的变量进行选择。选择纳入模型的原则我认为有三条。
SAS也是各种算法模块都有,这一点跟上跟python,r不会差多少,但sas的代码写法跟其他软件不太一样,用多了python,r去写那个还是感觉怪怪的。...sas可以对一个命令逐行去运行,不需要自己加循环,可以少写一点代码。SAS有一些函数的默认跟其他软件不太一致,比如merge函数,其他软件一般是默认内连接,只保留匹配到的,但sas默认是外连接的。...SAS有一个非常好的特性是提取数据库里的数据,create table生成的数据是存在硬盘里的,而python,r等等其他软件提到的数据都是存在内存里的,这样就没法用python和r提太多量的数据出来,...个人理解,excel中写vba尽量减少激活各个sheet单元格,跳来跳去,有什么东西可以直接往里写,这样会快一些。...STATA和EVIEWS是两个功能非常类似的软件,计量上用的比较多,我当时选择了stata,没有学EVIEWS,整体来说,都是计量软件,STATA对面板数据支持比较多,EVIEWS对时间序列支持更好,STATA
我认为 A 在我的检验中可能不成立,因为 C 或者,D 是判断 A 是否成立的条件; 又如,D 是另外一种计量方法但是并没有 A 这个假设前提....前期相关推文 专题:专题课程 公开重现资料时如何发布涉密数据 直播-我的甲壳虫-论文精讲与重现 专题:论文写作 连享会:论文重现网站大全 论文重现:外部命令的版本控制 Stata:论文重现代码模板 Stata-JPE...专题:Stata入门 Stata的版本兼容性问题:可重复研究 专题:Stata资源 可重复研究:Stata一马当先 会计期刊论文的结果可重现吗?...专题:数据处理 Stata:变量非重复值统计-distinct Stata结果重现:dependencies命令-外部命令的版本控制 Stata数据处理:统计组内非重复值个数 专题:结果输出 Stata...Stata:一行代码实现安慰剂检验-permute acreg:允许干扰项随意相关的稳健性标准误 aoeplacebo:地理安慰剂检验 专题:倍分法DID 多期DID之安慰剂检验、平行趋势检验 专题:内生性
但是由于,我的SPPS上个月删掉了,占用1.5g内存,而且没有破解。这次,我用最不怎么熟悉的Stata来做主成分分析和因子分析。...今天用Stata,点点就ok了 ,敲什么代码,再说了我又不会敲。 ? 在这里插入图片描述 ?...在这里插入图片描述 排名第一的就是净利润,企业就是要赚钱才有竞争力。这不是常识了,分析和没分析差不多,我还是干Spring,Django吧。 看看上面内容,没有代码怎么行?...虽说我专业都快挂完了,心痛一秒钟。 下面是百度百科给因子分析模型定义的,我抄了下。 的一种统计方法,是一种降维技术. 做因子分析的前提是自变量之间有相关关系....在这里插入图片描述 初始因子和Stata的结果一样 ? 在这里插入图片描述 在Stata中我们没有旋转变换, ? 在这里插入图片描述 旋转变换的后的 ? ? 答案是柳州化工,我听说柳州螺蛳粉,五菱。
这太简单了,让新手退出我vim就行了。 Linux老大:为啥? vim: 哈哈,因为新手不知道怎么才能退出vim,就会瞎胡乱按一通, 非常随机,这不就形成随机数了吗?!...不能预测,不能根据已经生成的随机数,推测出下一个随机数是啥 3. 不能重现, 无法重现和某一随机数列完全相同的数列 听到此处,大家都吸了一口冷气,这要求够高的!...人类通过掷骰子可以达到这个要求,但是计算机里都是确定的算法和程序,这该怎么办? C老头儿说:我提一个方案,我听说人类有个算法,叫做什么线性同余算法,似乎可以生成随机数。...但是,如果我知道了种子,那就可以生成和你一模一样的随机数列,所以不满足‘不可重现’的性质。” 看来生成真正的随机数太难了,大家都沉默了。...CPU阿甘马上接口:“然后,可以用个Hash算法对这个池子中的内容做个消息摘要,结果就是真随机数了!杂乱无章,无法预测,无法重现。” vim感觉有点不爽,这俩人也太会抢功劳了。
STATA与EVIEWS都是我偏好的计量软件。...在随机实验中,样本被随机分成两组,一组经历处理条件(进入干预组),另一组接受控制条件(进入对照组),然后比较两组样本的效果指标均值是否有差异。...关于拟合优度、变量选择原则及估计值绝对大小的意义。 在人人的“数据分析”小站中,某同学提出这样一个问题:“多元回归分析中,怎么选择自变量和因变量,可以使R方达到80%以上?”...根本没必要追求R方到80%的程度,一般来说,有个20%、30%就非常大了。...如果一定要增大R方,那么最应该做的的确是对纳入模型的变量进行选择。选择纳入模型的原则我认为有三条。
声明一下,我是在 Windows 10 (64-bit) + Stata/MP 16.0 + Python 3.7 环境下运行的代码,在不同系统和版本的程序下运行可能不太顺畅,如果你用自己的环境玩得不太顺畅...最后一种方式值得探索,winsound 中的 PlaySound 函数可以播放音乐,想着可以结合爬虫,爬取排行榜和下载音乐,最终用 Stata 封装一个播放最新排行榜音乐的命令。...我已经将过程中需要的资料打包备索,公众号后台回复【Beep】即可获取本文涉及的安装包、代码和其他文件。...\Music\ding.wav',winsound.SND_FILENAME) end 小结 以上就是我知道让 Stata 发出提示音的几种方法?你还有其他更好的方法吗?可以后台交流。...附录:配置环境变量 在设置中搜索环境变量,之后在编辑账户环境变量中找到 path 点击编辑-新建,然后粘贴进 Beep 所在的文件夹,点击确定。
STATA与EVIEWS都是我偏好的计量软件。...关于实验 在随机实验中,样本被随机分成两组,一组经历处理条件(进入干预组),另一组接受控制条件(进入对照组),然后比较两组样本的效果指标均值是否有差异。...关于拟合优度、变量选择原则及估计值绝对大小的意义 在人人的“数据分析”小站中,某同学提出这样一个问题:“多元回归分析中,怎么选择自变量和因变量,可以使R方达到80%以上?”...根本没必要追求R方到80%的程度,一般来说,有个20%、30%就非常大了。...如果一定要增大R方,那么最应该做的的确是对纳入模型的变量进行选择。选择纳入模型的原则我认为有三条。
计量经济学也有很多小的门类,请对号入座。 有很多软件,Stata, Matlab, R, Sas是相对来说用的比较多的。...但是我一般不做这方面,也不是很有发言权。 如果做理论计量,stata eviews是没有现成的包的,而且即便Stata可以编程,可编程能力也是很差的,而且不稳健。所以懂R和Matlab就非常顺手。...如果对速度要求高,特别是金融计量很多对速度有要求的,可以考虑C、Fortran等语言。C和Fortran肯定是最快的。还有一个叫做OX的,速度快,但是也很小众。...最后想起一句话来,关于这些软件的选择(stata除外,因为stata在应用计量领域的地位是几乎不可替代的)可以用两句话来概括:如果你自己的时间比计算机的时间宝贵的话,学习matlab、R、甚至Python...我觉着stata、matlab是比较推荐的,一个方便,一个灵活,都很强大,而且学习难度都不大,用的人都很多,交流起来相当方便。 PS:各个软件的速度比较, 可以看Julia的主页。
领取专属 10元无门槛券
手把手带您无忧上云