2017-10-20 1 views
0

처음으로 글을 작성 했으므로 누락 된 부분이 있으면 사과드립니다.카이 제곱 검정의 예상 데이터를 제공합니다.

저는 두 가지 DNA 정량법을 비교 중이며 어느 방법이 실제 결과 (내가 샘플에 추가 한 DNA)에 가까운 결과를 얻을 수 있는지 알아 보려고합니다.

각 방법마다 5 복제본이 있습니다. 1000 900 - 500, 400, 400, 500, 500

  • 방법 Q -

    • 방법 P : I는 다음의 정량 값 (관측 값)가, 각 샘플 (기대치)를 1000 세포에 넣어있어 1400, 700, 1000

    chisq() 함수를 사용하려고하면 함수에 예상 값이 무엇인지 알려주지 않고 예상 값을 계산하며 원하는 값이 아닌 것입니다. .

    > P<-c(500, 400, 400, 500, 500) 
    
    > Q<-c(1000, 900, 1400, 700, 1000) 
    
    > chisqQ <- chisq.test(Q) 
    Chi-squared test for given probabilities 
    
    data: Q 
    X-squared = 260, df = 4, p-value < 2.2e-16 
    
    > chisqP <- chisq.test(P) 
    
        Chi-squared test for given probabilities 
    
    data: P 
    X-squared = 26.087, df = 4, p-value = 3.039e-05 
    

    이 문제는 내 기대 값을 설정하지 않는 것이, 그리고 Q를 위해 자동으로 1000을 계산하면서 P를 위해이 페이지 인수

    > round(chisqQ$expected,2) 
    [1] 1000 1000 1000 1000 1000 
    
    > round(chisqP$expected,2) 
    [1] 460 460 460 460 460 
    

    이되지 않습니다 chisq 함수에서,하지만 그것은 내 경우가 아닌 확률이어야합니다.

    나는 chi square 값을 손으로 계산하여 비교하지만, 미래에 한 번 R에서 할 수 있는지 알고 싶은 몇 가지 기술과 세포 양을 갖게 될 것입니다.

    나는이 상황을 제대로 이해한다면, 당신은 이미 예상 값을 가지고 당신이 적합 얼마나 좋은 볼 카이 제곱을 사용하려면 사전에

    감사합니다,

    건배,

    조안나 그래서

  • +0

    재현 가능한 데이터 또는 예제를 제공하고 지금까지 작성한 코드를 표시하십시오. 그러면 답변을 얻을 수 있습니다 ... – Thai

    답변

    0

    .

    그렇다면 다음 해결책이 효과적입니다.

    obs <- c(500,400,400,500,500) 
    exp <- c(XX, XX, XX, XX, XX) 
    chisq.test(x = observed, p = expected) 
    
    관련 문제