FactoMineR을 사용하는 R을 사용한 다중 요인 분석 (MFA)

Question

FactoMineR의 MFA에 문제가 발생했습니다. 나는 토마토 식물에서 측정 된 물리적, 화학적 및 미생물 학적 연속 변수를 포함하는 데이터 세트로 작업하고 있는데, 2 가지 다른 치료법과 3 가지 시점에서 취한다. 나는이 같은 내 데이터를 수용 한 :

structure(list(row.names = structure(c(1L, 4L, 7L, 10L, 13L, 
16L), .Label = c("GBA1", "GBA2", "GBA3", "GBB1", "GBB2", "GBB3", 
"GBC1", "GBC2", "GBC3", "GBD1", "GBD2", "GBD3", "GBE1", "GBE2", 
"GBE3", "RWA1", "RWA2", "RWA3", "RWB1", "RWB2", "RWB3", "RWC1", 
"RWC2", "RWC3", "RWD1", "RWD2", "RWD3", "RWE1", "RWE2", "RWE3", 
"RWF1", "RWF2", "RWF3", "RWG1", "RWG2", "RWG3", "RWH1", "RWH2", 
"RWH3", "RWI1", "RWI2", "RWI3", "RWJ1", "RWJ2", "RWJ3"), class = "factor"), 
    Trt = structure(c(2L, 2L, 2L, 2L, 2L, 1L), .Label = c("Mineral", 
    "Organic"), class = "factor"), Status = structure(c(1L, 1L, 
    1L, 1L, 1L, 1L), .Label = c("H", "S"), class = "factor"), 
    Humidity = c(87.21704394, 80.29885961, 65.68047337, 85.9775641, 
    83.33333333, 85.98568282), pH = c(5.44, 5.94, 6.64, 6.19, 
    6.13, 5.45), Conductivity = c(837L, 867L, 752L, 871L, 699L, 
    406L), Nit.N = c(436.18, 433.92, 418.1, 458.78, 411.32, 167.24 
    ), Ammonia.N = c(3.8122, 2.6452, 1.945, 1.7116, 2.4896, 7.16 
    ), P = c(30.95, 15.2, 20.15, 16.1, 18.35, 48.2), K = c(135, 
    35, 95, 40, 145, 275), Ca = c(1287.5, 1427.5, 1610, 1570, 
    1640, 130), Mg = c(367.5, 575, 537.5, 532.5, 590, 42.5), 
    S = c(705L, 924L, 603L, 962L, 626L, 111L), Sodium = c(92.5, 
    170, 135, 127.5, 137.5, 35), Chlorides = c(15.1, 11.1, 15.4, 
    13.2, 13.8, 10.8), Fe = c(1.5, 2.2, 1.7, 2, 2.1, 3.1), Mn = c(1.1, 
    0.55, 0.7, 0.4, 0.65, 1.9), Rhizobium = c(0, 0, 0, 0, 0, 
    0), Total.bacteria = c(7207207.207, 5454545.455, 22727272.73, 
    18918918.92, 30630630.63, 64864864.86)), .Names = c("row.names", 
"Trt", "Status", "Humidity", "pH", "Conductivity", "Nit.N", "Ammonia.N", 
"P", "K", "Ca", "Mg", "S", "Sodium", "Chlorides", "Fe", "Mn", 
"Rhizobium", "Total.bacteria"), row.names = c(NA, 6L), class = "data.frame") 

난 후, 다른 (16)가 연속 범주 (처음 2)의 변수를 나누었다. 그러나 나는 2 개의 범주 형 변수를 따로 다루기를 원한다. 그래서 다음 코드를 작성했습니다 :

>res <- MFA(Oliver, group=c(1,1,3,11,2), type=c("n", "n","s", "s","s"),ncp=5,name.group=c("Sub","Stat", "Phys", "Chem", "Microbial")) 

그러나 작동하지 않는 것 같습니다.

>res=MFA(Oliver,group=c(2,16),type=c(rep("n",1),rep("s",1)),ncp=5,name.group=c("cat","cont")) 

를이 다른 : 따라서, 나는 다음과 같은 노력

>res=MFA(Oliver,group=c(2, 3, 11,2),type=c(rep("n",1),rep("s",3)), ncp=5,name.group=c("type","Phys", "Chem", "Microbial")) 

하지만 같은 문제 ("없습니다 편리 그룹 정의")를 가지고 있었다. 처음 2 개의 범주 그룹을 따로 유지하기 위해 할 수있는 일이 있습니까? 모델을 올바르게 실행하는 방법에 대한 귀하의 조언에 정말 감사드립니다!

행복을 빌며,

엠마

Answer 1

나는 문제가 모든 값에 "H"가 동일하기 때문에 변수하지 않은 당신의 변수 상태에서 오는 것이라고 생각합니다. 따라서이 "변수"로는 어떤 분석도 수행 할 수 없습니다. 이 열에 정보가 표시되지 않도록 설정할 수 있습니다. 그리고 나서 효과가 있습니다.

프랑소와