《随机效应模型理论与应用ppt课件.pptx》由会员分享,可在线阅读,更多相关《随机效应模型理论与应用ppt课件.pptx(16页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、我吓了一跳,蝎子是多么丑恶和恐怖的东西,为什么把它放在这样一个美丽的世界里呢?但是我也感到愉快,证实我的猜测没有错:表里边有一个活的生物随机效应模型理论及应用我吓了一跳,蝎子是多么丑恶和恐怖的东西,为什么把它放在这样一个美丽的世界里呢?但是我也感到愉快,证实我的猜测没有错:表里边有一个活的生物因子(Factor):影响研究变量的各个原因水平(Level):一个因子的不同状态效应(Effect):因子的各水平对所研究变量的影响1我吓了一跳,蝎子是多么丑恶和恐怖的东西,为什么把它放在这样一个美丽的世界里呢?但是我也感到愉快,证实我的猜测没有错:表里边有一个活的生物Example 1 : Respe
2、ct to simple linear regression model,advantages of random effects model and fixed effects model. n 个关于商品销售价格和销售量的数据点: (xk ,yk), k=1,n 1, , kkkyxkn 这些商品可以分为M 种。但是,传统的线性模型,忽视了商品的种别差异,即:2我吓了一跳,蝎子是多么丑恶和恐怖的东西,为什么把它放在这样一个美丽的世界里呢?但是我也感到愉快,证实我的猜测没有错:表里边有一个活的生物 但事实上,不同种商品之间存在很大的差异。按照不同商品将上述数据分组,记为(xij,yij),i
3、 =1,M, j =1,ni 。 (xij,yij)表示第i 种商品的第j 个数据点。将同种商品的价格和销售量的观测点连起来,大部分商品的变化趋势和简单线性回归得到的曲线并不相同。3我吓了一跳,蝎子是多么丑恶和恐怖的东西,为什么把它放在这样一个美丽的世界里呢?但是我也感到愉快,证实我的猜测没有错:表里边有一个活的生物假设每种商品都有自己特定的销售模型:假设截距i 是随机的,且可以表示为:这里 是总体的平均销售量,bi 是随机效应,反映特定商品销售与总体的商品平均销售之间的偏差,合并两式得:2(0,)ibN组间误差组内误差42220.99 1, 1,ijiijijiyxiMjn 1,iibiM
4、1, 1,ijiijijiybxiMjn我吓了一跳,蝎子是多么丑恶和恐怖的东西,为什么把它放在这样一个美丽的世界里呢?但是我也感到愉快,证实我的猜测没有错:表里边有一个活的生物Fixed effects or Random effects Fixed effects model or Random effects model:Cannot explain( ) 1, 1,ijiijyfiajbIndependent variable(单向分类模型)5我吓了一跳,蝎子是多么丑恶和恐怖的东西,为什么把它放在这样一个美丽的世界里呢?但是我也感到愉快,证实我的猜测没有错:表里边有一个活的生物( ) 1
5、, 1,ijijijyfiajb(双向分类模型)Random effects Fixed effects(混合效应模型)6我吓了一跳,蝎子是多么丑恶和恐怖的东西,为什么把它放在这样一个美丽的世界里呢?但是我也感到愉快,证实我的猜测没有错:表里边有一个活的生物Specific levels could be replaced by other levels : the chosen levels are arbitraryarbitrary or substitutablesubstitutable(水平是否随机)(水平是否随机)The conclusion of the experiment
6、is to be generalizedgeneralized(结论是否通用)(结论是否通用)Conclusions drawn for each separate level are not of interestnot of interest:not the specific items are important but the population theyare drawn from(水平是否感兴趣)(水平是否感兴趣)When should a factor be considered random?random?7我吓了一跳,蝎子是多么丑恶和恐怖的东西,为什么把它放在这样一个美丽的
7、世界里呢?但是我也感到愉快,证实我的猜测没有错:表里边有一个活的生物Example 2 : 假定农业基地引进a 个小麦品种进行试验。这时,小麦品种是我们主要感兴趣的因子。土质肥沃的程度基本一样的一大块田分成n块,其中ni 块用来种植第i 种的小麦,于是, 。yij 表示第i 种的小麦在它的第j 块田上的产量:1aiinnFixed effects term( ) 1, 1,ijiijiyfiajn8我吓了一跳,蝎子是多么丑恶和恐怖的东西,为什么把它放在这样一个美丽的世界里呢?但是我也感到愉快,证实我的猜测没有错:表里边有一个活的生物Example 3 : 通过b 个地震台站搜集到了a 次地震的
8、地震动数据,yij 表示第j个地震台站记录到的第i 次地震事件的地震动(PGA、Sa(T ).)( ) 1, 1,ijiijyfiajbRandom effects term9我吓了一跳,蝎子是多么丑恶和恐怖的东西,为什么把它放在这样一个美丽的世界里呢?但是我也感到愉快,证实我的猜测没有错:表里边有一个活的生物随机效应模型的极大似然估计法(、 2、 2):Step 1. Estimate the model parameter values,by linear regression procedure.Step 2. Given , estimate 2 and 2 by maximizing
9、 the likelihood given inequation blow (by simplex method). 222122222111111lnln(2 )()ln()ln()222()11 ()22iMiinMMiiiijiijiiLNNMnn YuyYn1111iiniijjiniijjiYynuunwhereObserved valuePredicted value10我吓了一跳,蝎子是多么丑恶和恐怖的东西,为什么把它放在这样一个美丽的世界里呢?但是我也感到愉快,证实我的猜测没有错:表里边有一个活的生物Step 3. Given , 2 and 2,estimate Yi by
10、equation blow.2122inijijjiiyunStep 4. Given i ,estimate new by using linear regression procedure for (ln yij - i).Step 5. Repeat steps 2,3,and 4 until the likelihood in step 2 is maximized.11我吓了一跳,蝎子是多么丑恶和恐怖的东西,为什么把它放在这样一个美丽的世界里呢?但是我也感到愉快,证实我的猜测没有错:表里边有一个活的生物在进行反应谱谱值或谱比的一系列随机效应模型回归时的问题初始回归之后的谱比设计反应谱
11、12我吓了一跳,蝎子是多么丑恶和恐怖的东西,为什么把它放在这样一个美丽的世界里呢?但是我也感到愉快,证实我的猜测没有错:表里边有一个活的生物112233 1, 1,ijiijiyxxxiajn13我吓了一跳,蝎子是多么丑恶和恐怖的东西,为什么把它放在这样一个美丽的世界里呢?但是我也感到愉快,证实我的猜测没有错:表里边有一个活的生物固定参数1:112233(new) 1, 1,ijijiijiyyxxxiajn 如此一个一个参数固定,直至剩最后一次参数,进行最后一次回归。但是,这很难保证一次就得到光滑的曲线,需要反复调节。14固定参数2:.我吓了一跳,蝎子是多么丑恶和恐怖的东西,为什么把它放在这样一个美丽的世界里呢?但是我也感到愉快,证实我的猜测没有错:表里边有一个活的生物谢谢