python 矩量法_矩量法:β二项分布
我想我已经能够用python准确地编写这个代码:import numpy as np
males = np.array([0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12], dtype=float)
fams = np.array([3, 24, 104, 286, 670, 1033, 1343, 1112, 829, 478, 181, 45, 7], dtype=float)
n = 12
k = 1
m_1 = (sum(fams*males**k))/(sum(fams))
k = 2
m_2 = (sum(fams*males**k))/(sum(fams))
alpha = (n*m_1-m_2)/(n*(m_2/m_1-m_1-1)+m_1)
beta = (n-m_1)*(n-m_2/m_1)/(n*(m_2/m_1-m_1-1)+m_1)
print "n =", n
print "m_1 =", m_1
print "m_2 =", m_2
print "alpha =", alpha
print "beta =", beta
哪些输出:
^{pr2}$
这与示例的结果相同。但是,如果我使用R包VGAM,它使用最大似然,α和β的估计值完全不同x = c(0,1,2,3,4,5,6,7,8,9,10,11,12)
y = c(3,24,104,286,670,1033,1343,1112,829,478,181,45,7)
library("VGAM")
fit=vglm(cbind(x, y) ~ 1, betabinomialff, trace = TRUE)
Coef(fit)
shape1 shape2
0.4241806 4.9069560
我做错什么了吗?在
python 矩量法_矩量法:β二项分布相关推荐
- python ui bs_Guibs的Python学习_列表
Guibs 的 Python学习_列表# 列表# 列表由一系列按特定顺序排列的元素组成, 其中元素和元素之间可以没有任何关系 # 在 Python 中, 用方括号 [] 来表示列表, 并用逗号 , 分 ...
- Python 进阶_生成器 生成器表达式
目录 目录 相关知识点 生成器 生成器 fab 的执行过程 生成器和迭代器的区别 生成器的优势 加强的生成器特性 生成器表达式 生成器表达式样例 小结 相关知识点 Python 进阶_迭代器 & ...
- python基础_格式化输出(%用法和format用法)(转载)
python基础_格式化输出(%用法和format用法) 目录 %用法 format用法 %用法 1.整数的输出 %o -- oct 八进制 %d -- dec 十进制 %x -- hex 十六进制 ...
- python判断题题库大数据技术_智慧树_大数据分析的python基础_搜题公众号
智慧树_大数据分析的python基础_搜题公众号 更多相关问题 社会公众可以查阅烟草专卖行政主管部门的监督检查记录.() 公民.法人或者其他组织不得利用自动售货机销售烟草制品.() 烟草广告中不得有下 ...
- python 切片_全面解读Python高级特性切片
大家好,欢迎来到Crossin的编程教室! 众所周知,我们可以通过索引值(或称下标)来查找序列类型(如字符串.列表.元组-)中的单个元素,那么,如果要获取一个索引区间的元素该怎么办呢? 切片(slic ...
- python爬虫_网易音乐歌单
小白flag7 python爬虫_网易音乐歌单 准备 import os #存放文件处理 import time #程序运行时间计算 import sys #错误信息返回 预留 from seleni ...
- Python基础_第2章_Python运算符与if结构
Python基础_第2章_Python运算符与if结构 文章目录 Python基础_第2章_Python运算符与if结构 Day02之`Python运算符与if结构` 一.昨日回顾 1.回顾昨天的课程 ...
- Python爬虫_案例分析(二)
Python爬虫_案例分析(二) 一.电影天堂案例 import scrapy from scrapy_movie.items import ScrapyMovieItem class MvSpide ...
- python学习_循环语句
python学习_循环语句 第1关:斐波那契数列 斐波那契数列(Fibonacci sequence),又称黄金分割数列. 因数学家莱昂纳多·斐波那契(Leonardoda Fibonacci)以兔子 ...
最新文章
- 进程——Windows核心编程学习手札系列之四
- Java小数中的四舍五入
- 国产服务器芯片对比分析,影速分析:通过国内云主机评测对比国内服务器异同...
- PyTorch框架学习三——张量操作
- 三分法解决凸(凹)函数极值问题
- 在JSP中应用JavaBean
- Flink笔记--深度全面总结
- 系统签名缺少libconscrypt_openjdk_jni.so解决
- Python3.x:定义一个类并且调用
- Android开发环境的安装与配置(两种方法)
- 小刘同学的第一百三十一篇博文
- RabbitMQ-Plugin configuration unchanged
- java毕业论文云笔记_《毕业设计指导的系统设计与实现》论文笔记(八)
- java每日一讲讲什么_每日一讲.ppt
- Mac下的spotlight无法搜索本地资源的解决办法
- Python:flatten()函数用法
- 电脑网线,电脑网线插了没反应怎么回事?如何解决?
- C#怎样解析.CSV文件同时把类转换成JSON
- JavaScript-事件高级
- Catia启动界面图片修改
热门文章
- Vue之设置背景图片
- android之uniapp自定义基座
- 和是java语言字符常量吗_在 JAVA 语言中,下列正确的字符型常量是_电路原理答案_学小易找答案...
- qwidget show 是否有信号_QT中的消息传递与函数回调机制:信号(signal)和槽(slot)...
- macbook快捷键_MacBook 键盘的「fn」键有什么用
- linux服务器知识学习:了解Linux系统的启动过程
- android常用窗口动画,android 自定义dialog,窗口动画,
- python new方法_Python中的__new__()方法的使用
- xilinx set up debug
- java list 赋值jsp,在Struts中使用JavaBean和List(多行数据)类型属性-JSP教程,Java技巧及代码...