c1684517745的博客

原创 python 模拟一个果园(给定规则的随机事件,最基本的程序设计方法学)

这是我初学python在咸鱼接的一个单子,所用的代码都是最简单入门的,大佬勿喷,当个乐子看,主要是一个思维方面的问题这应该是他们的一个作业,我当时学python才3个月作业,在咸鱼接单子,当时想着是写不出来不要钱,急单就不接,多写一点就会使自己的代码水平增长,当时抱着这个想法来着这篇文章就是基于那天的一个咸鱼的单子,要求模拟一个果园,然后计算若干年后的产量其实当时看到这里的时候是懵逼的,我以为要用一些什么高级的库才能完成这个事情,然后我问他,规则呢,他说,规则你自己定,我一寻思自己定,那...

原创 simplex法(单纯形) 并在python实现简单的应用

我最近在学习动物遗传育种GBLUP的一个算法-dfreml一般做gwas,这些传统软件如asreml gmat(我老师写的,有需要可以联系) dmu 等等,计算方差组分都是运用求导 用EM和AI迭代算法求出有没有一种不用求导的算法计算方差组分呢,那就是非求导的df-remlDf-reml基本包括两步第一步即为构建似然函数第二步为搜索未知值使似然函数极大化-simplex法simplex法的主要步骤Simplex法求似然函数的最大值假设有p个需要估计的参数: ...

原创 使用python发生此类错误 ValueError: endog must be in the unit interval.