楼主: perfect-秋
7536 15

[问答] 关于核主成分分析 [推广有奖]

  • 0关注
  • 0粉丝

大专生

95%

还不是VIP/贵宾

-

威望
0
论坛币
29 个
通用积分
0.0095
学术水平
0 点
热心指数
3 点
信用等级
0 点
经验
1610 点
帖子
92
精华
0
在线时间
20 小时
注册时间
2010-10-18
最后登录
2012-5-3

楼主
perfect-秋 发表于 2011-9-14 19:37:02 |AI写论文
5论坛币
我这边有一组数据,想做logistic回归,里面的自变量基本上都是分类变量,自变量有20多个,我想先提取主成分,但是PCA比较适合于线性数据,这种二分类的不适合,看过文献提示可用KPCA,但是没用过MATLAB编写这个程序,想请各位大侠帮帮忙哈,这个MATLAB的程序是咋样哈???如有需要我可以把数据提上来哈!

关键词:主成分分析 核主成分 主成分 logistic回归 logistic 分析 主成分分析法 spss主成分分析 逐步回归分析 多元回归分析 因子分析法 应用时间序列分析

本帖被以下文库推荐

沙发
perfect-秋 发表于 2011-9-15 18:39:25
没人顶啊,各位大侠帮帮忙啊!

藤椅
epoh 发表于 2011-9-16 20:49:10
Kernel Statistics toolbox
Table of Contents
=================
- Introduction
- Key Features
- Data Format
For classification
For regression
- Code Usage with Examples
KDR
UseKDR
KPCA
KSIR
- Dimension Reduction Using KPCA or KSIR with Examples
KPCA procedure
KSIR procedure for classification
KSIR procedure for regression

KStattoolbox.zip
   KStattoolbox.zip (224.4 KB) 本附件包括:
  • Housing_dataset.txt
  • Ionosphere_dataset.mat
  • KDR.m
  • KPCA.m
  • KSIR.m
  • README
  • README.pdf
  • srsplit.m
  • SVKernel_C.m
  • SVKernel_EX.dll
  • UseKDR.m

已有 2 人评分学术水平 热心指数 信用等级 收起 理由
zhangtao + 5 + 5 + 5 感谢epoh大师!向老师问好!
ywh19860616 + 1 + 1 + 1 热心助人

总评分: 学术水平 + 6  热心指数 + 6  信用等级 + 6   查看全部评分

板凳
perfect-秋 发表于 2011-10-20 17:41:50
epoh 发表于 2011-9-16 20:49
Kernel Statistics toolbox
Table of Contents
=================
不好意思,现在才看到哈
那请问大师,这样的数字可以用核主成分分析+logistic分析它的影响因素嘛?
还有就是,这个你给的程序我直接用数据套用嘛?

报纸
epoh 发表于 2011-10-20 22:04:17

参考范例

数据格式,直接参照Ionosphere_dataset.mat

将数据存成.txt也可以

他的label是1,-1

你要改成1,2也可以

readme.pdf page 6/7
Dimension Reduction Using KPCA or KSIR with Examples
KPCA procedure
%%%%%%%%
%load data
load Ionosphere_dataset.mat
%Construct PCs via KPCA
[Info] = KDR([], inst, '-s 0 -t 2 -g 0.1 -p 10');

%Get the projected inst form the PCs
[ProjInst] = UseKDR(inst, Info)

地板
perfect-秋 发表于 2011-10-21 16:50:07
弱弱的说句,跑不出来

7
epoh 发表于 2011-10-21 17:53:33
%1.get current working directory.
pwd
  %C:\Users\...\Documents\MATLAB
%2.put KStattoolbox folder in current working directory
%3.set path in matlab:
      File\Set Path\add folder:
      C:\Users\...\Documents\MATLAB\KStattoolbox
%4.in matlab command window
load Ionosphere_dataset.mat
[Info] = KDR([], inst, '-s 0 -t 2 -g 0.1 -p 10');
[ProjInst] = UseKDR(inst, Info)
已有 1 人评分学术水平 热心指数 信用等级 收起 理由
zhangtao + 5 + 5 + 5 精彩帖子

总评分: 学术水平 + 5  热心指数 + 5  信用等级 + 5   查看全部评分

8
zhangtao 发表于 2011-10-23 09:40:37

ProjInst =

  Columns 1 through 7

   -0.9648   -0.5063   -0.3472   -0.1082    0.3419    0.3317    0.1697
   -1.1787   -0.5114   -0.2184    0.0824    0.2076    0.1918    0.0523
   -0.7004   -0.6757   -0.2881   -0.1185    0.3505    0.2676    0.2273
   -0.5537   -0.3573   -0.3622   -0.2014    0.4023    0.4238    0.2828
   -0.8313   -0.2382   -0.4000   -0.0958    0.4347    0.4126    0.1361
   -0.6519   -0.2053   -0.3992   -0.1424    0.4533    0.4663    0.2050
   -0.9755   -0.2114   -0.3909   -0.0020    0.4011    0.3279    0.0236
   -1.0575   -0.2045    0.1558    0.3684   -0.1596   -0.0777   -0.0716
   -1.2431   -0.3549    0.0449    0.2952   -0.0456    0.0130   -0.0475
   -1.0042   -0.7466    0.1517    0.1521   -0.1417    0.0339    0.0571
   -1.0911   -0.2142   -0.3409    0.0873    0.3171    0.2358   -0.0591
   -1.2483   -0.5202   -0.0536    0.1703    0.0686    0.1281    0.0393
   -1.1743   -0.3400   -0.2738    0.1117    0.2491    0.2014   -0.0193
   -0.5734   -0.2450   -0.3345   -0.1331    0.3804    0.3772    0.1985
   -0.4795   -0.1771    0.4167    0.2798   -0.2415    0.1127    0.2328
   -0.8991   -0.0973   -0.3694    0.0595    0.3746    0.2601   -0.0829
   -0.9958   -0.3632   -0.3946   -0.0734    0.3905    0.3533    0.1175
   -0.9483   -0.1806   -0.1277    0.1728    0.0777   -0.0045   -0.1287
   -1.2194   -0.2367   -0.2073    0.2096    0.1680    0.0883   -0.1296
   -1.2535   -0.4128    0.0578    0.2699   -0.0439    0.0285   -0.0150
   -0.6738   -0.5346    0.5875    0.2848   -0.1522    0.3004    0.2590
   -0.6957   -0.7465    0.4728    0.2312   -0.0930    0.2911    0.1456
   -0.8378   -0.2184   -0.4197   -0.0768    0.4470    0.4131    0.1145
   -1.0995   -0.2210   -0.3467    0.0837    0.3267    0.2369   -0.0568
   -1.2289   -0.2922   -0.2298    0.1683    0.1933    0.1274   -0.0852
   -1.0145   -0.3310   -0.0968    0.1339    0.0613    0.0108   -0.0489
   -0.0180   -0.2720   -0.3217    0.4414   -0.1079    0.0850    0.2182
   -1.1372   -0.2148   -0.3124    0.1337    0.2689    0.1746   -0.1090
   -1.0672   -0.2004   -0.2094    0.1874    0.1826    0.1053   -0.1118
   -1.1814   -0.2868    0.1015    0.3423   -0.1132   -0.0346   -0.0671
   -1.0716   -0.4047    0.3059    0.3405   -0.2016    0.0549    0.0972
   -0.7762   -0.4875    0.5298    0.3098   -0.2155    0.2212    0.2473
   -1.0615   -0.2976   -0.1161    0.1941    0.0441    0.0856   -0.0729
   -1.1641   -0.2291   -0.0962    0.2762    0.0527    0.0164   -0.1447
   -0.0672   -0.4669   -0.4333    0.5728   -0.0657    0.0455    0.2771
   -0.7910   -0.2271    0.3753    0.3598   -0.2709    0.0205    0.1314
   -0.6092   -0.4453    0.5868    0.2884   -0.1877    0.2881    0.2834
   -0.8105   -0.3677    0.4763    0.3347   -0.2498    0.1419    0.2250
   -0.8913   -0.4065    0.4500    0.3368   -0.2361    0.1289    0.2004
   -0.2610   -0.2157    0.4403    0.2002   -0.1063    0.2827    0.2373
   -0.7057   -0.0800   -0.3129    0.0437    0.3097    0.2321   -0.0511
   -0.6235   -0.5900    0.6040    0.2592   -0.1137    0.3310    0.2370
   -1.2235   -0.6144   -0.1149    0.0920    0.0667    0.0943    0.0366
   -1.2493   -0.5192   -0.1487    0.1428    0.0242    0.0544   -0.0506
   -1.1050   -0.8722   -0.0085    0.0258   -0.0965    0.0490    0.0286
   -0.9640   -0.4307    0.3446    0.3045   -0.2407    0.0654    0.1312
   -1.0350   -0.8178   -0.1307    0.0579    0.0022    0.0939    0.0564
   -1.0717   -0.8425   -0.1150   -0.0310   -0.0090    0.0497    0.0397
   -0.1486   -0.6703    0.3950    0.3379    0.5542    0.1979   -0.3230
   -0.0864   -0.6111    0.2525    0.4021    0.6058    0.0499   -0.3340
   -1.2364   -0.6110    0.0445    0.1927   -0.0812    0.0555    0.0206
   -1.1879   -0.1903   -0.2240    0.2161    0.1771    0.0764   -0.1654
   -1.1429   -0.2383   -0.2905    0.1177    0.2422    0.1508   -0.1060
   -0.9730   -0.1365   -0.2637    0.1302    0.2417    0.1212   -0.1302
   -1.2233   -0.2777   -0.2082    0.1795    0.1661    0.0836   -0.1138
   -0.5167   -0.1023   -0.3416   -0.1113    0.3912    0.3960    0.1333
   -0.1761   -0.7636    0.3149    0.3866    0.5728    0.1174   -0.2975
   -0.1088   -0.6730    0.1933    0.4550    0.6016    0.0093   -0.2762
   -0.8617   -1.0920   -0.1045   -0.1161   -0.0193    0.0380    0.0003
   -1.2351   -0.5506    0.0797    0.2307   -0.0709    0.0533    0.0294
   -1.1964   -0.2356   -0.1922    0.2106    0.1438    0.0461   -0.1489
   -0.8521   -0.3242   -0.0533    0.1307   -0.0078    0.0161   -0.0320
   -1.0119   -0.9293   -0.1603   -0.0634    0.0375    0.0588    0.0548
   -0.8305   -0.4840    0.3994    0.2570   -0.1777    0.1389    0.1788
   -1.2204   -0.2196   -0.1375    0.2655    0.0867    0.0143   -0.1632
   -1.2115   -0.2009   -0.1155    0.2846    0.0648   -0.0083   -0.1757
   -0.8958   -1.0922   -0.0062   -0.0460   -0.0050    0.0793    0.0231
   -0.3002   -0.9636    0.2753    0.2193    0.4129    0.0804   -0.2472
   -0.1105   -0.7035    0.1088    0.4789    0.5852   -0.0573   -0.2220
   -1.0624   -0.8582   -0.0098   -0.0292   -0.0244    0.0235    0.0574
   -1.0881   -0.5570   -0.2171    0.0228    0.1769    0.1466    0.0588
   -1.2292   -0.2346   -0.1047    0.2799    0.0472   -0.0089   -0.1614
   -0.5928   -1.1654   -0.1033   -0.1179    0.0231   -0.0157   -0.0861
   -0.5608   -1.2109   -0.1018    0.0260    0.0354    0.0222   -0.1186
   -1.2030   -0.2041   -0.0803    0.2970    0.0374   -0.0266   -0.1693
   -1.2205   -0.2354   -0.0583    0.3015    0.0074   -0.0292   -0.1506
   -1.2414   -0.6154    0.0079    0.1505   -0.0736    0.0297    0.0074
   -0.1681   -0.8179   -0.4624    0.6772    0.0216    0.0190    0.1769
   -1.2844   -0.4619   -0.0784    0.1947    0.0139    0.0224   -0.0514
   -0.3399   -1.0937   -0.3019    0.3665    0.0582   -0.0121   -0.0409
   -0.2288   -0.6186   -0.3076    0.2002   -0.4028    0.4911   -0.4022
   -1.2737   -0.3913   -0.1162    0.2081    0.0473    0.0385   -0.0851
   -0.1528   -0.5511   -0.3476    0.2655   -0.4255    0.5315   -0.4109
   -0.1048   -0.4039   -0.3172    0.2593   -0.4148    0.5493   -0.3994
   -0.0953   -0.3814   -0.2319    0.1515   -0.3541    0.4441   -0.3921
   -1.1029   -0.4418    0.1456    0.2750   -0.1438    0.0121    0.0169
   -1.2205   -0.3227   -0.1289    0.2221    0.0567    0.0235   -0.1079

   

  Columns 8 through 10

    0.0347    0.0770   -0.1116
    0.0174   -0.0038   -0.0519
   -0.0035    0.1400   -0.1641
    0.0731    0.2640   -0.1467
    0.1781    0.1707   -0.0233
    0.1831    0.2348   -0.0588
    0.2198    0.0953    0.0259
    0.2095    0.1003   -0.2035
    0.1107    0.0462   -0.0997
   -0.0471   -0.0005    0.0172
    0.1914    0.0854    0.0619
   -0.0122    0.0635   -0.0184
    0.1115    0.0588    0.0233
    0.1211    0.0968   -0.1440
    0.3052    0.2209   -0.3049
    0.3175    0.0341    0.0532
    0.1080    0.1114   -0.0233
    0.1391    0.1044    0.1177
    0.1793    0.0246    0.0065
    0.0740    0.0382   -0.0868
    0.2265    0.0704   -0.0956
    0.1409   -0.0158   -0.0428
    0.2091    0.1504   -0.0158
    0.2095    0.0388    0.0234
    0.1439    0.0160    0.0009
    0.1001   -0.1797   -0.1503
    0.0868   -0.2023    0.1980
    0.2094    0.0031    0.0129
    0.1888    0.0089   -0.0087
    0.1416    0.1281   -0.0978
    0.1206    0.1579   -0.1607
    0.1961    0.1302   -0.1549
    0.0316    0.1285    0.0339
    0.1687   -0.0365   -0.1019
    0.1183   -0.1482    0.1503
    0.2431    0.2627   -0.2558
    0.2657    0.1145   -0.1331
    0.2192    0.1868   -0.2215
    0.1932    0.1591   -0.2125
    0.2901    0.0751   -0.1061
    0.2388    0.1454    0.1511
    0.2303    0.0315   -0.0411
   -0.0393   -0.0322   -0.0422
   -0.0200    0.0147    0.0294
   -0.0583   -0.0207   -0.0267
    0.0969    0.2227   -0.0988
   -0.1253    0.0092    0.0008
   -0.0783   -0.0416   -0.0040
   -0.0101   -0.0443    0.0353
   -0.1196    0.0779   -0.1283
   -0.0308   -0.0126   -0.0790
    0.2181    0.0151    0.0171
    0.1942   -0.0164    0.0053
    0.2604   -0.0806   -0.0271
    0.1818   -0.0658   -0.0640
    0.2389    0.1574   -0.0541
   -0.0354    0.0218   -0.0480
   -0.1268    0.1112   -0.1735
    0.0253    0.0373    0.0020
    0.0139    0.0402   -0.0625
    0.1752   -0.0120    0.0246
   -0.0006    0.0534    0.0340
   -0.0422   -0.1076   -0.0952
    0.1148    0.0816   -0.1407
    0.1795    0.0276   -0.0059
    0.1960    0.0366   -0.0128
    0.0275   -0.0220   -0.0416
    0.1390    0.0219    0.0169
   -0.1169    0.1536   -0.2289
   -0.0580   -0.0709   -0.0077
    0.0075   -0.0783   -0.0969
    0.1704    0.0315   -0.0228
    0.2103    0.0889   -0.0329
    0.1923    0.1109   -0.0102
    0.2078    0.0148   -0.0519
    0.1554    0.0688   -0.0135
   -0.0483    0.0377    0.0022
    0.1305   -0.0085    0.0199
    0.0276   -0.0148   -0.0313

>>
数学好就是要天天学

9
zhangtao 发表于 2011-10-23 09:46:11
以上结果是什么意思,看不懂?
数学好就是要天天学

10
perfect-秋 发表于 2011-11-26 14:34:12
epoh 发表于 2011-10-21 17:53
%1.get current working directory.
pwd
  %C:\Users\...\Documents\MATLAB
哎,我知识太浅薄了,需要现实中请教老师哈,不过这个论坛币给你哈

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群
GMT+8, 2026-1-8 04:34