画曼哈顿图和QQ plot 首推R包“qqman”,简约方便.下面具体介绍以下. 一.画曼哈顿图 install.packages("qqman") library(qqman) 1.准备包含SNP, CHR, BP, P的文件gwasResults(如果没有zscore可以不用管),如下所示: 2.上代码,如下所示: manhattan(gwasResults) 如果觉得不够美观,考虑添加一下参数: manhattan(gwasResults, main = "Manhat

QQplot 横坐标表示的是属性的其中一个测量值1,纵坐标表示另一个测量值2.散点是分位点.点的横纵坐标是这个测量值1和测量值2的分位点的取值. from scipy import stats from matplotlib import pyplot as plt import statsmodels.api as sm import numpy as np # example with the new ProbPlot class #对比两个不同测试值的分布,而不是看测试值满足什么既定分布

Matlab中实现QQ-plot的一个好工具gqqplot 26JUN June 26, 2013 这几天看了一下QQ-plot以及在Matlab中的实现,可是Matlab自带的qqplot函数不能满足我的使用,因此在网上搜索到了一个好工具:gqqplot.gqqplot可以于很多常见的分布进行比较,而qqplot仅仅可以比较正态分布.其zip文件下载链接:gqqplot gqqplot description: GQQPLOT(X,DIST) makes an plot of the quan

QQ-plot深入理解与实现 26JUN June 26, 2013 最近在看关于CSI(Channel State Information)相关的论文,发现论文中用到了QQ-plot.Sigh!我承认我是第一次见到这个名词,异常陌生.维基百科给出了如下定义: "在统计学中,QQ-plot(Q代表分位数Quantile)是一种通过画出分位数来比较两个概率分布的图形方法.首先选定区间长度,点(x,y)对应于第一个分布(x轴)的分位数和第二个分布(y轴)相同的分位数.因此画出的是一条含参数的曲线,参

多因素线性回归 系数由最小二乘法得到 R^2;adjusted R^2:变量变多之后,r^2自然变大,但是这不是反应客观事实,所以引入了adjusted R^2 使用散点图看独立性,也可以使用软件,car package: 任何一个变量显著便使得整个模型(y)显著. 要保证各变量之间相互独立,否则一个变量改变之后另一个变量改变,这两个变量都改变之后y必然改变,但是实际上是第一个变量导致的.所以要检查多元共线性,可使用膨胀系数,相关系数仅考查两个变量之间的关系,而膨胀系数考查一个变量与其他所有变量

转自:http://blog.csdn.net/colddie/article/details/7773278 函数名称 函数说明 调用格式 正态总体的参数检验 ztest 单样本均值的z检验 (总体服从正态分布) [h,sig,ci,zval] = ztest(x,mu0,sigma,alpha,tail) ttest 单样本均值t检验 (总体服从正态分布) [h,sig,ci,tval] = ttest(x,mu0,alpha,tail) ttest2 双样本均值差t检验 (两个总体均服从正

https://zhuanlan.zhihu.com/p/23400450 首发于 R语言数据分析与可视化 关注专栏 登录 写文章     R 语言矩阵散点图 EasyCharts· 15 天前 散点图矩阵因包含多种多副图片,因此含有更多的信息,特别是不同变量之间相关关系.本文介绍了四种绘制散点图矩阵的函数.graphics(pairs),car( spm),GGally(ggscatmat), GGally(ggpairs)四种. ① graphics 包 pairs 函数 data(tree

用Cadence的pad designer制作pad的时候会遇到为thermal relief和anti pad设计尺寸的问题 Thermal relief:正规的中文翻译应该叫做防散热PAD.它主要起一个防止焊接时焊盘散热太快不好焊的作用,在非整层都是铜的情况下它可以做成环形,大小跟Anti pad一样就成了.当在电源层或GND层用时,它的还有减少热冲击的作用,防止焊盘与铜层连接完整的面积过大,因板材与铜皮之间膨胀系数的差异而造成板翘.浮离,或起泡等毛病.这个时候如果就得做成那种镂空的. An

Allegro 总结: 1.防焊层(Solder Mask):又称绿油层,PCB非布线层,用于制成丝网印板,将不需要焊接的地方涂上防焊剂.在防焊层上预留的焊盘大小要比实际的焊盘大一些,其差值一般为10-20mil(这里1mil = 0.0254mm). 2.Anti pad:起一个绝缘的作用,使焊盘和该层铜之间形成一个电气隔离,同时在电路板中证明一下焊盘所占的电气空间.当这个值比焊盘尺寸小时,在负片静态铺铜时焊盘无法避开铜,就会形成短路. 3.Thermal relief:正规的中文翻译应该叫做

1 #f1=scan(file="f1.txt");  2 #f2=scan(file="f2.txt");  3 f=scan(file="5.2_topscreen.txt");                                                                                   4 png(file="myplot2.png", bg="transp

inflate: 胀, 膨, 通货膨胀, 膨胀系数

MEMS器件在射频比如无线通信上有很好的应用.RF MEMS谐振器和诱导器品质因子在微波上有大幅度提高.MEMS开关极大地改进了高频性能和降低了能耗.本篇概要介绍MEMS开关. 自从1979年彼特森(K. E. Petersen)提出微加工继电器以来,许多种类的MEMS开关被开发出来.随着微加工技术的发展,MEMS开关的价格因为材料成本的降低和使用寿命的延长而下降. MEMS开关是用微机械加工技术集成在普通硅片上的开关.它用在射频到毫米波(0.1到1000GHz)的通讯上.与传统半导体器件比如双

关于Tcp封包 很多朋友已经对此作了不少研究,也花费不少心血编写了实现代码和blog文档.当然也充斥着一些各式的评论,自己看了一下,总结一些心得. 首先我们学习一下这些朋友的心得,他们是: http://blog.csdn.net/stamhe/article/details/4569530 http://www.cppblog.com/tx7do/archive/2011/05/04/145699.html //……………… 当然还有太多,很多东西粘来粘区也不知道到底是谁的原作,J 看这些朋友

一般 状态 在人类文明的开始,并探讨了空间和时间的混乱从来没有停止过.马跑得更快.鱼下潜深.鸟振翅高飞,但是,人类并没有很深的不满潜飞不高.为什么?其原因是,马跑得更快,但它不会不知道他们为什么会跑得更快.鱼儿潜的深但它不会想自己为什么会潜的深:鸟儿飞的高但它不会想自己为什么会飞的高,而人类却会思考. 当每天面对着太阳从东方升起从西方落下会想:太阳是不是绕着地球转呢?困惑使人们幻想着一个又一个诗意的神话,同一时候不懈的探寻与求索又使人类跳出一个又一个思想的藩篱向着下一个奥林匹斯山进发.渴望与现实

本文对应<R语言实战>第8章:回归 回归是一个广义的概念,通指那些用一个或多个预测变量(也称自变量或解释变量)来预测响应变量(也称因变量.效标变量或结果变量)的方法.通常,回归分析可以用来挑选与相应变量相关的解释变量,可以描述两者的关系,也可以生成一个等式,通过解释变量来预测响应变量. 回归分析的各种变体 回归类型 用途 简单线性 用一个量化的解释变量预测一个量化的响应变量 多项式 用一个量化的解释变量预测一个量化的响应变量,模型的关系是n阶多项式 多元线性 用两个或多个量化的解释变量预测一个

本文对应<R语言实战>第9章:方差分析:第10章:功效分析 ==================================================================== 方差分析: 回归分析是通过量化的预测变量来预测量化的响应变量,而解释变量里含有名义型或有序型因子变量时,我们关注的重点通常会从预测转向组别差异的分析,这种分析方法就是方差分析(ANOVA).因变量不只一个时,称为多元方差分析(MANOVA).有协变量时,称为协方差分析(ANCOVA)或多元协方差分析

前言 最近迷上了看黄片(大雾).每次总是去搜索想看的片子,什么asian porn anal pussy 什么的我都不知道. 搜索着搜索着我手也累了,而且我喜欢的片子也是有一定的特征的,我不想把所有的时间花费在重复劳动上,于是决定让机器帮我--找出喜欢的片子. (所有代码都在Github上了,文中不放出全部代码了,实在太冗长了) 代码在此:https://github.com/TsingJyujing/xhamster_analysis 我连表结构(CREATE TABLE的SQL语句)都放上去

对于初学R语言的人,最常见的方式是:遇到不会的地方,就跑到论坛上吼一嗓子,然后欣然or悲伤的离去,一直到遇到下一个问题再回来.当然,这不是最好的学习方式,最好的方式是——看书.目前,市面上介绍R语言的书籍很多,中文英文都有.那么,众多书籍中,一个生手应该从哪一本着手呢?入门之后如何才能把自己练就成某个方面的高手呢?相信这是很多人心中的疑问.有这种疑问的人有福了,因为笔者将根据自己的经历总结一下R语言书籍的学习路线图以使Ruser少走些弯路. 本文分为6个部分,分别介绍初级入门,高级入门,绘图与可

R支持4种图形类型: base graphics, grid graphics, lattice graphics,  ggplot2. Base graphics是R的默认图形系统. 一.  基本图形函数plot() plot()命令中的type参数用于明确图形如何绘制,具体type值使用如下: "p" for "points" "l" for "lines" "o" for "overlaid

无论是火电还是核电,将能量转化为电能的方式主要还是烧开水,即加热水产生高压蒸汽驱动汽轮机做功再发电.在进行热力循环分析.流动传热计算时,需获得水和水蒸汽的物性参数.网上主流的水蒸汽物性计算程序是上海成套所的杨宇教授开发的,有Fortran.C.C#.VB等多个语言版本,还有桌面版本,被本专业学生和研究人员广泛使用.可以说杨教授为同行的便利做了很大贡献,本想贴一下他的个人博客的,但发现他的博客链接挂了. 最近打算使用混合编程,写个Python版本的水蒸汽物性计算的接口程序,搞个在线查询的web,然