2013-12-04 3 views
0

У меня есть небольшое сомнение в R, относящееся к LDA. Как и в spss, когда я пытался получить коэффициенты классификации классификаторов для линейки с линейным дискриминантным анализом в R с пакетом MASS, я получать только коэффициенты линейного дискриминанта, как следующие:Коэффициенты классификации Fisher для нескольких классов в LDA в R

коэффициенты линейных дискриминантами:

     LD1  LD2    LD3   LD4 
    Var1  0.018952518 0.010766163  0.534278507  -2.32E-02 
    Var2  -0.000827315 -0.012934214 -0.013986988  -3.14E-01 
    Var3  -3.616088667 -0.18684861  -2.962979702  1.36E-01 
    Var4  0.000139365 -0.003802969  0.000313853  9.33E-05 
    Var5  0.007675119 0.006891405  0.05536683  -9.69E-02 


       LD5    
    Var1  5.64E-02   
    Var2  6.04E-02   
    Var3  -3.69E-01   
    Var4  -3.40E-05   
    Var5  -4.92E-01 

Но у меня есть 7 групп в моей группировке переменной, мне нужно, чтобы получить коэффициенты функции рыболовы классификации для каждого кластера всех пяти переменные, чтобы я мог использовать их для дальнейшего анализа в Excel. Я намерен получить похожую таблицу добавить ниже:

Классификация функций Коэффициенты

 Cluster    
       7   8   9  10  11   12  13 
    Var1  .630  .580  .555 .571  .598  .714 .642 
    Var2  .025  .028  .028 .029  .026  .029 .029 
    Var3  .685  .684  .752 .681  .678  .695 .700 
    Var4  1.019  .997  1.015 .998  1.023  1.033 1.033 
    Var5 17.331  21.253 21.457 21.347  9.166  8.850 8.860 
(Constant)-44.687 -57.762 -59.353-58.928 -36.337 -42.367 -42.744 

Так что я хочу такой же результат, как и выше в R, пожалуйста, помогите мне:

ответ

1

Если я правильно понимаю ваш вопрос, это может быть решение вашей проблемы:

Classification functions in linear discriminant analysis in R

сообщение содержит скрипт, который генерирует коэффициенты функции классификации из дискриминантных функций и добавляет их к результатам вашей функции lda() в виде отдельной таблицы.

+0

Я запустил код, который был найден в ur link на моем наборе данных, но я не смог увидеть результат. Мог бы направить меня, чтобы получить таблицу коэффициентов классификации с кодом ур. –

+0

Как вы его запустили? Вы вызывали функцию с требуемыми параметрами ('legendre.lda (x, groups)')? Вы уверены, что данные обработаны правильно? Вы получили сообщение об ошибке? Если таблица называется 'table1', результаты добавляются к исходным результатам как' table1 $ class.funs'. Если вы все сделали правильно, вы сможете увидеть их, если вы наберете 'print (table1 $ class.funs)'. – Fiete

+0

u знаю, что я получил таблицу, но дело в том, что здесь мне нужно получить соответствующие номера кластеров в виде имен col и соответствующих имен переменных в виде имен строк. и здесь, в таблице, я получаю NA для двух последних переменных каждой группы. И результат не должен содержать переменную группировки. я не знаю, почему, мог бы помочь мне. –

3

Вы можете изменить код, предложенный Fiete, или вы можете использовать linDA от DiscriMine package;

library(DiscriMiner) 
mylda = linDA(iris[,1:4], iris[,5]) 
summary(mylda) 
mylda$functions 

       setosa versicolor virginica 
constant  -86.30847 -72.852607 -104.36832 
Sepal.Length 23.54417 15.698209 12.44585 
Sepal.Width 23.58787 7.072510 3.68528 
Petal.Length -16.43064 5.211451 12.76654 
Petal.Width -17.39841 6.434229 21.07911 
+0

Спасибо Soooooooo Muchhhh Friend ........ Это то, что я хочу, наконец. –

Смежные вопросы