prdata.html

一个关于数据聚类和模式识别的程序,在生物化学,化学中因该都可以用到.希望对大家有用,谢谢支持

0046         DS.dataName=<span class="string">'abalone'</span>;
0047         inputName={<span class="string">'Sex'</span>, <span class="string">'Length'</span>, <span class="string">'Diameter'</span>, <span class="string">'Height'</span>, <span class="string">'Whole weight'</span>, <span class="string">'Shucked weight'</span>, <span class="string">'Viscera weight'</span>, <span class="string">'Shell weight'</span>};
0048         [DS.inputName]=deal(inputName);
0049         DS.input=abalone(:, 1:end-1)';
0050         DS.output=abalone(:, end)';
0051         <span class="keyword">if</span> nargout==2
0052             dataNum=size(DS.input, 2);
0053             TS=DS;
0054             DS.input= DS.input(:, 1:2:dataNum);
0055             DS.output=DS.output(:, 1:2:dataNum);
0056             TS.input= TS.input(:, 2:2:dataNum);
0057             TS.output=TS.output(:, 2:2:dataNum);
0058         <span class="keyword">end</span>
0059     <span class="keyword">case</span> <span class="string">'random2'</span>        <span class="comment">% 2D random data</span>
0060         n=100;
0061         dim=2;
0062         c1 = [0.125 0.25]'; data1 = randn(dim,n)/8 + c1*ones(1,n); out1 = 1*ones(1,n);
0063         c2 = [0.625 0.25]'; data2 = randn(dim,n)/8 + c2*ones(1,n); out2 = 2*ones(1,n);
0064         c3 = [0.375 0.75]'; data3 = randn(dim,n)/8 + c3*ones(1,n); out3 = 3*ones(1,n);
0065         c4 = [0.875 0.75]'; data4 = randn(dim,n)/8 + c4*ones(1,n); out4 = 4*ones(1,n);
0066         data = [data1, data2, data3, data4];
0067         out = [out1, out2, out3, out4];
0068         DS.input=data;
0069         DS.output=out;
0070         DS.dataName=<span class="string">'random2'</span>;
0071         <span class="keyword">if</span> nargout==2
0072             dataNum=size(DS.input, 2);
0073             TS=DS;
0074             DS.input= DS.input(:, 1:2:dataNum);
0075             DS.output=DS.output(:, 1:2:dataNum);
0076             TS.input= TS.input(:, 2:2:dataNum);
0077             TS.output=TS.output(:, 2:2:dataNum);
0078         <span class="keyword">end</span>
0079     <span class="keyword">case</span> <span class="string">'random3'</span>        <span class="comment">% 3D random data</span>
0080         dataNum=100;
0081         mean1=[0 0 0]';
0082         input1=randn(3, dataNum)+mean1*ones(1,dataNum);
0083         output1=1*ones(1, dataNum);
0084         mean2=[0 5 5]';
0085         input2=randn(3, dataNum)+mean2*ones(1,dataNum);
0086         output2=2*ones(1, dataNum);
0087         mean3=[3 2 4]';
0088         input3=randn(3, dataNum)+mean3*ones(1,dataNum);
0089         output3=3*ones(1, dataNum);
0090         DS.input=[input1, input2, input3];
0091         DS.output=[output1, output2, output3];
0092     <span class="keyword">case</span> <span class="string">'random6'</span>        <span class="comment">% 6D random data</span>
0093         n=50;
0094         dim=6;
0095         c1 = [1 0 0 0 0 0]'; data1 = randn(dim,n)/4 + c1*ones(1,n); out1 = 1*ones(1,n);
0096         c2 = [0 1 0 0 0 0]'; data2 = randn(dim,n)/4 + c2*ones(1,n); out2 = 2*ones(1,n);
0097         c3 = [0 0 1 0 0 0]'; data3 = randn(dim,n)/4 + c3*ones(1,n); out3 = 3*ones(1,n);
0098         c4 = [0 0 0 1 0 0]'; data4 = randn(dim,n)/4 + c4*ones(1,n); out4 = 4*ones(1,n);
0099         data = [data1, data2, data3, data4];
0100         out = [out1, out2, out3, out4];
0101         DS.input=data;
0102         DS.output=out;
0103         DS.dataName=<span class="string">'random6'</span>;
0104         <span class="keyword">if</span> nargout==2
0105             dataNum=size(DS.input, 1);
0106             TS=DS;
0107             DS.input= DS.input(:, 1:2:dataNum);
0108             DS.output=DS.output(:, 1:2:dataNum);
0109             TS.input= TS.input(:, 2:2:dataNum);
0110             TS.output=TS.output(:, 2:2:dataNum);
0111         <span class="keyword">end</span>
0112     <span class="keyword">otherwise</span>
0113         disp(<span class="string">'Unknown method.'</span>)
0114 <span class="keyword">end</span>
0115 
0116 <span class="keyword">if</span> ~isfield(DS, <span class="string">'inputName'</span>)
0117     <span class="keyword">for</span> i=1:size(DS.input,1)
0118         DS.inputName{i}=int2str(i);
0119         <span class="keyword">if</span> exist(<span class="string">'TS'</span>)
0120             TS.inputName{i}=int2str(i);
0121         <span class="keyword">end</span>
0122     <span class="keyword">end</span>
0123 <span class="keyword">end</span>
0124 
0125 <span class="comment">% ======= Self demo</span>
0126 <a name="_sub1" href="#_subfunctions" class="code">function selfdemo</a>
0127 subplot(2,2,1);
0128 DS=feval(mfilename, <span class="string">'random2'</span>); <a href="dcprDataPlot.html" class="code" title="function dcprDataPlot(DS, plotTitle, displayAnnotation)">dcprDataPlot</a>(DS);
0129 subplot(2,2,2);
0130 DS=feval(mfilename, <span class="string">'iris'</span>); <a href="dcprDataPlot.html" class="code" title="function dcprDataPlot(DS, plotTitle, displayAnnotation)">dcprDataPlot</a>(DS);
0131 subplot(2,2,3);
0132 DS=feval(mfilename, <span class="string">'wine'</span>); <a href="dcprDataPlot.html" class="code" title="function dcprDataPlot(DS, plotTitle, displayAnnotation)">dcprDataPlot</a>(DS);
0133 subplot(2,2,4);
0134 DS=feval(mfilename, <span class="string">'abalone'</span>); <a href="dcprDataPlot.html" class="code" title="function dcprDataPlot(DS, plotTitle, displayAnnotation)">dcprDataPlot</a>(DS);</pre></div>
<hr><address>Generated on Thu 30-Oct-2008 12:53:56 by <strong><a href="http://www.artefact.tk/software/matlab/m2html/">m2html</a></strong> &copy; 2003</address>
</body>
</html>