중첩 데이터 프레임

Question

기술적 인 문제가있어 직접 해결할 수 없습니다. mcmcglmm 패키지로 견적을 냈습니다. 결과 $ Sol에 의해 추정 된 사후 분포에 접근 할 수 있습니다. class()을 적용하면 객체의 클래스가 "mcmc"임을 알 수 있습니다. as.data.frame()을 사용하면 다른 데이터 프레임 (많은 다른 데이터 프레임을 포함하는 하나의 데이터 프레임)을 포함하는 중첩 데이터 프레임이됩니다. 모든 후방 분배의 모든 값과 (보조) 데이터 프레임의 이름을 rowname으로하여 하나의 데이터 프레임 (또는 오히려 벡터)을 생성하기 위해 주 데이터 프레임 내의 모든 데이터 프레임을 rbind()으로 보내고 싶습니다. ., 어떤 아이디어? 나는 모든 힌트에 감사 할 것입니다!

업데이트 : 이러한 모든 샘플링 체인을 사용하여 stackoverflow를 위해 유용한 데이터 세트를 생성하지 못했습니다. 이러한 데이터 세트는 항상 너무 큽니다. 당신이 나를 돕고 싶다면, 다음 (exemplaric) 모델 모든 체인이 저장되는 위치를 $ 솔 (AN MCMC 개체)입니다

require(MCMCglmm) 
data(PlodiaPO) 
result <- MCMCglmm(PO ~ plate + FSfamily, data = PlodiaPO, nitt = 50, thin = 2, burn = 10, verbose = FALSE) 

결과를 실행하는 데 고려하시기 바랍니다. 나는 모든 사후 분포의 모든 값과 변수 이름을 rownames로 갖는 벡터를 가지기 위해 (또는 중복 된 rownames가 추가 문자 벡터로 허용되지 않기 때문에) 모든 사슬을 역전시키고 싶다.

Answer 1

(MCMCglmm의 예제 코드를 사용하여) as.data.frame (model $ Sol)이 데이터 프레임의 데이터 프레임을 제공하는 예제를 구성 할 수 없습니다. 그래서 아마 간단한 대답 일지라도 그것을 매우 쉽게 확인할 수는 없습니다.

그렇다면 여기에 도움이 될만한 예가 있습니다. 자식 데이터 프레임에 동일한 이름이없는 경우에는 작동하지 않습니다.

# create a nested data.frame example to work on 
a.df <- data.frame(c1=runif(10),c2=runif(10)) 
b.df <- data.frame(c1=runif(10),c2=runif(10)) 
full.df <- data.frame(1:10) 
full.df$a <- a.df 
full.df$b <- b.df 
full.df <- full.df[,c("a","b")] 

# the solution 
res <- do.call(rbind,full.df) 

편집

좋아, 당신의 새로운 예를 들어,

require(MCMCglmm) 
data(PlodiaPO) 
result<- MCMCglmm(PO ~ plate + FSfamily, data=PlodiaPO,nitt=50,thin=2,burn=10,verbose=FALSE) 
melt(do.call(rbind,(as.data.frame(result$Sol))))