針對不同群體設計不同特性的相機產品-資料下載頁
2025-06-26 04:11本頁面
【正文】 1:4)for(j in 1:5)[i,j]=sum(d3==i amp。 c16==j)par(mfrow=c(2,2))for(i in 1:4)barplot([i,],main=rownames()[i])=matrix(rep(0,10),nrow=2,byrow=T,dimnames=list(c(男,女),c(=1000,10012000,20013000,30016000,6000)))for(i in 1:2)for(j in 1:5)[i,j]=sum(d1==i amp。 a03==j)for(i in 1:2)barplot([i,],main=rownames()[i])par(mfrow=c(1,1))boxplot(a03[d1==1],a03[d1==2])=matrix(rep(0,20),nrow=4,dimnames=list(c(學生,企業(yè)大堂經(jīng)理,一般行政,專業(yè)技術),c(1000,10002000,20003000,30006000,6000)))for(i in 1:4) 職業(yè)、價格列聯(lián)表for(j in 1:5){[i,j]=sum(d3==i amp。 a03==j)}par(mfrow=c(2,2))for(i in 1:4)barplot([i,],main=rownames()[i])boxplot(a03[d3==1],a03[d3==2],a03[d3==3],a03[d3==4])Ridit檢驗部分：Ridit檢驗函數(shù)，主要是算出W值，求出p值，并且輸出置信區(qū)間，看重復程度，檢驗函數(shù)如下：=function(x){O=sum(x) 先求出幾個O.., r和sr=nrow(x)s=ncol(x)H=*apply(x,2,sum) 求H、RjC=rep(0,s)for(i in 2:s){C[i]=2*sum(H[1:i1])}N=H+C=N/O=apply(x,1,sum) 求行和=apply(x,2,sum) 求列和計算Ri R=matrix(rep(0,r*s),nrow=r,ncol=s)for(i in 1:r)for(j in 1:s){R[i,j]=x[i,j]/[j]*[j]}=apply(R,1,sum)計算Ridit得分矩陣Rid=matrix(rep(0,r*s),nrow=r)for(i in 1:r){Rid[i,]=*x[i,]}計算各處理的秩平均Ri_bar=rep(0,s)for(i in 1:r){Ri_bar[i]=sum(Rid[i,])/[i]} 求平均Ridit的置信區(qū)間=matrix(rep(0,2*r),nrow=r)for(i in 1:r){[i,]=c(Ri_bar[i]1/sqrt(3*[i]),Ri_bar[i]+1/sqrt(3*[i]))}dimnames()=list(NULL,c(置信下限,置信上限))作圖plot(rep(1,100),seq([1,1],[1,2],length=100),type=l,main=Ridit value confidence inetrval,xlab=group,ylab=R score,xlim=c(1,r),ylim=c(0,1))for(i in 2:r){lines(rep(i,100),seq([i,1],[i,2],length=100))}lines(c(0,r+1),c(,))W=12*sum(*()^2) 計算檢驗統(tǒng)計量pvalue=1pchisq(W,r1)list(Ridit Score=Rid,秩平均=Ri_bar,置信區(qū)間=,pvalue=pvalue)}找出需要檢驗的部分，以矩陣形式輸出：camera=(e:\\學習\\課件\\大三上\\非參數(shù)統(tǒng)計\\各章數(shù)據(jù) \\,header=T) malec01=femalec01=matrix(rep(0,20),4,5)for(i in 1:4){for(j in 1:5) {malec01[i,j]=sum(camera$d1==1amp。camera$d3==iamp。camera$c01==j)}}for(i in 1:4){for(j in 1:5) {femalec01[i,j]=sum(camera$d1==2amp。camera$d3==iamp。camera$c01==j)}}以c01和c02為例做Ridit檢驗mc01=malec01+femalec01(malec01)(femalec01)(mc01)關聯(lián)規(guī)則部分：camera=(39。C:\\Documents and Settings\\Administrator\\桌面\\39。,head=T)cameranew1=camera[,c(1:5,18)]ln1=length(cameranew1[1,])for(i in 1:(ln12)){cameranew1[[i]][cameranew1[[i]]=3]=0。cameranew1[[i]][cameranew1[[i]]3]=1。}cameranew1[[ln11]][cameranew1[ln11]==1]=0。cameranew1[[ln11]][cameranew1[ln11]==2]=1。cameranew1[[ln1]][cameranew1[ln1]==1]=1。cameranew1[[ln1]][cameranew1[ln1]=2]=0。cameranew1cameranew1=(cameranew1)。 trans2 as(cameranew1, transactions)。 summary(trans2)summary( trans2)itemFrequencyPlot(trans2, support = , = ) rules = apriori(trans2, parameter = list(support = ,confidence = ,maxlen=2))求規(guī)則inspect(rules,subset=lift)names(rules)x=subset(rules, subset = lhs %in% d3 )length(x)inspect(x)inspect(sort(x, by = support)[1:length(x)]) inspect(sort(x, by = confidence)[1:length(x)]) inspect(sort(x, by = lift)[1:length(x)])35 / 35
醫(yī)療健康相關推薦
鄂ICP備17016276號-1