getScore: Get differential score for each observation

Description Usage Arguments Value Examples

View source: R/getScore.R

Description

Get differential score for each observation

Usage

1
getScore(object, type = "weight")

Arguments

object

Object of class inheriting from "mpk".

type

String indicating whether the score is associate to a difference in mixture weights type = "weight", or a difference in the mean parameters type = "shift".

Value

A vector with a score for each observation.

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
n = c(250, 250)
p = 4

Y1 = rbind( matrix( rnorm( n[1]*p), ncol = p), matrix( rnorm(n[2]*p) + 3, ncol = p))
Y2 = rbind( matrix( rnorm( n[1]*p), ncol = p), matrix( rnorm(n[2]*p) + 4, ncol = p))
Y = rbind(Y1, Y2)
C = c( rep(1,sum(n)), rep(2,sum(n)))

ans = mpk(Y, C)
score = getScore(ans, type = "shift")
hist(score)

jacsor/MPK documentation built on May 17, 2017, 6:49 a.m.