找到导致低p值的条目

时间:2015-10-09 09:52:44

标签: r statistics table-statistics

在R中有2个向量

public class MyIntentService extends IntentService {

    private static final String ACTION_FOO = "bubblebearapps.co.uk.scratchpad.action.FOO";

    // TODO: Rename parameters
    private static final String EXTRA_PAREM_TTL = "bubblebearapps.co.uk.scratchpad.extra.TTL";

    private static final long SERVICE_EXPIRE_TIME_MILLIS = 3000;

    /**
     * Starts this service to perform action Foo with the given parameters. If
     * the service is already performing a task this action will be queued.
     *
     * @see IntentService
     */
    // TODO: Customize helper method
    public static void startActionFoo(Context context) {
        Intent intent = new Intent(context, MyIntentService.class);
        intent.setAction(ACTION_FOO);
        intent.putExtra(EXTRA_PAREM_TTL, System.currentTimeMillis() + SERVICE_EXPIRE_TIME_MILLIS);
        context.startService(intent);
    }

    public MyIntentService() {
        super("MyIntentService");
    }

    @Override
    protected void onHandleIntent(Intent intent) {
        if (intent != null) {
            final String action = intent.getAction();
            if (ACTION_FOO.equals(action)) {
                final long param1 = intent.getLongExtra(EXTRA_PAREM_TTL, 0);
                handleActionFoo(param1);
            }
        }
    }

    private void handleActionFoo(long expiryTime) {

        if(System.currentTimeMillis() > expiryTime){
            //skip work and return
        }

    }
}

我想测试H:u <- c(109, 77, 57, 158, 60, 63, 42, 20, 139, 15, 64, 18) v <- c(734, 645, 1001, 1117, 1071, 687, 162, 84, 626, 64, 218, 79) u是独立的,所以我进行了卡方检验:

v

并获得非常低的p值意味着我可以拒绝H,这意味着chisq.test( as.data.frame( rbind(u,v) ) ) u不是独立的。

但是当我输入

v

我得到了0.23的p值,这意味着我可以接受H. 我应该选择这两个测试中的哪一个?

此外,我想在这些向量中找到导致这种低p值的条目。任何想法如何做到这一点?

1 个答案:

答案 0 :(得分:0)

测试统计使用标准化残差的平方和。您可以查看这些值以了解特定值的重要性

m = chisq.test(u, v)

residuals(m)
m$stdres