RxLinMod() 실행 "에서" 그룹에 대 한 몇 가지 재발 하는 기능이 있습니까?
rxLinMod()이이 기능이 없으면 직접. 이렇게 하려면 한 가지 방법은 각 고유 그룹 값에 대해 하나씩 별도 XDF 파일을 만들고 분할 XDF 파일 각각에 대해 다른 회귀를 계산 하는 루프에서 rxLinMod()를 호출 하 여 rxSplit() 함수를 사용 하는 것입니다.
예를 들어, http://packages.revolutionanalytics.com/datasets/ 여기 큰 항공 데이터를 사용 하 여
7 다른 XDF 파일을 각 요일 마다 하나씩 만듭니다. 다음 주에 하루에 선형 회귀 모델을 계산 해야 합니다.
rxGetVarInfo( "C:/MyRData/AirlineData87to08.xdf") Var 1: Year
22 factor levels: 1987 1988 1989 1990 1991 ... 2004 2005 2006 2007 2008
Var 2: Month
12 factor levels: January February March April May ... August September October November December
Var 3: DayofMonth
31 factor levels: 1 2 3 4 5 ... 27 28 29 30 31
Var 4: DayOfWeek
7 factor levels: Monday Tuesday Wednesday Thursday Friday Saturday Sunday
Var 5: DepTime, Type: integer, Low/High: (0, 29)
Var 6: CRSDepTime, Type: integer, Low/High: (0, 24)
Var 7: ArrTime, Type: integer, Low/High: (0, 29)
Var 8: CRSArrTime, Type: integer, Low/High: (0, 24)
Var 9: UniqueCarrier
29 factor levels: PS TW UA WN EA ... HA OH F9 YV 9E
Var 10: FlightNum
8160 factor levels: 1451 1453 1454 1455 1457 ... 9742 9743 6813 6913 6513
Var 11: TailNum
13537 factor levels: NA N7298U N7449U N7453U N7288U ... N516AS N763JB N766JB N75428 N75429
Var 12: ActualElapsedTime, Type: integer, Low/High: (-719, 1883)
Var 13: CRSElapsedTime, Type: integer, Low/High: (-1240, 1613)
Var 14: AirTime, Type: integer, Low/High: (-3818, 3508)
Var 15: ArrDelay, Type: integer, Low/High: (-1437, 2598)
Var 16: DepDelay, Type: integer, Low/High: (-1410, 2601)
Var 17: Origin
347 factor levels: SAN SFO BUR OAK LAX ... ROW GCC RKS MKG OTH
Var 18: Dest
352 factor levels: SFO RNO OAK BUR LAX ... PIR GCC RKS MKG OTH
Var 19: Distance, Type: integer, Low/High: (0, 4983)
Var 20: TaxiIn, Type: integer, Low/High: (0, 1523)
Var 21: TaxiOut, Type: integer, Low/High: (0, 3905)
Var 22: Cancelled, Type: logical, Low/High: (0, 1)
Var 23: CancellationCode
4 factor levels: A B C D
Var 24: Diverted, Type: logical, Low/High: (0, 1)
Var 25: CarrierDelay, Type: integer, Low/High: (0, 2580)
Var 26: WeatherDelay, Type: integer, Low/High: (0, 1510)
Var 27: NASDelay, Type: integer, Low/High: (-60, 1392)
Var 28: SecDelay, Type: integer, Low/High: (0, 533)
Var 29: LateAircraftDelay, Type: integer, Low/High: (0, 1407)
Var 30: Late, Type: logical, Low/High: (0, 1)
Var 31: LateNew
2 factor levels: 1 0
rxSplit(inData = "C:/MyRData/AirlineData87to08.xdf", splitByFactor = "DayOfWeek", outFilesBase = "")
lmresults <- vector("list", length = 0L)
for (i in 1:length(myinfo$DayOfWeek)[[5]]
{
dataFile = paste("C:/MyRData/AirlineData87to08-ByDay.DayOfWeek.",myfac[i],".xdf", sep = "")
lmresults[[i]] <- rxLinMod(ArrDelay ~ CRSDepTime, data = dataFile)
assign("lmresults", lmresults, pos = 1)
}