2014-12-24 7 views
1

를 벡터 변환을 "부팅", "객체의 이름을 수 있습니다 mymodel_boot ". 내가경험 t-값은, 내가 "bootsempls"클래스의 객체를 얻을 내가) (R의 semPLS 패키지의 기능을 bootsempls를 실행하면 콘솔 출력이

>mymodel_boot 

     Estimate  Bias Std.Error Lower Upper 
lam_1_1  0.923 -0.000188 0.01400 0.892 0.946 
lam_1_2  0.918 -0.000971 0.01778 0.876 0.945 
lam_2_1  0.860 -0.001325 0.02647 0.799 0.903 
lam_2_2  0.897 0.000579 0.01604 0.864 0.926 
lam_2_3  0.799 -0.000253 0.03476 0.723 0.859 
lam_3_1  0.793 -0.003462 0.06742 0.636 0.897 
lam_3_2  0.942 0.001877 0.00892 0.927 0.963 
beta_1_3 0.429 0.003359 0.07722 0.280 0.579 
beta_2_3 0.294 -0.000880 0.07339 0.147 0.435 

을 실행하면

나는 위의 표를 얻을.

는 지금 해당 표준 오차에 의해 원래의 계수를 나눌, t - 값을 계산한다.

원래 계수

> mymodel_boot$t0 
[1] 0.9227120 0.9180390 0.8599884 0.8973823 0.7993398 0.7927142 0.9424423 0.4292315 0.2941801 

이전 코드 블록에 인쇄 할 수처럼 Std.Error을 대상으로 어떤 방법이 있나요 대상이 될 수있다, 그래서 난 그냥이 나누어의 벡터를 얻을 수 t- 값?

bootsempls의 도움말 페이지를 읽었지만 Std.Err 만 반환하는 구성 요소를 찾을 수 없어서 혼란 스럽습니다. 그 이유는 t- 값을 추정하는 일반적인 절차를 생각하고 똑똑하지 않기 때문입니다. 콘솔에서 각 Std.Err을 복사하여 복사하는 것보다? 나는 그것을 믿을 수 없다. Std.Err은 $ 연산자를 통해 타겟으로 할 수없는 경우

, I는 벡터 내로 콘솔 출력 Std.Err에서 열 변형시킬 수있는 방법이 있는가?

답변

2

좋은 질문입니다. 일반적인 방법은 알 수없는 클래스 (난 당신이 t0을 발견하는 방법이 추측?) 이러한 속성을 찾기 위해 attributesstr를 사용하는 것입니다, 그러나 이것은 나를 위해 아무것도 결과. 그래서 인쇄 방법을 확인하기로 결정했습니다.

methods(print) 

얼핏 보여 주었다 저 print.bootsempls*

다음 단계이었다 신비를 해결

A single object matching ‘print.bootsempls’ was found 
It was found in the following places 
    registered S3 method for print from namespace semPLS 
    namespace:semPLS 
with value 

function (x, digits = 3, ...) 
{ 
    t <- x$t 
    t0 <- x$t0 
    result <- data.frame(Estimate = t0, Bias = colMeans(t, ...) - 
     t0, Std.Error = apply(t, 2, sd, ...)) 
    rownames(result) <- attr(t, "path") 
    cat("Call: ") 
    dput(x$call) 
    cat("\n") 
    print(result, digits = digits, ...) 
    invisible(x) 
} 
<environment: namespace:semPLS> 

: 결과

getAnywhere(print.bootsempls) 

. 그래서 apperantly print 메서드는 S.E.를 계산합니다. 어디에도 저장되지 않습니다. 당신은 단순히 다음을 수행해야합니다 당신이 보여 출력을 복제하고 (SE를 계산하기 위해도 및)를 data.frmae 오브젝트로 저장하려면 :

t <- mymodel_boot$t 
t0 <- mymodel_boot$t0 
result <- data.frame(Estimate = t0, 
        Bias = colMeans(t) - t0, 
        Std.Error = apply(t, 2, sd)) 
+1

어려운 작업을 수행하는 (+1) ...이 반환하지 않는 인쇄 방법은 사소한 자극입니다;),하지만 당신은'summary'에서 견적을 얻을 수 있다고 생각합니다. '? bootsempls'의 첫 번째 예제를 사용하면 이것이 효과가있는 것으로 보입니다. 'STR (들 <-.. 요약 (ecsiBoot)'그런 다음'의 $ table' – user20650

+1

user20650 @에서 견적을 얻을 수있는, 덕분에 나는 하나를 보지 못했지만, 내 생각은 그 OP 만 SE를 필요로하는 경우, 아마도'summary'를 계산하는 것보다'apply (mymodel_boot $ t, 2, sd)'를하는 것이 더 효율적일 것입니다. 요약문은 아마 많은 추가/불필요한 계산을 수행 할 것입니다. –

+1

어느 쪽이든 완벽하게 괜찮습니다. 너 둘 다. –