關(guān)于乳腺癌術(shù)后生存情況與患者年齡的研究生存分析結(jié)課論文

《生存分析結(jié)課論文》——關(guān)于乳腺癌術(shù)后生存情況與患者年齡的研究 班級(jí)：姓名：學(xué)號(hào)：

目錄

摘要本文討論45歲以上乳腺癌患者的術(shù)后生存狀況。對(duì)44名45歲以上的乳腺癌患者的資料進(jìn)行回顧性分析，按年齡分為兩組，其中A組（<50歲,25例），B組（≥50歲，19例），探討乳腺癌患者術(shù)后生存情況與患者年齡間的關(guān)系。結(jié)果有統(tǒng)計(jì)學(xué)意義(P<0．01)。年齡是乳腺癌的一個(gè)獨(dú)立預(yù)后變量，但乳腺癌的其他影響患者生存狀況因素如：臨床分期、淋巴結(jié)轉(zhuǎn)移、病理類(lèi)型、手術(shù)方式對(duì)乳腺癌患者的影響也是不容忽視的。關(guān)鍵詞生存分析乳腺癌年齡Kaplan-Meier估計(jì)Nelson-Aalen估計(jì)Cox模型

1.問(wèn)題的提出乳腺癌是女性最常見(jiàn)的惡性腫瘤之一。且發(fā)病率呈逐年上升的趨勢(shì)，在歐美國(guó)家，乳腺癌占女性惡性腫瘤的25%-30%.乳腺癌常發(fā)病于停經(jīng)婦女，我國(guó)則常見(jiàn)于絕經(jīng)前婦女，45—50歲發(fā)病率較高。中老年婦女是乳腺癌發(fā)病的主要對(duì)象。發(fā)病年齡較歐美國(guó)家年輕10歲左右。由文獻(xiàn)報(bào)道年齡是一個(gè)對(duì)復(fù)發(fā)率有影響的獨(dú)立因素，年齡在45-50歲的患者復(fù)發(fā)率增加，為比較不同年齡乳腺癌術(shù)后生存狀況的差別。本文從生存狀況變化的角度做生存性分析，探討乳腺癌術(shù)后生存情況與患者間年齡關(guān)系。2.數(shù)據(jù)的來(lái)源選取患乳腺癌的44名婦女，初治均為手術(shù)治療，分為兩組。A組為年齡在45歲到50歲的患者，B組為年齡在50歲以上的患者。5年后得到下列復(fù)發(fā)時(shí)間。時(shí)間（月）數(shù)據(jù)來(lái)源于《生存數(shù)據(jù)分析的統(tǒng)計(jì)方法》A組45916121310232829313247414157627410013920+258+269+B組8101012142048707599105162169195220161+199+217+245+3.模型方法介紹和總結(jié)3.1生存時(shí)間函數(shù)描述生存時(shí)間分布規(guī)律的函數(shù)主要有生存函數(shù)、死亡概率函數(shù)、概率密度函數(shù)和危險(xiǎn)率函數(shù)。為了后文敘述方便，這里主要介紹生存函數(shù)和危險(xiǎn)率函數(shù)。3.1.1生存函數(shù)在描述生存規(guī)律的數(shù)量指標(biāo)中，以往常用的指標(biāo)是某個(gè)特定時(shí)間的生存率（例如：3年生存率、5年生存率）。這一指標(biāo)的主要缺陷為不能反映整個(gè)生存規(guī)律，一個(gè)理想的指標(biāo)應(yīng)該是任意時(shí)間的生存率，即生存率是任意時(shí)刻t的函數(shù)。其意義是研究個(gè)體生存時(shí)間長(zhǎng)于t的概率。若令T為生存期，s(t)為任意時(shí)刻t的生存率，得 （3.1）則稱為生存率函數(shù)，簡(jiǎn)稱生存函數(shù)。從T的分布函數(shù)可知，（3.2）將對(duì)在直角坐標(biāo)系作圖畫(huà)出生存率曲線，陡峭的生存曲線表示低的生存率或的生存時(shí)間。另外，從圖上可粗略估計(jì)出中位生存時(shí)間，即生存曲線上取生存率為50%時(shí)所對(duì)應(yīng)的時(shí)間。3.1.2危險(xiǎn)率函數(shù)（HazardFunction）如果我們考慮已活到t時(shí)刻的患者，在時(shí)間t附近的瞬間死亡危險(xiǎn)性，根據(jù)數(shù)學(xué)上極限性質(zhì)，可表示為：h(t)= (3.3)則稱h(t)為危險(xiǎn)率函數(shù)。即相當(dāng)于條件瞬間死亡率。3.1.3S(T)與h(t)的關(guān)系其中，稱H(t)為累計(jì)危險(xiǎn)率函數(shù)。生存函數(shù)和危險(xiǎn)率函數(shù)在不同的生存時(shí)間分布中（例如Weibull、loglogistic等）有著特定的函數(shù)形式。3.2.1生存函數(shù)的估計(jì)非參數(shù)法對(duì)生存數(shù)據(jù)的分布型沒(méi)有相應(yīng)的要求，因而適用面比較廣。醫(yī)學(xué)研究中，大量的生存數(shù)據(jù)其分布是不規(guī)則、不確定或未知分布，因而常用非參數(shù)法估計(jì)生存率。根據(jù)樣本含量的大小可分別選擇壽命表法或乘積限估計(jì)法。Nelson-Aalen估計(jì)在有刪失的情況下，可以根據(jù)累積死亡率與生存函數(shù)的關(guān)系來(lái)估計(jì)累積死亡力函數(shù)H(t)。這時(shí)估計(jì)式為：。另外有一個(gè)累積死亡力估計(jì)式，它與以乘積限估計(jì)式為基礎(chǔ)的估計(jì)式相比，具有更好的小樣本性質(zhì)，這一估計(jì)式由Nelson建議，然后由Aalen重新發(fā)現(xiàn)并加以改進(jìn)，這就是Nelson-Aalen估計(jì)式，即在最大的時(shí)間觀察范圍內(nèi)的定義如下：該估計(jì)式的方差可以從下式得到：以累積死亡率的Nelson-Aalen估計(jì)式為基礎(chǔ)，生存函數(shù)的另一個(gè)估計(jì)式為：。Nelson-Aalen估計(jì)式在分析數(shù)據(jù)時(shí)主要有以下兩方面的應(yīng)用，其一是在選擇事件發(fā)生時(shí)間的參數(shù)模型方面的應(yīng)用，其二是為死亡率h(t)提供粗估計(jì)，這些估計(jì)值是Nelson-Aalen估計(jì)式的斜率。乘積限估計(jì)法當(dāng)數(shù)據(jù)個(gè)體較少時(shí)，為充分利用每個(gè)數(shù)據(jù)的信息，必須采用更為精確的估計(jì)方法。這些估計(jì)方法中應(yīng)用最多、效率較高的是Kaplan-Meier在1958年提出的乘積限估計(jì)(Product-limitestimator)。因而此法又稱Kaplan-Meier法。乘積極限法適用于離散數(shù)據(jù)，它用于建立時(shí)刻t上的生存函數(shù)。它的原理是根據(jù)時(shí)刻及其之前各時(shí)間點(diǎn)上的條件生存率的乘積，來(lái)估計(jì)時(shí)刻的生存函數(shù)和它的標(biāo)準(zhǔn)誤。設(shè)代表個(gè)觀察對(duì)象的生存時(shí)間，設(shè)為時(shí)刻之前生存的個(gè)體數(shù)目，即危險(xiǎn)集的大小，再設(shè)表示生存時(shí)間的截尾性質(zhì)，。又令表示觀察對(duì)象在時(shí)刻的條件生存率，即對(duì)于，有：，其中那么，觀察對(duì)象在時(shí)刻時(shí)的條件死亡率如下：對(duì)于，Kaplan-Meier法定義時(shí)刻上的生存函數(shù)和它的標(biāo)準(zhǔn)誤的估計(jì)公式如下：該法的基本思想與壽命表法基本相同，所不同的是將生存時(shí)間（包括截尾數(shù)據(jù)）逐個(gè)由小到大依次排列，并對(duì)其中的每個(gè)死亡點(diǎn)進(jìn)行死亡概率、生存概率和生存率進(jìn)行估計(jì)。3.2生存率估計(jì)與組間比較3.2.2生存率的組間比較在醫(yī)學(xué)隨訪研究中，通常將病人按隨機(jī)化方法分配到兩種或多種治療組中，然后隨訪觀察和比較其生存時(shí)間的長(zhǎng)短和生存率的大小，以此來(lái)考察各種治療方案的優(yōu)劣；或者分析和比較同一治療方案下具有不同特征病人的生存率的大小，以此來(lái)探討影響這種療法的因素。因此，生存率組間比較實(shí)際上是兩條或多條生存曲線的比較。生存率的假設(shè)檢驗(yàn)方法有參數(shù)法和非參數(shù)法兩類(lèi)。參數(shù)法要求生存時(shí)間已知服從于某種概率分布，對(duì)實(shí)際資料擬合分布并求得其相應(yīng)的參數(shù)，然后通過(guò)比較不同組的分布參數(shù)來(lái)比較生存率是否相同。非參數(shù)法對(duì)資料的分布沒(méi)有要求，適用面比較廣。常見(jiàn)的有Log-rank檢驗(yàn)、Wilcoxon檢驗(yàn)（Gehan檢驗(yàn)）和似然比檢驗(yàn)，似然比檢驗(yàn)要求資料服從指數(shù)分布才有效。這里主要介紹Log-rank檢驗(yàn)和Wilcoxon檢驗(yàn)兩種方法。1.對(duì)數(shù)秩和檢驗(yàn)log-rank檢驗(yàn)Log-rank檢驗(yàn)是Mantel等人在1966年提出的，這種方法是在組間生存率相同的檢驗(yàn)假設(shè)（H0）下，對(duì)每組生存數(shù)據(jù)依據(jù)在各個(gè)時(shí)刻尚存活的患者數(shù)和實(shí)際死亡數(shù)計(jì)算期望死亡數(shù)，然后將期望死亡數(shù)與實(shí)際死亡數(shù)進(jìn)行比較，作假設(shè)檢驗(yàn)。這種方法可適合兩組或多組生存率比較。這種方法在兩組生存率比較時(shí)，計(jì)算比較簡(jiǎn)單。Log-rank檢驗(yàn)的漸進(jìn)平均值E和方差V：分別為實(shí)驗(yàn)組和對(duì)照組在第K時(shí)間間隔第i個(gè)病人死亡前生存的人數(shù)，分別為實(shí)驗(yàn)組和對(duì)照組在第K時(shí)間間隔第i個(gè)病人死亡前的死亡危險(xiǎn)率，為第k個(gè)時(shí)間間隔的死亡人數(shù)。設(shè)，則，其中，，，將log-rank統(tǒng)計(jì)量的分布視作N（E，1），有：根據(jù)區(qū)間上兩組概率分配向量中的治愈率，很容易求出所需總樣本含量：式中，為試驗(yàn)組和對(duì)照組的事件發(fā)生率。在隨訪研究中，樣本含量除受統(tǒng)計(jì)學(xué)要求及治療效果影響外，還有許多不確定性影響因素，例如患者入組、失訪、治愈時(shí)間的分布，患者在試驗(yàn)階段的依從性，以及是否滿足比例風(fēng)險(xiǎn)等等。Log-rank檢驗(yàn)除考慮最后結(jié)局，還考慮了出現(xiàn)結(jié)局的時(shí)間，并充分利用失訪資料所提供的不完全信息。對(duì)于具體的試驗(yàn)，本法都能擬合一個(gè)獨(dú)特的生存過(guò)程，較好反應(yīng)實(shí)際情況，應(yīng)用靈活，因此是一種有效、可行的樣本含量估計(jì)方法，能更好適應(yīng)臨床試驗(yàn)的復(fù)雜性和多樣性，巧妙解決多種復(fù)雜因素并存對(duì)樣本含量的影響問(wèn)題。2.Wilcoxon檢驗(yàn)法當(dāng)g=1,2,…,g,…,m時(shí),Wilcoxon檢驗(yàn)法χ2統(tǒng)計(jì)量計(jì)算公式仍可表示為:χ2=sˊv-1sυ=m-1其中sˊ=(s1,s2,…,sg,…,sm-1),sˊ為向量s的轉(zhuǎn)置。Sg的計(jì)算公式為：sg= (11.21)V為（m-1）(m-1)矩陣，記為V=[Vgh](m-1)(m-1)Vgh的計(jì)算公式為：Vgh= (11.22)上面sg和vgh計(jì)算公式中wi為權(quán)重,這里wi=nI。3.兩種檢驗(yàn)方法的比較Logrank檢驗(yàn)法和Wilcoxon檢驗(yàn)法實(shí)際上可以用統(tǒng)一公式來(lái)表示，即Wilcoxon檢驗(yàn)法的公式，公式中的權(quán)重wi=1時(shí)為L(zhǎng)ogrank檢驗(yàn)法，wi=ni時(shí)為Wilcoxon檢驗(yàn)。因而可以發(fā)現(xiàn)，Logrank檢驗(yàn)對(duì)生存時(shí)間較長(zhǎng)的個(gè)體在檢驗(yàn)中權(quán)重較大，對(duì)生存時(shí)間較短的個(gè)體在檢驗(yàn)中權(quán)重較小，在生存率（曲線）比較中，這種方法對(duì)尾部較為敏感；而Wilcoxon檢驗(yàn)則與Logrank檢驗(yàn)相反，對(duì)生存時(shí)間較短的個(gè)體在檢驗(yàn)中權(quán)重較大，比較中對(duì)數(shù)據(jù)的頭部的差別較為敏感。從理論和實(shí)踐中均發(fā)現(xiàn)，當(dāng)生存資料的各死亡點(diǎn)的危險(xiǎn)率在兩組或多組間成比例時(shí)，Logrank檢驗(yàn)的效率高于Wilcoxon檢驗(yàn)法，宜選用Logrank檢驗(yàn)。當(dāng)生存資料各時(shí)點(diǎn)的危險(xiǎn)率服從其他狀態(tài)時(shí)，Wilcoxon檢驗(yàn)法效率高于Logrank檢驗(yàn)，宜選用Wilcoxon檢驗(yàn)法。4.Cox模型像通常的回歸分析一樣，人們也希望能建立起生存時(shí)間（因變量或反映變量）隨危險(xiǎn)因素（自變量或協(xié)變量）變化的回歸過(guò)程，以便對(duì)危險(xiǎn)因素的作用大小有一個(gè)全面的了解和掌握，并根據(jù)危險(xiǎn)因素的不同取值對(duì)生存概率進(jìn)行預(yù)測(cè)。由于很難獲得準(zhǔn)確的生存時(shí)間，前述目的較難直接實(shí)現(xiàn)。1972年Cox提出了比例危險(xiǎn)模型，簡(jiǎn)稱Cox模型。由于此模型在表達(dá)形式上與參數(shù)模型相似，但在對(duì)模型中的各參數(shù)進(jìn)行估計(jì)時(shí)卻不依賴于特定的假設(shè)，所以又稱為半?yún)?shù)模型。設(shè)是影響生存時(shí)間t的k個(gè)危險(xiǎn)因素。設(shè)為i名受試者在時(shí)刻t的風(fēng)險(xiǎn)率，即t時(shí)刻外后一瞬間的死亡速率。又設(shè)表示不受危險(xiǎn)因素x的影響下，在時(shí)刻t的風(fēng)險(xiǎn)率，又稱為基準(zhǔn)風(fēng)險(xiǎn)率或基準(zhǔn)函數(shù)。其模型的具體形式如下：式中，為第i名受試者生存到t時(shí)刻的危險(xiǎn)率函數(shù)，是當(dāng)所有的危險(xiǎn)因素（即）不存在時(shí)的基礎(chǔ)危險(xiǎn)率函數(shù)，是可能與生存時(shí)間有關(guān)的m個(gè)危險(xiǎn)因素所構(gòu)成的向量。4SAS程序及結(jié)果輸出1.Kaplan-Meier和Nelson-Aalen程序：dataa;inputt@@;ift<0thencensor=1;elsecensor=0;/*如果時(shí)間小于0，為刪失變量賦值為1，否則賦值為0*/if_n_<20thengroup='A';elsegroup='B';/*前二十五個(gè)數(shù)據(jù)為A組，剩下的為B組*/t=abs(t);/*t的標(biāo)準(zhǔn)化*/cards;459161213102328293132474141576274100139-20-258-2698101012142048707599105162169195220-161-199-217-245;proclifetestmethod=plnelsonplots=(s,ls,lls);/*利用lifetest過(guò)程進(jìn)行生存分析并作生存函數(shù)圖，pl為Kaplan-Meier，nelson為Nelson-Aalen*/timet*censor(1);/*制定時(shí)間變量和刪失變量，指出刪失變量時(shí)刪失變量的取值*/stratagroup;/*指定分組變量*/run;結(jié)果第一列至第八列分別是生存時(shí)間、累積生存率、死亡概率、累計(jì)生存率標(biāo)準(zhǔn)誤、累積危險(xiǎn)率、累積危險(xiǎn)標(biāo)準(zhǔn)誤差、已觀測(cè)到的失效時(shí)間的例數(shù)、尚未觀測(cè)到的失效或截尾例數(shù)。有*號(hào)者表示截尾觀測(cè)值。給出A組生存時(shí)間四分位數(shù)、點(diǎn)估計(jì)及95%可信區(qū)間，生存時(shí)間均數(shù)及其標(biāo)準(zhǔn)誤。結(jié)果表示，A組患者的平均生存期為33.368個(gè)月。給出B組生存時(shí)間四分位數(shù)、點(diǎn)估計(jì)及95%可信區(qū)間，生存時(shí)間均數(shù)及其標(biāo)準(zhǔn)誤。結(jié)果表示，平均生存期為122.731個(gè)月。為兩者患者總?cè)藬?shù)、死亡數(shù)、截尾數(shù)和截尾百分比。層間等效檢驗(yàn)中，p值均小于0.0001，拒絕原假設(shè)（），認(rèn)為方差不齊性。以下是Kaplan-Meier法作出的生存函數(shù)圖：為兩組患者的生存分布函數(shù)曲線。兩條曲線在開(kāi)始時(shí)重疊，A組在100個(gè)月終止，B組在260個(gè)月終止。A組生存率下降顯著高于B組，說(shuō)明B組比A組生存時(shí)間長(zhǎng)。是-LogS(t)對(duì)生存時(shí)間T的散點(diǎn)圖，呈非直線趨勢(shì)，說(shuō)明生存時(shí)間不呈指數(shù)分布。是Log（-LogS(t)）對(duì)LogT的散點(diǎn)圖，兩條線分別近似直線，說(shuō)明生存時(shí)間近似呈Weibull分布。是各組生存函數(shù)曲線齊性檢驗(yàn)。依次給出秩次統(tǒng)計(jì)量、Logrank統(tǒng)計(jì)量的協(xié)方差矩陣、Wilcoxon統(tǒng)計(jì)量的協(xié)方差矩陣、各組生存函數(shù)一致性檢驗(yàn)結(jié)果等?？偨Y(jié)論：結(jié)果表明此資料不服從指數(shù)分布，近似服從Weibull分布，故宜選用logrank法的結(jié)果，兩條生存曲線分布有顯著性差異（p<0.05）。B組患者的生存時(shí)間顯著長(zhǎng)于A組患者。2.cox模型dataa;/*定義數(shù)據(jù)集*/inputgroupmonth;censor=(month<0);/*時(shí)間小于0的值賦給刪失變量*/months=abs(month);/*把變量標(biāo)準(zhǔn)化*/cards;141519116112113110123128129131132147141141157162174110011391-201-2581-26928210210212114220248270275299210521622169219522202-1612-1992-2172-245;procphregdata=a;/*利用PHREG過(guò)程做Cox模型*/modelmonths*censor(1)=group;/*建立生存時(shí)間months，截尾指示變量為censor的COX模型，censor取0時(shí)表示刪失*/run;結(jié)果SAS系統(tǒng)PHREG過(guò)程模型信息數(shù)據(jù)集WORK.A因變量months刪失變量censor刪失值1結(jié)值處理BRESLOW讀取的觀測(cè)數(shù)使用的觀測(cè)數(shù)4242事件和刪失值個(gè)數(shù)匯總合計(jì)事件刪失刪失

百分比4235716.67收斂狀態(tài)滿足收斂準(zhǔn)則(GCONV=1E-8)。模型擬合統(tǒng)計(jì)量準(zhǔn)則沒(méi)有

協(xié)變量具有

協(xié)變量-2LOGL214.126210.155AIC214.126212.155SBC214.126213.711為輸入數(shù)據(jù)信息。從這里可以看出截尾事件和終點(diǎn)時(shí)間。本例截尾數(shù)據(jù)14例。檢驗(yàn)全局零假設(shè):BETA=0檢驗(yàn)卡方自由度Pr

>

卡方似然比3.971110.0463評(píng)分4.000710.0455Wald3.869310.0492最大似然估計(jì)值分析參數(shù)自由度參數(shù)

估計(jì)值標(biāo)準(zhǔn)

誤差卡方Pr

>

卡方危險(xiǎn)

比group1-0.701100.356423.86930.04920.496為模型的檢驗(yàn)結(jié)果，從結(jié)果可以看出，模型較好的擬合了研究數(shù)據(jù)，有統(tǒng)計(jì)學(xué)意義（p<0.05）。為參數(shù)檢驗(yàn)結(jié)果以及幾個(gè)描述統(tǒng)計(jì)量。輸出結(jié)果包括參數(shù)估計(jì)值、標(biāo)準(zhǔn)誤、Wald值、p值、相對(duì)危險(xiǎn)比。變量group在統(tǒng)計(jì)學(xué)上有顯著性差異（p=0.0492）。HR=0.496。根據(jù)參數(shù)估計(jì)值，可寫(xiě)出一下Cox回歸方程：基于單片機(jī)和DSP的卷繞控制器數(shù)據(jù)采集和通訊設(shè)計(jì)基于MSP430單片機(jī)的柴油發(fā)電機(jī)監(jiān)控器的設(shè)計(jì)基于CPLD/FPGA和單片機(jī)的爆速儀設(shè)計(jì)基于單片機(jī)控制的晶閘管中頻感應(yīng)電源的研制基于十六位單片機(jī)的電力設(shè)備故障在線監(jiān)測(cè)裝置的設(shè)計(jì)與算法研究基于SPCE061A單片機(jī)的語(yǔ)音識(shí)別系統(tǒng)的研究基于PIC單片機(jī)的生物機(jī)能實(shí)驗(yàn)裝置的研究基于MotorolaMC68HC08系列單片機(jī)演示系統(tǒng)的設(shè)計(jì)與實(shí)現(xiàn)基于TCP/IP協(xié)議的單片機(jī)與INTERNET互連的設(shè)計(jì)與實(shí)現(xiàn)基于嵌入式實(shí)時(shí)操作系統(tǒng)和TCP/IP協(xié)議的單片機(jī)測(cè)控系統(tǒng)AVR8位嵌入式單片機(jī)在車(chē)載全球定位系統(tǒng)顯示終端中的應(yīng)用基于AVR單片機(jī)的250WHID燈電子鎮(zhèn)流器的研究基于單片機(jī)的TCP/IP技術(shù)研究及應(yīng)用基于P87C591單片機(jī)的CAN總線應(yīng)用層協(xié)議的研究基于單片機(jī)實(shí)現(xiàn)對(duì)二級(jí)倒立擺的控制C8051FXXX系列單片機(jī)仿真器的研制基于80C196MC單片機(jī)控制的變頻調(diào)速及配料控制系統(tǒng)的應(yīng)用研究基于單片機(jī)的膠印機(jī)控制系統(tǒng)開(kāi)發(fā)研究基于凌陽(yáng)單片機(jī)的二次壓降全自動(dòng)測(cè)量?jī)x的研制基于單片機(jī)的超聲測(cè)距系統(tǒng)基于MOTOROLA單片機(jī)的專(zhuān)用電池組智能充電儀全站儀動(dòng)態(tài)測(cè)量的研究以及其與單片機(jī)在軌道式龍門(mén)吊實(shí)時(shí)檢測(cè)中的應(yīng)用一種基于80C196KC單片機(jī)的新型電子負(fù)載的設(shè)計(jì)基于單片機(jī)的對(duì)講系統(tǒng)的研究開(kāi)發(fā)基于單片機(jī)的微波加熱瀝青路面再生修復(fù)機(jī)溫度控制器的開(kāi)發(fā)與研究基于單片機(jī)ATmega128的嵌入式工業(yè)控制器設(shè)計(jì)基于單片機(jī)的壓電閉環(huán)微位移控制系統(tǒng)的研究基于單片機(jī)的高壓靜電除塵整流設(shè)備的自動(dòng)監(jiān)控系統(tǒng)設(shè)計(jì)采用W78E58單片機(jī)的酸堿濃度檢測(cè)技術(shù)基于單片機(jī)的糧庫(kù)溫度監(jiān)控系統(tǒng)設(shè)計(jì)基于單片機(jī)控制的微型軸流式血泵外磁驅(qū)動(dòng)系統(tǒng)研究基于AVR單片機(jī)的電動(dòng)自行車(chē)控制系統(tǒng)研究基于PIC單片機(jī)的配電網(wǎng)綜合參數(shù)測(cè)控儀研究全自動(dòng)包裝機(jī)的單片機(jī)控制研究基于單片機(jī)系統(tǒng)的Java虛