2016-09-19 3 views
2

RmiscsummarySE 함수를 사용하여 통합 문서의 각 열에 대한 데이터 요약을 생성하려고합니다. 워크 시트의 첫 번째 열은 그룹화 변수이며 다른 열을 반복하고 싶습니다. 나는 다음과 같은 코드를 사용하고 있습니다 :for 루프를 사용한 데이터 요약

library(Rmisc) 
for(i in 2:ncol(file)){ 
    sum<-summarySE(file, measurevar = file[,i], groupvars = file[1]) 
} 

을하지만 같은 오류가 계속 :

'Error in UseMethod("as.quoted") : 
no applicable method for 'as.quoted' applied to an object of class  
"data.frame"' 

내가 file[1]list 것을 알고 있고 vector해야하지만, unlist 더 많은 문제가 발생하여. 어떤 아이디어?

은 데이터 :

structure(list(Treatment = c("SKELE", "SKELE", "SKELE", "SKELE", 
"SKELE", "SKELE", "SKELE", "SKELE", "SKELE", "SKELE", "SKELE", 
"SKELE", "TISSUE", "TISSUE", "TISSUE", "TISSUE", "TISSUE", "TISSUE", 
"TISSUE", "TISSUE", "TISSUE", "TISSUE", "TISSUE", "TISSUE"), 
`% lipid in skeleton` = c(21.8706902567934, 31.1736436075643, 
62.2246234617322, 86.6248675033794, 46.5607971373041, 34.7532319115317, 
32.7686161366371, 6.73685660233744, 33.7111477556584, 48.8970450055359, 
54.3687328279357, 48.9086732773318, 78.1293097432066, 68.8263563924357, 
37.7753765382678, 13.3751324966206, 53.4392028626959, 65.2467680884683, 
67.2313838633629, 93.2631433976626, 66.2888522443416, 51.1029549944641, 
45.6312671720643, 51.0913267226682), `% ash in skeleton` = c(97.370981485225, 
98.6169174273543, 99.2417548180554, 99.1330769035889, 98.5523872323069, 
98.0077944962001, 97.7848485294277, 98.0738823145836, 98.1567971208113, 
98.8047064451889, 97.1790753033603, 98.7503991978965, 2.62901851477497, 
1.38308257264571, 0.75824518194458, 0.866923096411125, 1.44761276769314, 
1.99220550379987, 2.2151514705723, 1.92611768541643, 1.84320287918869, 
1.19529355481109, 2.82092469663973, 1.24960080210352), `% tissue in skeleton` = c(55.2224357342865, 
70.022864703591, 77.5880978578982, 83.1168129092154, 67.3012504898307, 
62.1455896726595, 64.2488985210074, 67.3089347382539, 59.9276126303114, 
70.5681668501146, 67.717146912379, 68.8185249866557, 44.7775642657135, 
29.977135296409, 22.4119021421018, 16.8831870907846, 32.6987495101694, 
37.8544103273405, 35.7511014789926, 32.6910652617461, 40.0723873696886, 
29.4318331498854, 32.2828530876211, 31.1814750133443)), class = "data.frame", row.names = c(NA, 
-24L), .Names = c("Treatment", "% lipid in skeleton", "% ash in skeleton", 
"% tissue in skeleton")) 

답변

1

우리는 summarySE을 적용 할 data.frame이 필요합니다. 지정의 '파일을'부분 집합 lapply, 컬럼의 순서하지만 루프를 사용하여

lapply(2:ncol(file), function(j) summarySE(file[c(1, j)], measurevar = 2, groupvars = 1)) 
#[[1]] 
# Treatment N  2  sd  se  ci 
#1  SKELE 12 42.38324 20.53836 5.928913 13.04945 
#2 TISSUE 12 57.61676 20.53836 5.928913 13.04945 

#[[2]] 
# Treatment N   2  sd  se  ci 
#1  SKELE 12 98.306052 0.6567242 0.1895799 0.4172626 
#2 TISSUE 12 1.693948 0.6567242 0.1895799 0.4172626 

#[[3]] 
# Treatment N  2  sd  se  ci 
#1  SKELE 12 67.83219 7.442443 2.148448 4.728703 
#2 TISSUE 12 32.16781 7.442443 2.148448 4.728703 

또는 우리는 OP의 방법을 사용하는 경우 (컬럼의 인덱스 기준) measurevargroupvars

lst <- vector("list", ncol(file)-1) 
for(i in 2:ncol(file)){ 
    lst[[i]] <- summarySE(file, measurevar = i, groupvars = 1) 
} 

lst 
우리가 대신 인덱스

for(i in 2:ncol(file)){ 
lst[[i]] <- summarySE(file, measurevar = names(file)[i], groupvars = names(file)[1]) 
} 

0의 names를 지정할 수 있습니다

OP 코드에서 measurevargroupvars은 열 이름 대신 열 값을 사용합니다.

관련 문제