2016-03-07 2 views
1

Спарк 1,6Спарк 1,6 корреляции Пирсона

Если у меня есть набор данных, и я хочу, чтобы identifiy особенности в наборе данных с наибольшей предсказательной силой с помощью корреляции Пирсона, какие инструменты я должен использовать?

Наивный подход, который я использовал ... был:

val columns = x.columns.toList.filterNot(List("id","maxcykle","rul") contains) 
val corrVithRul = columns.map(c => (c,x.stat.corr("rul", c, "pearson"))) 

Output: 

    columns: List[String] = List(cykle, setting1, setting2, setting3, s1, s2, s3, s4, s5, s6, s7, s8, s9, s10, s11, s12, s13, s14, s15, s16, s17, s18, s19, s20, s21, label1, label2, a1, sd1, a2, sd2, a3, sd3, a4, sd4, a5, sd5, a6, sd6, a7, sd7, a8, sd8, a9, sd9, a10, sd10, a11, sd11, a12, sd12, a13, sd13, a14, sd14, a15, sd15, a16, sd16, a17, sd17, a18, sd18, a19, sd19, a20, sd20, a21, sd21) 
    corrVithRul: List[(String, Double)] = List((cykle,-0.7362405993070199), (setting1,-0.0031984575547410617), (setting2,-0.001947628351500473), (setting3,NaN), (s1,-0.011460304217886725), (s2,-0.6064839743782909), (s3,-0.5845203909175897), (s4,-0.6789482333860454), (s5,-0.011121400898477964), (s6,-0.1283484484732187), (s7,0.6572226620548292), (s8,-0.5639684065744165), (s9,-0.3901015749180319), (s10,-0.04924720421765515), (s11,-0.6962281014554186), (s12,0.6719831036132922), (s13,-0.5625688251505582), (s14,-0.30676887025759053), (s15,-0.6426670441973734), (s16,-0.09716223410021836), (s17,-0.6061535537829589), (s18,NaN), (s19,NaN), (s20,0.6294284994377392), (s21,0.6356620421802835), (label1,-0.5665958821050425), (label2,-0.548191636440298), (a1,0.040592887198906136), (sd1,NaN), (a2,-0.7364292... 

который, конечно, представив одну работу за карте итерации, Statistics.corr может быть то, что я ищу?

ответ

4

Statistics.corr выглядит как правильный выбор здесь. Другими параметрами, которые вы можете рассмотреть, являются RowMatrix.columnSimilarities (сходство с косинусом между столбцами, опционально с оптимизированной версией, которая использует выборку с порогом) и RowMatrix.computeCovariance. Так или иначе вам придется сначала собрать ваши данные в Vectors. Если предположить, что столбцы уже из DoubleType вы можете использовать VectorAssembler:

import org.apache.spark.ml.feature.VectorAssembler 
import org.apache.spark.mllib.linalg.Vector 

val df: DataFrame = ??? 

val assembler = new VectorAssembler() 
    .setInputCols(df.columns.diff(Seq("id","maxcykle","rul"))) 
    .setOutputCol("features") 

val rows = assembler.transform(df) 
    .select($"features") 
    .rdd 
    .map(_.getAs[Vector]("features")) 

Далее вы можете использовать Statistics.corr

import org.apache.spark.mllib.stat.Statistics 

Statistics.corr(rows) 

или преобразовать в RowMatrix:

import org.apache.spark.mllib.linalg.distributed.RowMatrix 

val mat = new RowMatrix(rows) 

mat.columnSimilarities(0.75) 
+0

Я не знаю, как interpit матрицу из S.corr, если столбец «x» - мой столбец предсказателя/цели – oluies

+0

К сожалению, вам придется идти по индексу, поэтому полезно поместить depende nt переменная. В противном случае это просто корреляционная матрица, как вы можете ожидать. – zero323

+0

должен, вероятно, нормализовать и до корреляции? – oluies

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