I can't get my head around the following odd thing:
I have a data.frame that I need to tidy up, and in this step I want to delete all entries that start with [ANONYMOUS] or a year in ().
Consider this example, a subset with two entries to illustrate my problem.
> print(allstacked2) #mysample- subset
values ind
711 [ANONYMOUS], 2010, COMMUNICATION. V1
1588 [ANONYMOUS], 2008, END UNLAWFUL ATTACKS. V2
1297 (1986) ACTIVE MEASURES, , U.S. DEPARTMENT OF STATE, AUGUST 1986 V1
3994 (1900) NEW YORK TIMES, P. 17. , 22 APRIL V3
> pat <- "(?:^)\\(.*"
> repl <- NA
> allstacked2$V3 <- gsub(pat, repl, allstacked2$values, perl=TRUE)
>
> pat <- "(^)\\[ANONYMOUS]"
> repl <- NA
> allstacked2$V3 <- gsub(pat, repl, allstacked2$V3, perl=TRUE)
>
> pat <- "[\\s]+"
> repl <- " "
> allstacked2$V3 <- gsub(pat, repl, allstacked2$V3, perl=TRUE)
>
> allstacked2[is.na(allstacked2)] = 'notavailable'
> allstacked2$V4 <- ifelse(allstacked2$V3==allstacked2$values,1,0) #to compare what has changed
> print(allstacked2)
values ind V3 V4
711 [ANONYMOUS], 2010, COMMUNICATION. V1 notavailable 0
1588 [ANONYMOUS], 2008, END UNLAWFUL ATTACKS. V2 [ANONYMOUS], 2008, END UNLAWFUL ATTACKS. 1
1297 (1986) ACTIVE MEASURES, , U.S. DEPARTMENT OF STATE, AUGUST 1986 V1 notavailable 0
3994 (1900) NEW YORK TIMES, P. 17. , 22 APRIL V3 (1900) NEW YORK TIMES, P. 17. , 22 APRIL 1
>
It seems to me that the problem has to do with the index colum, which I need to keep in order to unstack my data later.
I just don't get why. Thanks for any help here!!
EDIT:
Wiktor to the rescue - this is what happens when I use his pattern. It takes the stuff I want to replace, but won't replace the whole pattern. And I still wonder why it would return different results in my own regex.
> print(allstacked2)
values ind V3 V4
711 [ANONYMOUS], 2010, COMMUNICATION. V1 notavailable 0
1297 (1986) ACTIVE MEASURES, , U.S. DEPARTMENT OF STATE, AUGUST 1986 V1 notavailable 0
1588 [ANONYMOUS], 2008, END UNLAWFUL ATTACKS. V2 [ANONYMOUS], 2008, END UNLAWFUL ATTACKS. 1
3994 (1900) NEW YORK TIMES, P. 17. , 22 APRIL V3 (1900) NEW YORK TIMES, P. 17. , 22 APRIL 1
> allstacked2$V3 <- gsub("^[[:space:]]*(\\(\\d{4}\\)|\\[ANONYMOUS])[[:space:]]*", "", allstacked2$V3)
> allstacked2$V3 <- gsub("[[:space:]]+", " ", allstacked2$V3)
>
>
> print(allstacked2)
values ind V3 V4
711 [ANONYMOUS], 2010, COMMUNICATION. V1 notavailable 0
1297 (1986) ACTIVE MEASURES, , U.S. DEPARTMENT OF STATE, AUGUST 1986 V1 notavailable 0
1588 [ANONYMOUS], 2008, END UNLAWFUL ATTACKS. V2 , 2008, END UNLAWFUL ATTACKS. 1
3994 (1900) NEW YORK TIMES, P. 17. , 22 APRIL V3 NEW YORK TIMES, P. 17. , 22 APRIL 1
>
Edit: Here the dput!
> dput(allstacked2)
structure(list(values = c("[ANONYMOUS], 2010, COMMUNICATION.",
"(1986) ACTIVE MEASURES, , U.S. DEPARTMENT OF STATE, AUGUST 1986",
" [ANONYMOUS], 2008, END UNLAWFUL ATTACKS.", " (1900) NEW YORK TIMES, P. 17. , 22 APRIL"
), ind = structure(c(1L, 1L, 2L, 3L), .Label = c("V1", "V2",
"V3", "V4", "V5", "V6", "V7", "V8", "V9", "V10", "V11", "V12",
"V13", "V14", "V15", "V16", "V17", "V18", "V19", "V20", "V21",
"V22", "V23", "V24", "V25", "V26", "V27", "V28", "V29", "V30",
"V31", "V32", "V33", "V34", "V35", "V36", "V37", "V38", "V39",
"V40", "V41", "V42", "V43", "V44", "V45", "V46", "V47", "V48",
"V49", "V50", "V51", "V52", "V53", "V54", "V55", "V56", "V57",
"V58", "V59", "V60", "V61", "V62", "V63", "V64", "V65", "V66",
"V67", "V68", "V69", "V70", "V71", "V72", "V73", "V74", "V75",
"V76", "V77", "V78", "V79", "V80", "V81", "V82", "V83", "V84",
"V85", "V86", "V87", "V88", "V89", "V90", "V91", "V92", "V93",
"V94", "V95", "V96", "V97", "V98", "V99", "V100", "V101", "V102",
"V103", "V104", "V105", "V106", "V107", "V108", "V109", "V110",
"V111", "V112", "V113", "V114", "V115", "V116", "V117", "V118",
"V119", "V120", "V121", "V122", "V123", "V124", "V125", "V126",
"V127", "V128", "V129", "V130", "V131", "V132", "V133", "V134",
"V135", "V136", "V137", "V138", "V139", "V140", "V141", "V142",
"V143", "V144", "V145", "V146", "V147", "V148", "V149", "V150",
"V151", "V152", "V153", "V154", "V155", "V156", "V157", "V158",
"V159", "V160", "V161", "V162", "V163", "V164", "V165", "V166",
"V167", "V168", "V169", "V170", "V171", "V172", "V173", "V174",
"V175", "V176", "V177", "V178", "V179", "V180", "V181", "V182",
"V183", "V184", "V185", "V186", "V187", "V188", "V189", "V190",
"V191", "V192", "V193", "V194", "V195", "V196", "V197", "V198",
"V199", "V200", "V201", "V202", "V203", "V204", "V205", "V206",
"V207", "V208", "V209", "V210", "V211", "V212", "V213", "V214",
"V215", "V216", "V217", "V218", "V219", "V220", "V221", "V222",
"V223", "V224", "V225", "V226", "V227", "V228", "V229", "V230",
"V231", "V232", "V233", "V234", "V235", "V236", "V237", "V238",
"V239", "V240", "V241", "V242", "V243", "V244", "V245", "V246",
"V247", "V248", "V249", "V250", "V251", "V252", "V253", "V254",
"V255", "V256", "V257", "V258", "V259", "V260", "V261", "V262",
"V263", "V264", "V265", "V266", "V267", "V268", "V269", "V270",
"V271", "V272", "V273", "V274", "V275", "V276", "V277", "V278",
"V279", "V280", "V281", "V282", "V283", "V284", "V285", "V286",
"V287", "V288", "V289", "V290", "V291", "V292", "V293", "V294",
"V295", "V296", "V297", "V298", "V299", "V300", "V301", "V302",
"V303", "V304", "V305", "V306", "V307", "V308", "V309", "V310",
"V311", "V312", "V313", "V314", "V315", "V316", "V317", "V318",
"V319", "V320", "V321", "V322", "V323", "V324", "V325", "V326",
"V327", "V328", "V329", "V330", "V331", "V332", "V333", "V334",
"V335", "V336", "V337", "V338", "V339", "V340", "V341", "V342",
"V343", "V344", "V345", "V346", "V347", "V348", "V349", "V350",
"V351", "V352", "V353", "V354", "V355", "V356", "V357", "V358",
"V359", "V360", "V361", "V362", "V363", "V364", "V365", "V366",
"V367", "V368", "V369", "V370", "V371", "V372", "V373", "V374",
"V375", "V376", "V377", "V378", "V379", "V380", "V381", "V382",
"V383", "V384", "V385", "V386", "V387", "V388", "V389", "V390",
"V391", "V392", "V393", "V394", "V395", "V396", "V397", "V398",
"V399", "V400", "V401", "V402", "V403", "V404", "V405", "V406",
"V407", "V408", "V409", "V410", "V411", "V412", "V413", "V414",
"V415", "V416", "V417", "V418", "V419", "V420", "V421", "V422",
"V423", "V424", "V425", "V426", "V427", "V428", "V429", "V430",
"V431", "V432", "V433", "V434", "V435", "V436", "V437", "V438",
"V439", "V440", "V441", "V442", "V443", "V444", "V445", "V446",
"V447", "V448", "V449", "V450", "V451", "V452", "V453", "V454",
"V455", "V456", "V457", "V458", "V459", "V460", "V461", "V462",
"V463", "V464", "V465", "V466", "V467", "V468", "V469", "V470",
"V471", "V472", "V473", "V474", "V475", "V476", "V477", "V478",
"V479", "V480", "V481", "V482", "V483", "V484", "V485", "V486",
"V487", "V488", "V489", "V490", "V491", "V492", "V493", "V494",
"V495", "V496", "V497", "V498", "V499", "V500", "V501", "V502",
"V503", "V504", "V505", "V506", "V507", "V508", "V509", "V510",
"V511", "V512", "V513", "V514", "V515", "V516", "V517", "V518",
"V519", "V520", "V521", "V522", "V523", "V524", "V525", "V526",
"V527", "V528", "V529", "V530", "V531", "V532", "V533", "V534",
"V535", "V536", "V537", "V538", "V539", "V540", "V541", "V542",
"V543", "V544", "V545", "V546", "V547", "V548", "V549", "V550",
"V551", "V552", "V553", "V554", "V555", "V556", "V557", "V558",
"V559", "V560", "V561", "V562", "V563", "V564", "V565", "V566",
"V567", "V568", "V569", "V570", "V571", "V572", "V573", "V574",
"V575", "V576", "V577", "V578", "V579", "V580", "V581", "V582",
"V583", "V584", "V585", "V586", "V587", "V588", "V589", "V590",
"V591", "V592", "V593", "V594", "V595", "V596", "V597", "V598",
"V599", "V600", "V601", "V602", "V603", "V604", "V605", "V606",
"V607", "V608", "V609", "V610", "V611", "V612", "V613", "V614",
"V615", "V616", "V617", "V618", "V619", "V620", "V621", "V622",
"V623", "V624", "V625", "V626", "V627", "V628", "V629", "V630",
"V631", "V632", "V633", "V634", "V635", "V636", "V637", "V638",
"V639", "V640", "V641", "V642", "V643", "V644", "V645", "V646",
"V647", "V648", "V649", "V650", "V651", "V652", "V653", "V654",
"V655", "V656", "V657", "V658", "V659", "V660", "V661", "V662",
"V663", "V664", "V665", "V666", "V667", "V668", "V669", "V670",
"V671", "V672", "V673", "V674", "V675", "V676", "V677", "V678",
"V679", "V680", "V681", "V682", "V683", "V684", "V685", "V686",
"V687", "V688", "V689", "V690", "V691", "V692", "V693", "V694",
"V695", "V696", "V697", "V698", "V699", "V700", "V701", "V702",
"V703", "V704", "V705", "V706", "V707", "V708", "V709", "V710",
"V711", "V712", "V713", "V714", "V715", "V716", "V717", "V718",
"V719", "V720", "V721", "V722", "V723", "V724", "V725", "V726",
"V727", "V728", "V729", "V730", "V731", "V732", "V733", "V734",
"V735", "V736", "V737", "V738", "V739", "V740", "V741", "V742",
"V743", "V744", "V745", "V746", "V747", "V748", "V749", "V750",
"V751", "V752", "V753", "V754", "V755", "V756", "V757", "V758",
"V759", "V760", "V761", "V762", "V763", "V764", "V765", "V766",
"V767", "V768", "V769", "V770", "V771", "V772", "V773", "V774",
"V775", "V776", "V777", "V778", "V779", "V780", "V781", "V782",
"V783", "V784", "V785", "V786", "V787", "V788", "V789", "V790",
"V791", "V792", "V793", "V794", "V795", "V796", "V797", "V798",
"V799", "V800", "V801", "V802", "V803", "V804", "V805", "V806",
"V807", "V808", "V809", "V810", "V811", "V812", "V813", "V814",
"V815", "V816", "V817", "V818", "V819", "V820", "V821", "V822",
"V823", "V824", "V825", "V826", "V827", "V828", "V829", "V830",
"V831", "V832", "V833", "V834", "V835", "V836", "V837", "V838",
"V839", "V840", "V841", "V842", "V843", "V844", "V845", "V846",
"V847", "V848", "V849", "V850", "V851", "V852", "V853", "V854",
"V855", "V856", "V857", "V858", "V859", "V860", "V861", "V862",
"V863", "V864", "V865", "V866", "V867", "V868", "V869", "V870",
"V871", "V872", "V873", "V874", "V875", "V876", "V877", "V878",
"V879", "V880", "V881", "V882", "V883", "V884", "V885", "V886",
"V887", "V888", "V889", "V890", "V891", "V892", "V893", "V894",
"V895", "V896", "V897", "V898", "V899", "V900", "V901", "V902",
"V903", "V904", "V905", "V906", "V907", "V908", "V909", "V910",
"V911", "V912", "V913", "V914", "V915", "V916", "V917", "V918",
"V919", "V920", "V921", "V922", "V923", "V924", "V925", "V926",
"V927", "V928", "V929", "V930", "V931", "V932", "V933", "V934",
"V935", "V936", "V937", "V938", "V939", "V940", "V941", "V942",
"V943", "V944", "V945", "V946", "V947", "V948", "V949", "V950",
"V951", "V952", "V953", "V954", "V955", "V956", "V957", "V958",
"V959", "V960", "V961", "V962", "V963", "V964", "V965", "V966",
"V967", "V968", "V969", "V970", "V971", "V972", "V973", "V974",
"V975", "V976", "V977", "V978", "V979", "V980", "V981", "V982",
"V983", "V984", "V985", "V986", "V987", "V988", "V989", "V990",
"V991", "V992", "V993", "V994", "V995", "V996", "V997", "V998",
"V999", "V1000", "V1001", "V1002", "V1003", "V1004", "V1005",
"V1006", "V1007", "V1008", "V1009", "V1010", "V1011", "V1012",
"V1013", "V1014", "V1015", "V1016", "V1017", "V1018", "V1019",
"V1020", "V1021", "V1022", "V1023", "V1024", "V1025", "V1026",
"V1027", "V1028", "V1029", "V1030", "V1031", "V1032", "V1033",
"V1034", "V1035", "V1036", "V1037", "V1038", "V1039", "V1040",
"V1041", "V1042", "V1043", "V1044", "V1045", "V1046", "V1047",
"V1048", "V1049", "V1050", "V1051", "V1052", "V1053", "V1054",
"V1055", "V1056", "V1057", "V1058", "V1059", "V1060", "V1061",
"V1062", "V1063", "V1064", "V1065", "V1066", "V1067", "V1068",
"V1069", "V1070", "V1071", "V1072", "V1073", "V1074", "V1075",
"V1076", "V1077", "V1078", "V1079", "V1080", "V1081", "V1082",
"V1083", "V1084", "V1085", "V1086", "V1087", "V1088", "V1089",
"V1090", "V1091", "V1092", "V1093", "V1094", "V1095", "V1096",
"V1097", "V1098", "V1099", "V1100", "V1101", "V1102", "V1103",
"V1104", "V1105", "V1106", "V1107", "V1108", "V1109", "V1110",
"V1111", "V1112", "V1113", "V1114", "V1115", "V1116", "V1117",
"V1118", "V1119", "V1120", "V1121", "V1122", "V1123", "V1124",
"V1125", "V1126", "V1127", "V1128", "V1129", "V1130", "V1131",
"V1132", "V1133", "V1134", "V1135", "V1136", "V1137", "V1138",
"V1139", "V1140", "V1141", "V1142", "V1143", "V1144", "V1145",
"V1146", "V1147", "V1148", "V1149", "V1150", "V1151", "V1152",
"V1153", "V1154", "V1155", "V1156", "V1157", "V1158", "V1159",
"V1160", "V1161", "V1162", "V1163", "V1164", "V1165", "V1166",
"V1167", "V1168", "V1169", "V1170", "V1171", "V1172", "V1173",
"V1174", "V1175", "V1176", "V1177", "V1178", "V1179", "V1180",
"V1181", "V1182", "V1183", "V1184", "V1185", "V1186", "V1187",
"V1188", "V1189", "V1190", "V1191", "V1192", "V1193", "V1194",
"V1195", "V1196", "V1197", "V1198", "V1199", "V1200", "V1201",
"V1202", "V1203", "V1204", "V1205", "V1206", "V1207", "V1208",
"V1209", "V1210", "V1211", "V1212", "V1213", "V1214", "V1215",
"V1216", "V1217", "V1218", "V1219", "V1220", "V1221", "V1222",
"V1223", "V1224", "V1225", "V1226", "V1227", "V1228", "V1229",
"V1230", "V1231", "V1232", "V1233", "V1234", "V1235", "V1236",
"V1237", "V1238", "V1239", "V1240", "V1241", "V1242", "V1243",
"V1244", "V1245", "V1246", "V1247", "V1248", "V1249", "V1250",
"V1251", "V1252", "V1253", "V1254", "V1255", "V1256", "V1257",
"V1258", "V1259", "V1260", "V1261", "V1262", "V1263", "V1264",
"V1265", "V1266", "V1267", "V1268", "V1269", "V1270", "V1271",
"V1272", "V1273", "V1274", "V1275", "V1276", "V1277", "V1278",
"V1279", "V1280", "V1281", "V1282", "V1283", "V1284", "V1285",
"V1286", "V1287", "V1288", "V1289", "V1290", "V1291", "V1292",
"V1293", "V1294", "V1295", "V1296", "V1297", "V1298", "V1299",
"V1300", "index", "index2"), class = "factor")), .Names = c("values",
"ind"), row.names = c(711L, 1297L, 1588L, 3994L), class = "data.frame")
>
You have some whitespaces at the start of the input. You need to make sure you match those with [[:space:]]*:
allstacked2$V3 <- gsub("^[[:space:]]*(\\(\\d{4}\\)|\\[ANONYMOUS],)[[:space:]]*", "", allstacked2$V3)
allstacked2$V3 <- gsub("[[:space:]]+", " ", allstacked2$V3)
See the R demo online.
Details
^ - start of string
[[:space:]]* - 0+ whitespaces
(\\(\\d{4}\\)|\\[ANONYMOUS],) - one of the two alternatives:
\\(\\d{4}\\) - a (, 4 digits, )
| - or
\\[ANONYMOUS], - [ANONYMOUS], substring
[[:space:]]* - 0+ whitespaces
The [[:space:]] in a TRE regex will match Unicode whitespaces, so it is preferable than \s.
Related
I want to remove rows of exp.normal that ends with ".1" substring.
Then, I want to retain only rows where the rownames of exp.normal match another dataframe exp.kirp. I want to combine the dataframes column-wise.
My code returns the two dataframe with different row lengths.
exp.normal <- exp.normal[!is.infinite(rowSums(exp.normal)),]
exp.normal <- na.omit(exp.normal)
exp.normal <- exp.normal[!grepl('.1$', rownames(exp.normal)),]
exp.kirp.samp <- exp.kirp[rownames(exp.kirp) %in% rownames(exp.normal),]
exp.norm <- exp.normal[rownames(exp.normal) %in% rownames(exp.kirp.samp),]
Output:
> dim(exp.normal)
[1] 19947 32
> dim(exp.kirp)
[1] 12097 202
Traceback:
Error in cbind(...) : number of rows of matrices must match (see arg 2)
Example data:
dput(exp.norm)
structure(c(45.7005, 14525.5304, 2691.0051, 3648.1196, 3785.6462,
508.7428, 3386.262, 1189.0624, 375.0458, 1767.0259, 27.3361,
17196.2434, 2821.7784, 3730.9721, 8095.7046, 955.9156, 2899.3971,
1115.2977, 457.7995, 1821.7784, 45.3806, 19112.246, 3016.1901,
4261.1092, 9791.2504, 967.2683, 2105.4082, 795.0396, 419.566,
1941.4399, 50.9688, 14891.6723, 3558.9722, 3323.9259, 4598.989,
451.9671, 1407.3294, 410.278, 435.5518, 1647.4305, 108.6162,
9145.2729, 2705.0238, 2702.0338, 991.7092, 820.278, 3857.5328,
2015.6448, 332.5631, 1643.2139, 57.4382, 16482.7822, 2320.9426,
2881.0338, 5242.4173, 434.1923, 2427.6701, 985.4808, 369.4413,
1952.1095, 17.378, 16061.5305, 2530.1829, 4110.9756, 16517.6829,
1559.0915, 1917.0732, 607.9268, 558.2317, 1896.3415, 20.626,
16527.3632, 2925.5966, 4156.8309, 10196.4714, 635.9169, 2173.8045,
781.4243, 469.7868, 1987.3547, 9.4012, 21390.7875, 2698.0142,
4198.1283, 15663.0906, 1081.0287, 1256.9429, 591.6457, 585.8632,
1796.5457, 70.6982, 36833.3219, 2334.0517, 3118.3636, 3214.0421,
1119.1551, 2489.5293, 1626.5344, 372.9573, 1927.1034, 27.0724,
41553.7462, 2221.6031, 3961.8517, 7251.1787, 1841.8174, 1867.5525,
1033.0047, 366.9096, 2048.0069, 93.6277, 20812.3777, 3167.1196,
3068.6141, 1033.9674, 238.3764, 3262.2283, 927.9891, 368.2065,
1723.5054, 61.5543, 29531.5934, 3283.48, 3599.218, 531.7693,
224.8289, 2014.6628, 1123.1672, 321.6031, 1888.563, 28.7045,
14833.5053, 2878.5984, 3181.6764, 840.2164, 340.3734, 2266.8694,
1306.6222, 430.3965, 1614.8444, 32.3228, 20450.3576, 2638.1234,
3181.1024, 7459.6457, 791.6043, 2367.126, 982.2835, 444.2257,
1665.3543, 9.6095, 16474.3224, 2779.0965, 4393.5681, 14546.3247,
1160.0268, 1175.3446, 641.2711, 529.4793, 1929.9387, 33.5985,
13559.3697, 2940.1548, 3312.0163, 741.4452, 493.1579, 4346.1672,
1557.4905, 370.7226, 1643.4799, 26.9466, 26391.888, 2614.1732,
3922.3097, 11886.2642, 1467.3491, 2908.8364, 740.1575, 431.4961,
1638.1452, 51.9065, 20112.0953, 2840.2768, 3421.635, 1817.6237,
567.9274, 2244.2088, 867.8801, 385.6534, 1525.9852, 91.0713,
21614.0244, 2868.2575, 2628.6353, 408.8988, 486.2043, 1470.0472,
2422.3217, 404.9217, 1386.5275, 28.0177, 16412.3814, 2963.8349,
4044.6946, 7816.1037, 625.0426, 3087.001, 940.9758, 393.0399,
1762.8796, 45.7438, 30459.0464, 2345.8452, 3155.2153, 9420.2711,
1107.0516, 1614.3134, 818.9615, 484.3678, 1613.9445, 17.0691,
14674.209, 3037.4688, 3741.8818, 5809.3256, 688.5928, 4134.4713,
828.4763, 341.3822, 1690.2581, 44.1648, 12330.3918, 2558.2269,
3040.4316, 2229.2957, 686.3165, 4200.9857, 959.8567, 420.3641,
1570.6747, 17.6861, 19647.681, 2496.6604, 4041.9243, 6503.5269,
749.796, 2689.1762, 1000.3566, 423.1116, 1590.8995, 41.1018,
40847.6969, 2306.5671, 3025.5213, 3386.5546, 861.7896, 1909.7417,
1051.6651, 394.6467, 1680.361, 53.3787, 12392.7638, 2879.6926,
3503.1038, 7530.2986, 559.5507, 2631.9834, 796.9258, 437.1859,
1609.8138, 20.4024, 15847.5751, 2082.0815, 2966.7382, 1372.3176,
962.441, 1013.412, 810.0858, 478.0043, 1715.6652, 44.1964, 27787.4338,
2517.8366, 3247.7944, 3787.8788, 680.8477, 2439.9693, 955.1208,
431.5305, 2040.2762, 38.2464, 11097.4617, 3136.9164, 3515.104,
1173.7936, 411.0357, 3917.6148, 1790.5061, 351.9027, 1664.5743,
41.0752, 12219.0158, 3045.3933, 3753.4768, 371.6285, 612.8901,
2541.1433, 1990.7641, 445.4353, 1827.0053, 29.0631, 22662.2977,
2744.6112, 4235.6337, 13412.0641, 1127.8605, 1697.7296, 904.7316,
513.0433, 1923.9901), dim = c(10L, 32L), dimnames = list(c("A1BG",
"A2M", "AAMP", "AARS", "ABAT", "ABCA1", "ABCA2", "ABCA3", "ABCB7",
"ABCF1"), c("TCGA-BQ-5887-11A-01R-1965-07", "TCGA-DZ-6133-11A-01R-1965-07",
"TCGA-BQ-5884-11A-01R-1592-07", "TCGA-BQ-7044-11A-01R-1965-07",
"TCGA-BQ-5888-11A-01R-1592-07", "TCGA-BQ-7051-11A-02R-1965-07",
"TCGA-BQ-5879-11A-01R-1592-07", "TCGA-BQ-5882-11A-01R-1592-07",
"TCGA-BQ-5894-11A-01R-1592-07", "TCGA-DZ-6132-11A-01R-1965-07",
"TCGA-BQ-7045-11A-01R-1965-07", "TCGA-GL-A59R-11A-11R-A26U-07",
"TCGA-P4-A5E8-11A-12R-A28H-07", "TCGA-P4-A5ED-11A-11R-A28H-07",
"TCGA-BQ-7059-11A-01R-1965-07", "TCGA-BQ-5877-11A-01R-1592-07",
"TCGA-A4-A4ZT-11A-11R-A26U-07", "TCGA-BQ-5891-11A-01R-1592-07",
"TCGA-A4-A57E-11A-11R-A26U-07", "TCGA-BQ-7061-11A-01R-1965-07",
"TCGA-GL-6846-11A-01R-1965-07", "TCGA-BQ-5890-11A-01R-1592-07",
"TCGA-DZ-6131-11A-01R-1965-07", "TCGA-BQ-7055-11A-01R-1965-07",
"TCGA-B9-4115-11A-01R-1758-07", "TCGA-BQ-5875-11A-01R-1592-07",
"TCGA-BQ-7046-11A-01R-1965-07", "TCGA-GL-7966-11A-01R-2204-07",
"TCGA-DZ-6134-11A-01R-1965-07", "TCGA-GL-A9DE-11A-11R-A37K-07",
"TCGA-Y8-A8RY-11A-11R-A37K-07", "TCGA-BQ-5878-11A-01R-1592-07"
)))
> dput(exp.kirp)
structure(c(7.65342121905285, 14.3511850042327, 10.3737643425674,
10.0819596419255, 9.44832324553207, 5.36085937172008, 9.78880184184623,
10.3776687573505, 11.16757118884, 9.53872845925388, 9.59256168492636,
11.6467250199966, 9.64995723240483, 7.72893066783674, 7.8938008505495,
10.2355345297148, 4.7113572383413, 10.307474405626, 5.1591164988591,
6.82029258613417, 11.8163747078537, 10.6949102196217, 11.2422803547626,
8.59352109955772, 3.36085586100493, 9.85271363828624, 8.10771088195776,
11.3942834720292, 7.22380769346318, 10.1931428909004, 8.38100686984167,
10.16821542222, 7.89538927254103, 9.15267449482738, 10.3931817074315,
12.7214881212179, 10.1435311751243, 11.5928015984057, 9.61257956084687,
6.32433538607651, 9.96819071167712, 10.0874192149778, 10.14614948871,
9.68818857795196, 10.4090748876599, 10.6687443412299, 10.8848143975244,
11.5872726722286, 6.84628770347992, 10.2373774088459, 5.09389393824392,
12.4136523086721, 11.1918237390263, 10.1912122382252, 9.9623840324273,
6.22565668754941, 10.3398477765017, 10.3103072842012, 11.1287210937383,
9.12552428380532, 10.0220514952103, 11.5514020529076, 9.33154847897224,
7.15215907641367, 8.30013562492319, 10.2772713571793, 3.07471093384149,
9.90686083971644, 5.84486014570242, 7.85453119568792, 10.1727165187711,
11.3218051787965, 11.6207019570471, 8.14964399360302, 2.57221059331231,
9.770250762185, 8.2978654585063, 11.084297043801, 5.60196387752462,
8.88994408176988, 8.11137409958519, 10.6640911447013, 7.63002607924622,
8.98301980527478, 9.95231536364664, 9.50590581220822, 12.3983662709548,
9.74297479204394, 9.19670588373063, 5.03164662255944, 9.10229620506499,
10.4765987884097, 10.4203162794891, 9.17212793704954, 10.6963363394643,
10.517658858193, 10.6383042190661, 9.89639568127639, 4.96453263672961,
10.153415608879, 4.70168212029528, 10.1689338100564, 9.96839262629172,
9.87305770150294, 9.75535162798677, 6.170389794965, 10.238532641469,
9.94050095178643, 11.0690397931313, 9.18120494347434, 10.1026705963882,
12.2134644534149, 9.76558563096394, 7.17400985460863, 8.51390402420206,
10.6841464369114, 4.46450490913191, 10.2220147915752, 6.26732229936107,
8.34111937214434, 9.4563488702199, 11.3165579369515, 11.3692227748971,
7.56443152168505, 3.28865386360517, 9.3642795739989, 7.97047579449743,
4.39174710202512, 1.99320334863612, 9.09364811306327, 7.99728901996415,
11.0296194917623, 7.76961109298212, 9.33304596974955, 9.89296486273821,
13.2154969925887, 10.5733689401427, 10.0743768517653, 8.12622125684463,
4.20618996652279, 9.22288882589791, 10.7973270736212, 10.0356016423526,
9.9017250769121, 10.2083981309069, 11.6540022938525, 11.3739522748552,
8.46892945061225, 3.85569076543445, 10.0627068729804, 7.99645936536463,
13.6832285748428, 10.3714563849361, 10.4176870383992, 10.0652444551968,
5.98935248863472, 10.1079093507719, 11.2050505161752, 11.6645692817891,
9.14131578384568, 10.4504026669882, 11.7286541625055, 9.0418201112886,
7.18088382351455, 8.02888128245278, 10.4406656190814, 4.29335576875801,
9.41768714387146, 5.17869870850833, 6.88790782664042, 10.6188892562694,
10.3842123858415, 11.7658466970944, 7.8638220194434, 4.34994552290862,
9.22504206592907, 7.63249652719248, 10.3723102883537, 5.34994198553451,
9.48564780343169, 8.04610739496541, 10.4463535299365, 7.68066614055097,
9.04610766787879, 9.95556447653129, 11.8525448023587, 10.4640851987958,
10.4308625115247, 9.91159766583895, 5.37742225475561, 9.4487843137183,
9.99933682779041, 10.2770366214296, 9.0418201112886, 10.5578953924552,
10.4808099887997, 10.9764939558956, 10.4242897631097, 4.93490330712903,
10.9389453617464, 5.13719199914349, 12.0367193976262, 10.8202555636581,
10.3262700402849, 9.91216810777653, 5.52284042813955, 10.0653680664815,
10.5954686012028, 11.2355920880251, 9.10909645558853, 9.94472260124272,
11.5783590097537, 9.29504775099388, 7.53455423776441, 9.73144632203062,
10.3357853667922, 2.69972929338667, 11.0827581319988, 5.12711734972685,
7.20385010562547, 10.2823150418909, 9.98084618229907, 10.5982618153731,
8.87118417188149, 4.52284042813955, 8.94910562068286, 8.00490358291707,
11.0657044429112, 7.06536766254031, 8.7728136979997, 8.38298319671723,
10.0327018868117, 7.20326283373569, 9.47196961605789, 9.90466577812603,
9.49413435666035, 10.19403569598, 10.4724772111563, 11.1448925980171,
7.28004970633421, 9.06939981937328, 10.8754120177161, 10.2352930230415,
9.35902847180107, 10.7267739242029, 10.5529081501417, 11.2956215011293,
8.85725977829128, 4.55387039485638, 11.0056049014623, 6.95117512427229,
10.9656928148969, 10.5991523452113, 10.4556415168452, 9.46537845450025,
6.72337288854289, 10.0139441477751, 9.28408724134641, 11.4833270722276,
9.49617650025352, 10.4425557194467, 11.449658606754, 8.69606512512137,
6.7730033396692, 9.79804905625499, 10.6006257403548, 6.41043269493023,
9.33218159403931, 6.79787807595744, 3.66448284036468, 9.36947738932311,
12.1844424384667, 11.1649886100581, 8.30833493415398, 1.07949758402178,
8.50996853808744, 8.25444192773356, 11.0514188710826, 7.25444192773355,
10.2092692396731, 8.12117551878819, 10.7242737250263, 7.07511328163257,
9.24190737941017, 9.38899238984847, 11.5242722006977, 10.8138023505891,
9.39692881604932, 7.33690372047785, 5.99837831118633, 9.20104960820531,
10.271192322104, 10.3421982587672, 9.9583334671237, 10.8104108086575,
11.3588364978417, 10.3011032299907, 9.79976033549459, 4.03371371384131,
9.4502035528435, 3.61712213221935, 11.9975977242282, 9.91379851626213,
9.95999222715916, 9.90779350021794, 7.85831302551685, 10.3997246047534,
11.7402171909708, 11.7246448152361, 9.75167259046219, 9.68598603763708,
11.889881014998, 10.3276744540453, 8.04315308658688, 9.09516358798289,
10.1351204181247, 2.58736499093646, 10.6921133041207, 6.09727578619146,
10.2388193351428, 11.1387141417635, 9.60973686632844, 9.91029869099744,
7.0969238361025, 1.22478119439559, 7.85207251858217, 8.9132990393591,
11.4041748741886, 7.92325557174226, 11.9555213595621, 8.66299657133128,
9.99675988081849, 6.52046939465285, 9.19469924394784, 11.4453499696318,
13.7764385227361, 10.9582985177207, 9.73202730526663, 11.8793070895146,
5.674830656291, 9.30694800222568, 11.6556994451055, 10.4213132544326,
10.4922544904181, 10.5534738814693, 10.6996637140919, 10.2332275555211,
9.52531486289245, 6.37164075239763, 11.3566558699168, 3.15639661659767,
10.9506421354115, 10.8015070949819, 9.51027580369917, 10.0114476969219,
5.9232585434528, 10.1410206692489, 10.9093204661635, 11.1601119970792,
9.35889972298825, 9.89445931817658, 11.7097666903094, 10.3085136675886,
7.97175793055237, 9.39274501543653, 10.3787624907905, 6.94163806584425,
11.8483258820978, 5.28766436585954, 6.24271872469928, 10.0904144275087,
10.7658517928209, 10.6194247493219, 8.72532024089938, 6.12970957941129,
7.37932464409351, 8.97175821789909, 13.5632882613554, 6.31028194625214,
8.64822705939532, 9.18732138776524, 10.4159346076083, 6.46711478729057,
8.74638197005123, 9.80945126361765, 11.0458332956411, 12.6596471469673,
7.83425486655716, 10.9194011668402, 4.66760264911311, 8.44306371542639,
11.6735729859861, 10.592053792505, 10.0394481089966, 10.738081785684,
10.6804287993454, 10.4484290791818, 8.73413344315171, 6.43659995275074,
10.6685246199734, 8.20162111992077, 12.7954861179578, 10.0094938620897,
9.93400880935778, 9.75330735360058, 7.0082672553098, 9.74081003982032,
10.7382235152475, 11.6720072357516, 8.44494819471576, 9.99223068991282,
11.0703963722836, 9.47126368282115, 7.81421970385892, 9.50816775056871,
11.4983748936735, 6.23485957395859, 11.5903647810385, 6.26655044400704,
5.21798676951157, 10.0355914707199, 10.2338111446746, 11.8068330986116,
7.98850462265062, 2.26801511865652, 8.75403900144463, 8.99470954056443,
13.8050246222504, 9.39637273263616, 9.25174607188694, 8.59323012958596,
9.64227593422842, 7.73613603148148, 8.71848835304267, 11.3433845026852,
11.2390478373181, 10.8940677306605, 10.8800576375586, 10.6070670903737,
5.6912382257558, 9.46054395023954, 9.9164661944151, 9.35162700413489,
9.86387775885121, 10.0058120312308, 10.7880204671486, 10.3793997040699,
9.6038232649825, 4.26803009485916, 10.2782371039023, 6.26475409489153,
12.8198023802381, 11.7916439373724, 9.53606689182685, 10.3591574288036,
5.71406413888836, 10.0630462305102, 10.2195932783632, 11.455724780085,
8.42239448735832, 9.24745049099849, 12.1236985893816, 9.92518871297862,
7.89993146943836, 8.72944026793844, 10.4251666934438, 3.90670784266208,
9.20041190448898, 6.65613736826174, 10.2389705995139, 9.53818563335118,
11.7495320050765, 12.0646887131882, 8.44896824561834, 3.71406963555001,
9.8415805435622, 8.31270283399943, 5.24299494486063, 6.49167354977479,
9.98614168480194, 8.4582261368972, 9.59531026428392, 8.57655917226128,
8.95437954428976, 10.6495249285462, 14.9290954211312, 10.5422470205379,
9.17443338491354, 3.94233621005718, 5.67224428652735, 8.7869060376421,
10.8632316688763, 9.68413003969761, 10.0472843996015, 10.3973606067328,
10.2873877943623, 10.1185985217259, 9.54984491532867, 3.90670784266208,
8.67224499381939, 5.31084177712261, 13.1978705809921, 9.83862580682156,
9.54417855835282, 10.3479785104067, 6.95803180123955, 10.5275135616911,
10.7557094705532, 11.5723066760841, 9.32468663898682, 9.95090752477796,
11.895394953164, 10.1730449599892, 7.96865297090106, 8.25269807010003,
10.3946717952085, 3.87765591897716, 9.72523791794449, 6.55356786048423,
6.95803180123955, 10.5069909263905, 10.2877565830483, 10.5166977513248,
8.4726051481255, 4.38016145559557, 8.85875002019905, 8.67423888569003,
8.68286505123176, 5.31506259026344, 10.3333282396821, 9.31228692654063,
10.4142296315944, 7.85981553329836, 8.44751431678991, 11.0858929239232,
9.28987549215634, 11.1939401219005, 9.81036549777261, 8.34254452326977,
5.73563860631975, 9.00522312099486, 10.3688436329726, 9.94464471906402,
10.1522308353126, 10.7920437057173, 10.4456146711682, 10.5376544440661,
9.5106538108757, 6.00694202058687, 10.5239172935582, 6.07015744755466,
13.1209957412573, 10.9576837919022, 9.8945869439465, 9.23740931858995,
4.16326524153571, 10.523017510386, 10.532865308348, 10.619622367405,
9.45558963540716, 11.5136546699209, 11.209027561542, 8.50311772962086,
6.94462773667725, 9.43239470040379, 10.9762474117231, 3.74822505823549,
8.22396403818439, 5.28179467379438, 5.07015744755466, 9.40882053452885,
9.15195265459219, 11.7917134419519, 8.19669113197608, 4.74823042627325,
10.2400835226029, 7.76318049103673, 10.5252116829711, 5.07015744755466,
9.41356636990191, 8.69464935705204, 10.5448108847847, 6.4670479262821,
8.25602512546126, 10.0954826060933, 9.24807628369884, 10.2848014332547,
9.6103510748559, 7.88436713970017, 6.04591033802914, 7.98344684046261,
9.83392420760322, 10.426537121557, 8.77797787921888, 11.1932452952465,
10.889486091222, 10.329430966151, 9.26655583156462, 3.62269612816767,
10.5609425545211, 5.59521719360418, 13.9452588877727, 11.136365018286,
10.0338699199786, 9.01178744914112, 7.04236745871339, 9.77739489226773,
10.1944306108943, 11.0829417734315, 9.2483306049673, 10.2526052400654,
11.2392043296817, 9.78826123015354, 7.37009277188134, 9.28978490856306,
10.770835486621, 4.07112766139071, 8.79617819177588, 6.50274907789938,
4.95248267443048, 10.0124053945822, 10.8791588020481, 10.1400457232732,
9.64018220112663, 7.9855897990336, 8.44341851204054, 8.1851425389976,
11.3009567982918, 7.87032217794857, 9.12571179465448, 8.10396912418597,
10.0952450966319, 7.11418480239357, 9.3240403215499, 9.55790821677093,
12.5506935970456, 11.2890199649612, 9.74578287914873, 8.83302700041379,
6.80619256918401, 8.91595554395471, 10.4010993451482, 10.426361210739,
9.65687923171916, 10.3180541086121, 10.2515587448026, 10.4679560946132,
9.3710569683502, 5.63055399711973, 9.66161458306801, 3.49639843291562,
12.046947416625, 7.48250063297666, 10.3197777792499, 9.47365522397164,
6.95629018836077, 9.81146637438849, 9.67304984423569, 11.9690923038785,
9.99940969863312, 10.4101699932371, 12.2230772873117, 9.90523370112278,
7.10116469832952, 8.84409058546722, 11.2100883041446, 2.89400283481279,
12.1723192712031, 7.02813039834545, 7.76654547216671, 9.33005660517673,
9.71268261052322, 10.1695833227844, 7.83721098924632, 1.12849137851042,
7.14639353071401, 7.90194088386521, 10.2858186675831, 4.09194519835092,
9.82683604431664, 9.28935150377015, 10.3990006824976, 7.03632419457704,
9.25259305937555, 9.96009564626238, 13.6495813760682, 11.2087562804992,
8.62991100391626, 10.8093586506308, 4.36112967113519, 8.66452479510135,
10.680372183305, 10.4092425218598, 10.166842081886, 11.1148826551454,
11.5218621620828, 10.8830252212276, 8.50698349076537, 6.57472823565416,
10.546746035965, 7.53012147669141, 15.4681548425702, 9.79693513247838,
9.38793395115019, 8.98389796929759, 5.4990039364996, 9.68780798972713,
10.2221443965538, 11.2238327266469, 8.86631579791574, 9.0839648421511,
11.7688762222962, 9.43339244300495, 7.75707059749995, 7.30848010698361,
10.3788698682696, 5.10854118358023, 9.397938764198, 5.14397700761083,
5.54294517859488, 11.3477092783888, 10.1821779329165, 10.5068822633094,
7.95165772727702, 6.07899829224459, 9.17521394310747, 8.44018606106649,
10.2131062858423, 5.98655201209625, 9.76585079853419, 8.60644889080884,
9.71847037220014, 7.04446188071599, 9.27521498813589, 10.9901943954469,
11.5209125420426, 9.65147026322989, 10.6316810922306, 5.74458092333422,
8.26977703134258, 10.0141332026215, 9.03736135461269, 9.96719613631738,
9.93388860627356, 10.7248434525252, 10.2215812006449, 10.3621030533842,
11.4668080776894, 7.99184662847717, 9.33580121861008, 5.66820140316967,
11.9947030005322, 9.64794224178718, 10.1925126545798, 10.2674620603105,
7.5327081297796, 9.69101099541322, 11.2548391298699, 11.7689577047578,
9.19050300947907, 9.56006949444913, 11.2747943605218, 10.0095406093405,
8.00727063981421, 9.97974599674609, 10.3602156709497, 3.94458674388554,
12.5277668843542, 5.61301589155362, 10.2348465241763, 11.2408106474539,
10.7817980294952, 11.3929731417988, 8.56239088440517, 2.39709102017243,
8.66820104853111, 9.02757458085238, 12.6567500587125, 6.6095466615844,
11.4664830513311, 8.87780988126623, 10.3650994493445, 8.35292135920633,
7.93983129663749, 9.4812924824758, 12.99012924981, 10.4898544816469,
8.58540087815148, 10.8618786972288, 7.71067834148438, 8.93864071660805,
11.4031536223651, 10.7979767770093, 9.84830839236378, 9.47801724034487,
10.6468489299639, 10.3611048811716, 8.96812071185483, 5.42457950821431,
10.7176375280321, 5.99997294921438, 11.249476392449, 10.9207261360998,
10.3485991467441, 9.01994649367039, 6.53868559158395, 9.41478373733666,
10.432374028565, 11.0171355028297, 9.38017184771789, 10.5614675309011,
11.5465057188577, 9.17787519818794, 7.04065904752138, 11.0984136358699,
10.7076152049804, 2.37928824997193, 12.014117542008, 6.92057199739813,
4.98912992235826, 10.120606346187, 10.2084736088023, 10.4873799256012,
9.27923169328214, 1.52972101613862, 9.22178985762886, 7.95596194920281,
12.8045414410067, 7.92201464494711, 9.47902214483086, 7.4885399127133,
10.1807490909157, 6.70138409679205, 9.69962242189537, 10.1778751981879,
9.210584414225, 10.547619289895, 10.7420721677751, 11.6561438053443,
5.89426967756534, 9.09193991076902, 10.4509432609399, 10.2168981211661,
9.50697734395297, 10.1421973186731, 10.289252816822, 11.0579561962947,
8.44378357961471, 4.8801957289431, 11.4458754454719, 4.65344988200111,
12.0454913320053, 10.7230307474822, 9.54770815725829, 9.91346136214138,
6.41298301444827, 10.2791008470934, 11.2856294138809, 11.2550839389669,
9.39663175119314, 10.3644682533465, 12.3720956896476, 9.96448281969001,
7.26859346154833, 9.52846036210416, 10.6382351536235, 4.31017286300409,
11.0690011128305, 5.82641628715103, 6.80493638177069, 9.89786451034248,
9.94481833785949, 11.6774877172589, 9.01061202172265, 6.17266952106643,
9.5752544613585, 8.43194131360949, 11.0206656082034, 7.72110534236959,
10.5081439758347, 8.68994452434197, 10.8685929396885, 6.5444359272107,
8.68679113793189, 10.8525006396094, 11.1533795693109, 11.5873539730162,
8.99794551516942, 9.87102611145061, 3.93820109105871, 9.06141196830279,
10.3168170965086, 10.0968873453464, 9.78516864210074, 10.6089638880433,
10.6772891469749, 10.6027030400136, 10.5354892373923, 5.69623034672501,
10.4314703380043, 3.77147362287765, 11.9414284466395, 8.62911943026167,
9.47885474304884, 9.32603460205856, 7.52941181265447, 10.1897935554257,
10.6869529617348, 10.9824945325978, 9.54925041577186, 9.24348510543789,
12.6197561225204, 10.4229927556855, 7.78641143939821, 10.3518043181888,
10.2176900900617, 2.72462844562171, 13.9420477990178, 6.05669840185354,
8.84668667779751, 10.976770882955, 11.2310420929738, 11.9128532309879,
8.82897202130395, 2.92626515369416, 9.02841610083176, 9.70477480883974,
13.4396963975546, 6.34522248637767, 10.2353559766448, 8.78641143939821,
10.5578417013582, 6.64287771079513, 9.06181129266044, 10.6360273009065,
12.9150155941383, 11.9719694740862, 8.95056581144122, 10.0640399701312,
7.2355968744756, 9.21728961537032, 11.5561009223631, 9.95201054312933,
10.0229270264038, 10.7583335783419, 10.7564054571028, 10.3137246220883,
8.76465102394471, 5.79502400957859, 9.85443348382268, 5.2076547248541,
12.2695992772649, 10.7969615679355, 9.78033864553782, 9.39835742542764,
6.02727350425616, 9.5278442796299, 11.038246992221, 11.1457431696807,
9.00232447215102, 9.03646624386732, 12.062072892785, 10.0526651781329,
7.29576250712518, 7.86933843816656, 9.61223637364334, 3.32380238041871,
8.70355167900769, 5.73415681989236, 8.69454663545621, 11.3998442086989,
11.3703586767908, 12.3235644563292, 9.05367152133834, 10.196743500022,
9.91601721982026, 8.03544796058228, 12.1704251930224, 5.87162190579696,
8.07993112355535, 8.03136652207373, 9.91767640083016, 5.98915036355025,
9.10990689967612, 11.1609999554059, 11.726925238744, 12.2657872278982,
9.73588750953117, 11.4370505772681, 5.57906848944612, 9.27950936832482,
9.65861638042757, 9.17054127127384, 9.78579998093652, 10.5632561528303,
10.254342831578, 10.148100924125, 10.2726105411389, 4.18254929202166,
10.6197309312023), dim = c(50L, 20L), dimnames = list(c("A1BG",
"A2M", "A4GALT", "AAAS", "AACS", "AADAT", "AAGAB", "AAK1", "AAMP",
"AARS2", "AARSD1", "AARS", "AASDHPPT", "AASDH", "AASS", "AATF",
"AATK", "ABAT", "ABCA11P", "ABCA12", "ABCA1", "ABCA2", "ABCA3",
"ABCA5", "ABCA6", "ABCA7", "ABCB10", "ABCB1", "ABCB4", "ABCB6",
"ABCB7", "ABCB8", "ABCB9", "ABCC10", "ABCC1", "ABCC3", "ABCC4",
"ABCC5", "ABCC6", "ABCC9", "ABCD1", "ABCD3", "ABCD4", "ABCE1",
"ABCF1", "ABCF2", "ABCF3", "ABCG1", "ABCG2", "ABHD10"), c("TCGA.2K.A9WE.01A",
"TCGA.2Z.A9J1.01A", "TCGA.2Z.A9J3.01A", "TCGA.2Z.A9J5.01A", "TCGA.2Z.A9J6.01A",
"TCGA.2Z.A9J7.01A", "TCGA.2Z.A9J8.01A", "TCGA.2Z.A9JD.01A", "TCGA.2Z.A9JI.01A",
"TCGA.2Z.A9JJ.01A", "TCGA.2Z.A9JO.01A", "TCGA.2Z.A9JQ.01A", "TCGA.4A.A93W.01A",
"TCGA.4A.A93X.01A", "TCGA.4A.A93Y.01A", "TCGA.5P.A9JU.01A", "TCGA.5P.A9JY.01A",
"TCGA.5P.A9KE.01A", "TCGA.A4.7288.01A", "TCGA.A4.7583.01A")))
We may use intersect
nm1 <- intersect(row.names(exp.normal), row.names(exp.kirp))
exp.kirp.samp <- exp.kirp[nm1,]
exp.norm <- exp.normal[nm1,]
dim(exp.norm)
#[1] 8 32
dim(exp.kirp.samp)
#[1] 8 20
In the sample data showed, there are no duplicates for row names. It may be better to check for any duplicates with frequency count on the row names
> table(row.names(exp.normal))
A1BG A2M AAMP AARS ABAT ABCA2 ABCA3 ABCB7
1 1 1 1 1 1 1 1
> table(row.names(exp.kirp))
A1BG A2M A4GALT AAAS AACS AADAT AAGAB AAK1 AAMP AARS AARS2 AARSD1 AASDH AASDHPPT
1 1 1 1 1 1 1 1 1 1 1 1 1 1
AASS AATF AATK ABAT ABCA1 ABCA11P ABCA12 ABCA2 ABCA3 ABCA5 ABCA6 ABCA7 ABCB1 ABCB10
1 1 1 1 1 1 1 1 1 1 1 1 1 1
ABCB4 ABCB6 ABCB7 ABCB8 ABCB9 ABCC1 ABCC10 ABCC3 ABCC4 ABCC5 ABCC6 ABCC9 ABCD1 ABCD3
1 1 1 1 1 1 1 1 1 1 1 1 1 1
ABCD4 ABCE1 ABCF1 ABCF2 ABCF3 ABCG1 ABCG2 ABHD10
1 1 1 1 1 1 1 1
If there are duplicates, it is not clear about the logic for processing based on the post.
I was wondering if anyone knows of a way to get the dates to appear on the x axis of an stl plot.
res<- (stl(ts(data[,variable],frequency=52,start=as.Date(data[1,date]) ),s.window="per",robust=TRUE) )
plot(res)
However this is just producing some unicode labels at the bottom of the chart, I have also tried chaning as.Date to as.character however this didn't work.
dput of some data:
structure(c("2007-01-01", "2007-01-08", "2007-01-15", "2007-01-22",
"2007-01-29", "2007-02-05", "2007-02-12", "2007-02-19", "2007-02-26",
"2007-03-05", "2007-03-12", "2007-03-19", "2007-03-26", "2007-04-02",
"2007-04-09", "2007-04-16", "2007-04-23", "2007-04-30", "2007-05-07",
"2007-05-14", "2007-05-21", "2007-05-28", "2007-06-04", "2007-06-11",
"2007-06-18", "2007-06-25", "2007-07-02", "2007-07-09", "2007-07-16",
"2007-07-23", "2007-07-30", "2007-08-06", "2007-08-13", "2007-08-20",
"2007-08-27", "2007-09-03", "2007-09-10", "2007-09-17", "2007-09-24",
"2007-10-01", "2007-10-08", "2007-10-15", "2007-10-22", "2007-10-29",
"2007-11-05", "2007-11-12", "2007-11-19", "2007-11-26", "2007-12-03",
"2007-12-10", "2007-12-17", "2007-12-24", "2007-12-31", "2008-01-07",
"2008-01-14", "2008-01-21", "2008-01-28", "2008-02-04", "2008-02-11",
"2008-02-18", "2008-02-25", "2008-03-03", "2008-03-10", "2008-03-17",
"2008-03-24", "2008-03-31", "2008-04-07", "2008-04-14", "2008-04-21",
"2008-04-28", "2008-05-05", "2008-05-12", "2008-05-19", "2008-05-26",
"2008-06-02", "2008-06-09", "2008-06-16", "2008-06-23", "2008-06-30",
"2008-07-07", "2008-07-14", "2008-07-21", "2008-07-28", "2008-08-04",
"2008-08-11", "2008-08-18", "2008-08-25", "2008-09-01", "2008-09-08",
"2008-09-15", "2008-09-22", "2008-09-29", "2008-10-06", "2008-10-13",
"2008-10-20", "2008-10-27", "2008-11-03", "2008-11-10", "2008-11-17",
"2008-11-24", "2008-12-01", "2008-12-08", "2008-12-15", "2008-12-22",
"2008-12-29", "2009-01-05", "2009-01-12", "2009-01-19", "2009-01-26",
"2009-02-02", "2009-02-09", "2009-02-16", "2009-02-23", "2009-03-02",
"2009-03-09", "2009-03-16", "2009-03-23", "2009-03-30", "2009-04-06",
"2009-04-13", "2009-04-20", "2009-04-27", "2009-05-04", "2009-05-11",
"2009-05-18", "2009-05-25", "2009-06-01", "2009-06-08", "2009-06-15",
"2009-06-22", "2009-06-29", "2009-07-06", "2009-07-13", "2009-07-20",
"2009-07-27", "2009-08-03", "2009-08-10", "2009-08-17", "2009-08-24",
"2009-08-31", "2009-09-07", "2009-09-14", "2009-09-21", "2009-09-28",
"2009-10-05", "2009-10-12", "2009-10-19", "2009-10-26", "2009-11-02",
"2009-11-09", "2009-11-16", "2009-11-23", "2009-11-30", "2009-12-07",
"2009-12-14", "2009-12-21", "2009-12-28", "2010-01-04", "2010-01-11",
"2010-01-18", "2010-01-25", "2010-02-01", "2010-02-08", "2010-02-15",
"2010-02-22", "2010-03-01", "2010-03-08", "2010-03-15", "2010-03-22",
"2010-03-29", "2010-04-05", "2010-04-12", "2010-04-19", "2010-04-26",
"2010-05-03", "2010-05-10", "2010-05-17", "2010-05-24", "2010-05-31",
"2010-06-07", "2010-06-14", "2010-06-21", "2010-06-28", "2010-07-05",
"2010-07-12", "2010-07-19", "2010-07-26", "2010-08-02", "2010-08-09",
"2010-08-16", "2010-08-23", "2010-08-30", "2010-09-06", "2010-09-13",
"2010-09-20", "2010-09-27", "2010-10-04", "2010-10-11", "2010-10-18",
"2010-10-25", "2010-11-01", "2010-11-08", "2010-11-15", "2010-11-22",
"2010-11-29", "2010-12-06", "2010-12-13", "2010-12-20", "2010-12-27",
"2011-01-03", "2011-01-10", "2011-01-17", "2011-01-24", "2011-01-31",
"2011-02-07", "2011-02-14", "2011-02-21", "2011-02-28", "2011-03-07",
"2011-03-14", "2011-03-21", "2011-03-28", "2011-04-04", "2011-04-11",
"2011-04-18", "2011-04-25", "2011-05-02", "2011-05-09", "2011-05-16",
"2011-05-23", "2011-05-30", "2011-06-06", "2011-06-13", "2011-06-20",
"2011-06-27", "2011-07-04", "2011-07-11", "2011-07-18", "2011-07-25",
"2011-08-01", "2011-08-08", "2011-08-15", "2011-08-22", "2011-08-29",
"2011-09-05", "2011-09-12", "2011-09-19", "2011-09-26", "2011-10-03",
"2011-10-10", "2011-10-17", "2011-10-24", "2011-10-31", "2011-11-07",
"2011-11-14", "2011-11-21", "2011-11-28", "2011-12-05", "2011-12-12",
"2011-12-19", "2011-12-26", "442573", "452832", "452785", "459228",
"479509", "477631", "465619", "462001", "485567", "462381", "456059",
"457401", "474094", "468766", "456945", "539126", "545640", "511801",
"486619", "484430", "481428", "470622", "479677", "486755", "477617",
"483656", "497479", "493436", "480080", "481527", "516029", "532349",
"503939", "472171", "461550", "504532", "489715", "480032", "488376",
"470772", "467395", "488155", "455712", "474456", "471237", "482943",
"459320", "456956", "465056", "461930", "441201", "451255", "464508",
"449199", "455411", "476323", "521761", "513416", "521070", "497596",
"485461", "485593", "461148", "429938", "441207", "459484", "462099",
"469285", "454395", "456729", "469251", "517727", "526719", "477985",
"484538", "469766", "472399", "481162", "479039", "481071", "485068",
"462108", "459079", "452410", "488574", "502166", "504323", "514529",
"524206", "504315", "462935", "461899", "461551", "455491", "456857",
"443314", "454586", "458943", "450555", "454311", "442808", "440126",
"414876", "415787", "413352", "423864", "448319", "415440", "431948",
"433313", "448509", "436400", "454154", "454183", "447735", "452220",
"451433", "472808", "446767", "426595", "463693", "435673", "452704",
"456828", "465069", "448685", "457353", "443859", "463972", "480139",
"493342", "500982", "529602", "527365", "512005", "484585", "455007",
"470006", "477522", "443956", "459038", "488877", "466669", "476242",
"470862", "457298", "466438", "454062", "460216", "466645", "445113",
"457255", "451553", "451504", "447991", "435100", "416140", "481390",
"489041", "496984", "486245", "478191", "455201", "448085", "459097",
"480859", "490248", "462523", "489755", "468391", "460229", "481276",
"472845", "481099", "476435", "487314", "475043", "476847", "461928",
"483488", "479379", "456732", "461538", "480773", "471101", "459898",
"482129", "464356", "459420", "457850", "456764", "438152", "467928",
"464732", "458767", "470256", "449488", "433263", "428548", "435099",
"429163", "424251", "434723", "425841", "399498", "418791", "405051",
"385037", "425144", "418296", "397644", "414283", "431907", "429117",
"424862", "439664", "432791", "443588", "434985", "442418", "445760",
"449290", "451412", "456247", "444372", "441390", "458192", "456435",
"450670", "447609", "439083", "464513", "462784", "439423", "450857",
"442374", "447753", "440207", "435254", "430841", "437233", "426523",
"430127", "431305", "470244", "508878", "511064", "504182", "462076",
"452218", "426535", "436892", "459008", "441449", "438783", "427497",
"432275", "436745", "423068", "429574", "416074"), .Dim = c(261L,
2L))
The problem is with how you have created the time series object.
?ts will give you
start : the time of the first observation. Either a single number or a
vector of two integers, which specify a natural time unit and a
(1-based) number of samples into the time unit.
Since your data starts at 2007-01-01 and is a weekly data, you have to specify c(2007,1), Note that the "1" here denotes that data start from the 1st week of 2007, as the data is weekly. Here the frequency should be set to 52 because, after 52 cycles/weeks a year will change from 2007 to 2008.
This should work
res<- (stl(ts(data[,2],frequency=52,start=c(2007,1 )),s.window="per",robust=TRUE) )
I have two dfs - benchmark return and allreturns (samples below), and I'd like to create a third df that returns a subtracted value allreturns-benchmarkreturns, matched by a "startmonth" value for each row of allreturns, as well as the colname (V1:V12) that is consistent in both dfs.
This problem is akin to an index match with two criteria in Excel, where I'd match by startmonth on the row side and "V month" in the column side (colnames here).
So for instance, if in allreturns I had a January 2005 startmonth and the value for that row in V3 was 5, and the benchmarkreturn for row January 2005 and column V3 was 3, the value for that row of allreturns in the new df would be 3. And so on.
The actual allreturns dataframe is quite large, so I am concerned about merging before subtracting, and was wondering if that is the most viable method or if I can subtract the elements in one step.
Thank you for your help!
All returns:
structure(list(startmonth = structure(c(2012, 2012.08333333333,
2012.16666666667, 2012.25, 2012.33333333333, 2012.41666666667,
2012.5, 2012.58333333333, 2012.66666666667, 2012.75), class = "yearmon"),
V1 = c(0.0284247303, 0.0010087300099999, 0.0022169700300001,
-0.01468004942, 0.00102045998000011, 0.00591265976000011,
0.0199267303899999, 0.0417483615900001, 0.00335606009, 0.00334482998000007
), V2 = c(0.0294621331884797, 0.00322793636420049, -0.012495624619603,
-0.0136745698429375, 0.0069391533726606, 0.0259572101269254,
0.0625070003256281, 0.0452445316901551, 0.00671211552040374,
0.0314807661806089), V3 = c(0.0317444198847785, -0.0114994993211506,
-0.0114879159244523, -0.00784276316178301, 0.0270041584010525,
0.0687892427111718, 0.0660728376647666, 0.0487406969361834,
0.0349424776135903, 0.0430893809414579), V4 = c(0.0165983608120608,
-0.0104907741199978, -0.00564318030266497, 0.0119276866011795,
0.0698798993604133, 0.0723761636332563, 0.0696386700530516,
0.0781496303938529, 0.0465900515831383, 0.0514300069951692
), V5 = c(0.0176357587550031, -0.00464014273798818, 0.0141710999549016,
0.0541740095643377, 0.0734704805917499, 0.0759630795752142,
0.0996336273987892, 0.090283471493128, 0.054958669279074,
0.0744177901533509), V6 = c(0.0236526927561311, 0.0151941247787009,
0.056511081749947, 0.057711880885752, 0.0770610568378784,
0.106135386898253, 0.112009256278245, 0.099001465266642,
0.0780236008188706, 0.0925232124743112), V7 = c(0.0440507439776798,
0.0575768161840056, 0.0600567964260508, 0.0612497472951408,
0.107264153793045, 0.118584188660589, 0.120900971163225,
0.123029319908411, 0.0961897860250946, 0.0862168299189658
), V8 = c(0.0876381519555687, 0.06112610752891, 0.0636025061792393,
0.0910094612943682, 0.119725659059378, 0.127528477230367,
0.145407622386291, 0.141953913583335, 0.0898622388709431,
0.0799104502736558), V9 = c(0.0912883309497081, 0.0646753939459335,
0.0934281965725043, 0.103288031047544, 0.128679074917691,
0.152180027943971, 0.164709321337116, 0.135362201702614,
0.0835346946365931, 0.076045245502748), V10 = c(0.0949385048758928,
0.0945311704081673, 0.105733987546834, 0.11211001031858,
0.153355781552244, 0.171595851273483, 0.157986258190939,
0.128770492863592, 0.0796565180062141, 0.069524282329458),
V11 = c(0.125642923731963, 0.10684937460315, 0.114575524881519,
0.136424461680214, 0.172791417952441, 0.164833036942378,
0.151263198147072, 0.124730408712938, 0.0731136700886696,
0.0680519003357034), V12 = c(0.138311269558968, 0.115699830661878,
0.138943880653118, 0.155574781977547, 0.166021702439959,
0.158070225731925, 0.147142608328795, 0.117914407614703,
0.0716363466924688, 0.086773098425986)), .Names = c("startmonth",
"V1", "V2", "V3", "V4", "V5", "V6", "V7", "V8", "V9", "V10",
"V11", "V12"), class = "data.frame", row.names = c(NA, 10L))
Benchmark returns:
structure(list(startmonth = structure(c(2005, 2005.08333333333,
2005.16666666667, 2005.25, 2005.33333333333, 2005.41666666667,
2005.5, 2005.58333333333, 2005.66666666667, 2005.75, 2005.83333333333,
2005.91666666667, 2006, 2006.08333333333, 2006.16666666667, 2006.25,
2006.33333333333, 2006.41666666667, 2006.5, 2006.58333333333,
2006.66666666667, 2006.75, 2006.83333333333, 2006.91666666667,
2007, 2007.08333333333, 2007.16666666667, 2007.25, 2007.33333333333,
2007.41666666667, 2007.5, 2007.58333333333, 2007.66666666667,
2007.75, 2007.83333333333, 2007.91666666667, 2008, 2008.08333333333,
2008.16666666667, 2008.25, 2008.33333333333, 2008.41666666667,
2008.5, 2008.58333333333, 2008.66666666667, 2008.75, 2008.83333333333,
2008.91666666667, 2009, 2009.08333333333, 2009.16666666667, 2009.25,
2009.33333333333, 2009.41666666667, 2009.5, 2009.58333333333,
2009.66666666667, 2009.75, 2009.83333333333, 2009.91666666667,
2010, 2010.08333333333, 2010.16666666667, 2010.25, 2010.33333333333,
2010.41666666667, 2010.5, 2010.58333333333, 2010.66666666667,
2010.75, 2010.83333333333, 2010.91666666667, 2011, 2011.08333333333,
2011.16666666667, 2011.25, 2011.33333333333, 2011.41666666667,
2011.5, 2011.58333333333, 2011.66666666667, 2011.75, 2011.83333333333,
2011.91666666667, 2012, 2012.08333333333, 2012.16666666667, 2012.25,
2012.33333333333, 2012.41666666667, 2012.5, 2012.58333333333,
2012.66666666667, 2012.75, 2012.83333333333, 2012.91666666667,
2013, 2013.08333333333, 2013.16666666667, 2013.25, 2013.33333333333,
2013.41666666667, 2013.5, 2013.58333333333, 2013.66666666667,
2013.75, 2013.83333333333, 2013.91666666667, 2014, 2014.08333333333,
2014.16666666667, 2014.25, 2014.33333333333, 2014.41666666667,
2014.5, 2014.58333333333, 2014.66666666667, 2014.75, 2014.83333333333
), class = "yearmon"), V1 = c(0.83740290789865, 0.175849379040301,
0.115745179820806, 0.588703967630863, 0.663440374424681, 0.122854199260473,
0.329111734172329, 0.386456901207566, 0.328031520592049, 0.91476626554504,
0.838957223575562, 0.482753940392286, 0.203118846053258, 0.832898694789037,
0.569874424720183, 0.85476457583718, 0.873173910891637, 0.599100088467821,
0.990275882417336, 0.360726963961497, 0.366468257969245, 0.213221889222041,
0.736634878907353, 0.770827260101214, 0.78438491653651, 0.443004491971806,
0.0431681957561523, 0.362984816078097, 0.933778994018212, 0.890268561197445,
0.0288882539607584, 0.704298130935058, 0.0533608365803957, 0.965311451349407,
0.797033908311278, 0.702818523859605, 0.615030416054651, 0.783928481163457,
0.0931849081534892, 0.693762929877266, 0.709550492931157, 0.527024435112253,
0.463098735082895, 0.975416449131444, 0.418448116164654, 0.716492718551308,
0.900618054671213, 0.937315494287759, 0.371745373355225, 0.876447483897209,
0.638362225377932, 0.230460224673152, 0.84510417538695, 0.709547182312235,
0.135460949270055, 0.86736886552535, 0.0408105747774243, 0.59634031355381,
0.378888477338478, 0.400274467887357, 0.73448179056868, 0.0201477019581944,
0.786700234282762, 0.145882265642285, 0.722347423899919, 0.495645943330601,
0.8204965388868, 0.944398637628183, 0.920839294558391, 0.244003586936742,
0.770501866471022, 0.608119423966855, 0.599260926479474, 0.316019677789882,
0.594015010166913, 0.862004647962749, 0.248981321696192, 0.251351157668978,
0.86709617683664, 0.477806584443897, 0.308428246993572, 0.930116105359048,
0.820579398889095, 0.407806032337248, 0.00576947652734816, 0.205415812321007,
0.503058285918087, 0.421452875249088, 0.171408428112045, 0.927357386332005,
0.736080843256786, 0.0634949475061148, 0.628076422959566, 0.825837341602892,
0.860992948058993, 0.637153514893726, 0.496678973315284, 0.532001776387915,
0.742080342955887, 0.279888727236539, 0.702802252722904, 0.20291490922682,
0.761654698755592, 0.413414426380768, 0.0555722850840539, 0.125917585333809,
0.545490033458918, 0.19256273820065, 0.242046363884583, 0.0174444266594946,
0.380734569625929, 0.735217718640342, 0.346396475099027, 0.997578310314566,
0.544119880301878, 0.610004446236417, 0.643721028929576, 0.0819557739887387,
0.666767429094762), V2 = c(0.325795059790835, 0.0280552506446838,
0.563847373938188, 0.968374367104843, 0.0397681496106088, 0.409115779446438,
0.568448322825134, 0.498488057404757, 0.725102471420541, 0.745882571907714,
0.956671225838363, 0.732607456855476, 0.723114193649963, 0.715586860431358,
0.740428047254682, 0.447626604698598, 0.136078626848757, 0.343137193005532,
0.994878244586289, 0.569102957844734, 0.264083825284615, 0.0172849236987531,
0.902412769151852, 0.237230452476069, 0.342157111503184, 0.940635218983516,
0.355561858741567, 0.0913543966598809, 0.730372012127191, 0.424820594489574,
0.862009244738147, 0.673514985479414, 0.0558020728640258, 0.769041493767872,
0.429111589910462, 0.814037810778245, 0.485316230216995, 0.407698386581615,
0.0245183343067765, 0.38149883877486, 0.939248957671225, 0.206260631326586,
0.773903579683974, 0.955869939411059, 0.930366629268974, 0.497170565184206,
0.232823575846851, 0.119068749016151, 0.188956700265408, 0.769491335609928,
0.487886165967211, 0.249603169038892, 0.61157290590927, 0.967520525446162,
0.669223197968677, 0.0509379308205098, 0.630447809584439, 0.93567297863774,
0.296445045853034, 0.231306056724861, 0.664361900882795, 0.570044811815023,
0.988849612651393, 0.96225642086938, 0.567738547222689, 0.562166472664103,
0.382391084451228, 0.87047664122656, 0.553916300181299, 0.348980325739831,
0.861658359179273, 0.824328954331577, 0.602889677742496, 0.108061951119453,
0.510826006531715, 0.787226195679978, 0.209260584786534, 0.467020642245188,
0.916463258676231, 0.810527657391503, 0.259206298971549, 0.383183149620891,
0.399636579910293, 0.517378240358084, 0.277383098145947, 0.0540329110808671,
0.334274391178042, 0.762725146021694, 0.624413190409541, 0.664978531654924,
0.877504705451429, 0.133810454979539, 0.303127287654206, 0.401904000202194,
0.333468962227926, 0.880315525457263, 0.693340804660693, 0.104616170749068,
0.726646555587649, 0.787352371029556, 0.805738570401445, 0.842183391330764,
0.203002500347793, 0.4008568006102, 0.764609909849241, 0.68458701716736,
0.826388106215745, 0.458343503531069, 0.681572974193841, 0.816049969289452,
0.99422895279713, 0.586194961098954, 0.106009811395779, 0.564460329711437,
0.763485615840182, 0.641962635098025, 0.556110577424988, 0.624630474718288,
0.0865604616701603), V3 = c(0.103272917214781, 0.145110326353461,
0.139097063103691, 0.704820184968412, 0.499401112552732, 0.183894460555166,
0.142426387872547, 0.654179619625211, 0.231400778284296, 0.07779724849388,
0.639861593721434, 0.329660887364298, 0.207785484148189, 0.349164958344772,
0.918049042113125, 0.561889169970527, 0.801518984371796, 0.374053293140605,
0.279966498259455, 0.211821462959051, 0.988277795957401, 0.208132355008274,
0.358667319407687, 0.489217024296522, 0.235113901318982, 0.0665953925345093,
0.682081600185484, 0.923956993501633, 0.951206840341911, 0.153028624830768,
0.457171397982165, 0.497927478514612, 0.464426794555038, 0.700641778064892,
0.350401868345216, 0.239867100957781, 0.495585883036256, 0.659481293288991,
0.121178397443146, 0.966434765374288, 0.781230272026733, 0.191542416810989,
0.144150232896209, 0.649679786758497, 0.990280882688239, 0.182715938193724,
0.990455492865294, 0.0560185567010194, 0.258451806847006, 0.745103346649557,
0.389899141388014, 0.0482800761237741, 0.0293514989316463, 0.393075271975249,
0.53582366858609, 0.522070667007938, 0.802251521963626, 0.292497304268181,
0.523711683694273, 0.606568831950426, 0.198962268885225, 0.391696797916666,
0.059201990487054, 0.0645288161467761, 0.808837213553488, 0.0917617143131793,
0.698996715713292, 0.330968038411811, 0.843029006384313, 0.9479373190552,
0.514315341366455, 0.325804576743394, 0.3519700136967, 0.537765094777569,
0.871437344467267, 0.667971923248842, 0.0599088978487998, 0.560071154031903,
0.902071484364569, 0.193852530093864, 0.629741749959067, 0.0646166354417801,
0.238403898198158, 0.737684201914817, 0.0984734564553946, 0.314926577033475,
0.667433348251507, 0.200710391392931, 0.510746134677902, 0.0941306182648987,
0.408070910954848, 0.349233734188601, 0.0384121341630816, 0.469132867408916,
0.968564240029082, 0.250593126285821, 0.0388190811499953, 0.301829227479175,
0.62604419933632, 0.981997888535261, 0.699610964162275, 0.319871086394414,
0.0948076455388218, 0.0173410377465189, 0.733540053712204, 0.304378144675866,
0.593588132876903, 0.85377266141586, 0.272565411170945, 0.151829193811864,
0.951305341674015, 0.21787527250126, 0.728856807108968, 0.932554795406759,
0.530237913830206, 0.146221859846264, 0.220554714323953, 0.133143884129822,
0.845931079704314), V4 = c(0.828090658877045, 0.244675925467163,
0.617208154639229, 0.675795117858797, 0.921914914855734, 0.777150499867275,
0.841221024747938, 0.492060288786888, 0.326705329120159, 0.686669388087466,
0.121503615053371, 0.367948720930144, 0.106932319700718, 0.363004767335951,
0.350841429783031, 0.204824617831036, 0.233977065887302, 0.7365177560132,
0.0454314323142171, 0.865328938700259, 0.328468273859471, 0.187198075000197,
0.0820242392364889, 0.516697632148862, 0.481362534221262, 0.620970563730225,
0.817455855663866, 0.98652182915248, 0.980461026541889, 0.0709870441351086,
0.536431754939258, 0.976550402119756, 0.40865416219458, 0.890750366030261,
0.662430537166074, 0.388162778690457, 0.959246737649664, 0.931342099327594,
0.770709854317829, 0.279344801558182, 0.305849055293947, 0.370182184502482,
0.562077804701403, 0.240711282240227, 0.630282753147185, 0.170173112303019,
0.313635889207944, 0.694733579177409, 0.319224840262905, 0.34643409284763,
0.695703814737499, 0.39429243397899, 0.189475272549316, 0.516221221769229,
0.777665859088302, 0.558720602886751, 0.229818729683757, 0.484828221844509,
0.481090292567387, 0.0106754100415856, 0.684031964046881, 0.815850237151608,
0.484693359583616, 0.362359361024573, 0.116397734964266, 0.350968763930723,
0.150200203061104, 0.91834011208266, 0.52454294054769, 0.148975162534043,
0.510950653115287, 0.803136837668717, 0.918003916973248, 0.100692971842363,
0.0494276131503284, 0.526557822246104, 0.0758368712849915, 0.269769576610997,
0.101892914390191, 0.424805645598099, 0.110284168738872, 0.709118485450745,
0.784806538838893, 0.0519805552903563, 0.398219058057293, 0.34969964553602,
0.695627773646265, 0.549499853979796, 0.273849666584283, 0.151610748609528,
0.874867760110646, 0.198461376130581, 0.569974714424461, 0.725238451501355,
0.594061649870127, 0.287728592054918, 0.167462798533961, 0.126261120429263,
0.557413443457335, 0.464988547144458, 0.552185143576935, 0.0831529255956411,
0.333197562023997, 0.574352362193167, 0.4974719658494, 0.190328923054039,
0.82773098279722, 0.356611662777141, 0.803166383411735, 0.320444498909637,
0.366018867818639, 0.809054122306406, 0.553693874971941, 0.293376821791753,
0.97094266442582, 0.879881211323664, 0.308937143767253, 0.579827836947516,
0.496185180498287), V5 = c(0.923095905454829, 0.762332424288616,
0.972339176572859, 0.422411200823262, 0.940925979521126, 0.540574067272246,
0.985336700454354, 0.8358859855216, 0.452674505067989, 0.919172912370414,
0.177077513653785, 0.484894159948453, 0.228302743285894, 0.588133673649281,
0.120827580103651, 0.320581880165264, 0.132633685600013, 0.244240331463516,
0.0881568337790668, 0.0182364885695279, 0.619532398879528, 0.763255036668852,
0.360307943308726, 0.291531737195328, 0.637159892357886, 0.423571469495073,
0.144552590325475, 0.232900225091726, 0.431558183860034, 0.624517237534747,
0.660041945520788, 0.0605304769705981, 0.894808513810858, 0.766446133842692,
0.708146936027333, 0.980518834665418, 0.831859190016985, 0.63390138419345,
0.0906804192345589, 0.863758355844766, 0.81235816469416, 0.754735963186249,
0.96086330129765, 0.379362520063296, 0.86472422350198, 0.434156836709008,
0.921033434569836, 0.930532508296892, 0.919923858484253, 0.0185073432512581,
0.190673178527504, 0.389489938504994, 0.251335079083219, 0.398261951981112,
0.323421294335276, 0.0406879091169685, 0.0622786141466349, 0.511600321624428,
0.761615104507655, 0.267784109571949, 0.593046235619113, 0.923055295599625,
0.481184788281098, 0.539936339249834, 0.71660251962021, 0.525291009107605,
0.878714059945196, 0.127967000938952, 0.616779315285385, 0.174148306716233,
0.737354298587888, 0.931568023283035, 0.51300498447381, 0.860497320536524,
0.776048321742564, 0.0440801647491753, 0.0390385524369776, 0.765144015429541,
0.334515287308022, 0.312304191756994, 0.874887532787398, 0.495062640635297,
0.772741753375158, 0.847717159660533, 0.465354095678777, 0.988151066936553,
0.933549190638587, 0.493211443535984, 0.604139760835096, 0.11511691333726,
0.955369158880785, 0.00772590702399611, 0.396924669155851, 0.360593063989654,
0.44654577318579, 0.775453717214987, 0.103126123780385, 0.0639552725479007,
0.382581982528791, 0.221267726272345, 0.569829758256674, 0.189883300568908,
0.343288340838626, 0.372111623873934, 0.235726190032437, 0.540518779074773,
0.706296919612214, 0.939652142813429, 0.0611454977188259, 0.428340997546911,
0.625002122484148, 0.828457030234858, 0.038334577344358, 0.986632968764752,
0.838095314567909, 0.635932896286249, 0.375527458963916, 0.63829215336591,
0.54482639557682), V6 = c(0.370218734024093, 0.419946043053642,
0.545280206250027, 0.277894535567611, 0.00736652850173414, 0.790516802109778,
0.313465759390965, 0.809552269754931, 0.842788263224065, 0.499014808563516,
0.623795239022002, 0.548406617250293, 0.05658471561037, 0.965675859013572,
0.472378655336797, 0.524981135968119, 0.46187190245837, 0.387807991588488,
0.938760220305994, 0.50433289213106, 0.233549498720095, 0.00805983133614063,
0.127320852363482, 0.672702707583085, 0.415367575129494, 0.250937715405598,
0.0297437864355743, 0.406187633285299, 0.675050344318151, 0.292069549206644,
0.61900219717063, 0.616054858779535, 0.250656338408589, 0.518912112340331,
0.183584470301867, 0.0910171526484191, 0.546375629259273, 0.621274302946404,
0.0846014332491904, 0.621578456601128, 0.576011084485799, 0.504746064078063,
0.221732606412843, 0.69746913574636, 0.164020692463964, 0.158075062790886,
0.847464352613315, 0.605109791038558, 0.596407983452082, 0.0114535829052329,
0.063888885313645, 0.353724629618227, 0.612920761108398, 0.957113454118371,
0.741108618909493, 0.0766351209022105, 0.641778253018856, 0.988549496047199,
0.554500735830516, 0.490220351610333, 0.814872513525188, 0.473424580646679,
0.210200395900756, 0.70044357329607, 0.208990931045264, 0.680869682924822,
0.932328240247443, 0.718588735442609, 0.658000691328198, 0.185951039893553,
0.753725055372342, 0.123452978907153, 0.212786188349128, 0.228114538360387,
0.00365104852244258, 0.814572377130389, 0.124759285245091, 0.11509175482206,
0.3558327362407, 0.487612128257751, 0.344766818918288, 0.25334310438484,
0.712012055562809, 0.453843667870387, 0.938410549657419, 0.718660024693236,
0.325095203239471, 0.156200079247355, 0.175217281095684, 0.551201074384153,
0.830405764747411, 0.281535709509626, 0.403883191291243, 0.545596988638863,
0.282760329078883, 0.755758260609582, 0.310922948876396, 0.819937098305672,
0.0865648619364947, 0.936576428823173, 0.880421161418781, 0.454629295505583,
0.14402421656996, 0.840076276101172, 0.655401981435716, 0.569122347747907,
0.438042796216905, 0.0409106772858649, 0.944646059069782, 0.825017294846475,
0.0786613943055272, 0.0726817650720477, 0.588666438357905, 0.588839216856286,
0.444881158880889, 0.721039973665029, 0.534136849688366, 0.852300873491913,
0.281278134556487), V7 = c(0.141886763740331, 0.532325896434486,
0.588067500153556, 0.640411838656291, 0.946903814794496, 0.0120973803568631,
0.952789547387511, 0.107461815234274, 0.236731971381232, 0.79851196357049,
0.312241489300504, 0.458273995434865, 0.275089886970818, 0.966928232926875,
0.153488726355135, 0.250537436921149, 0.262579692527652, 0.622048743302003,
0.497019267408177, 0.081929708365351, 0.995605900418013, 0.385106679983437,
0.176982956239954, 0.696908152429387, 0.732609438244253, 0.187662536511198,
0.936794911511242, 0.168030542787164, 0.569703943328932, 0.604032611008734,
0.115415274398401, 0.807121691294014, 0.139127567177638, 0.530566544737667,
0.27088794275187, 0.702946513192728, 0.755043589742854, 0.811690449714661,
0.114138029748574, 0.233510547783226, 0.654041203204542, 0.441771358018741,
0.969852354610339, 0.964909290196374, 0.541466456372291, 0.0193931085523218,
0.201263489434496, 0.33053227327764, 0.687692743493244, 0.18348641670309,
0.418929822510108, 0.383909322321415, 0.137825107900426, 0.217210126807913,
0.453193115768954, 0.0340990910772234, 0.61968918191269, 0.780179363442585,
0.80423841974698, 0.491154700284824, 0.351554640568793, 0.587180010043085,
0.239602132700384, 0.644263347610831, 0.969654035987332, 0.547199918190017,
0.397612665314227, 0.722763631725684, 0.72788112401031, 0.795322017977014,
0.686627405928448, 0.762283272109926, 0.926369941327721, 0.894270734628662,
0.257203514222056, 0.0720642674714327, 0.985171274049208, 0.660412719240412,
0.648350277217105, 0.917640027124435, 0.694397774524987, 0.917451658984646,
0.524428785312921, 0.0992080934811383, 0.100496702594683, 0.975102295167744,
0.492288384120911, 0.0227630862500519, 0.813120178645477, 0.71656278334558,
0.913717438234016, 0.636773794423789, 0.929305002791807, 0.892870238050818,
0.780625773593783, 0.850657029775903, 0.204829484922811, 0.640264326473698,
0.663275289349258, 0.789437669329345, 0.485021750209853, 0.0749609414488077,
0.655508011812344, 0.74328700476326, 0.31854355102405, 0.0728525989688933,
0.382023223210126, 0.215404908405617, 0.656497932272032, 0.92760404618457,
0.781216817675158, 0.833905356237665, 0.361858946271241, 0.300276404013857,
0.373293370939791, 0.853469310561195, 0.133900812361389, 0.986530470428988,
0.411433706991374), V8 = c(0.849042785586789, 0.539557616924867,
0.564340554876253, 0.755421123001724, 0.984947801567614, 0.112984226550907,
0.403615512186661, 0.0420174030587077, 0.611705861985683, 0.659787082113326,
0.876130681717768, 0.948364004027098, 0.203974137781188, 0.309610587079078,
0.984034060733393, 0.534302587388083, 0.971052328590304, 0.812090678606182,
0.12758506112732, 0.4455144808162, 0.415974463103339, 0.518267513252795,
0.0146076162345707, 0.24067388032563, 0.0458081206306815, 0.366667992202565,
0.0854597128927708, 0.0380873722024262, 0.663318049395457, 0.873134173220024,
0.715449047507718, 0.83272837777622, 0.269359193509445, 0.936196858529001,
0.00233305594883859, 0.890992778819054, 0.980090971337631, 0.483554403064772,
0.535345176234841, 0.51444540428929, 0.802370720077306, 0.772680272115394,
0.102163553936407, 0.00971000618301332, 0.902426401618868, 0.880951991537586,
0.284606313100085, 0.0439106288831681, 0.937477112514898, 0.912608936196193,
0.998304838780314, 0.745025476673618, 0.751772192539647, 0.278778153471649,
0.162150844698772, 0.837488571181893, 0.388283311389387, 0.335549421375617,
0.650775108020753, 0.990717119770125, 0.0961552720982581, 0.302420736989006,
0.276211500866339, 0.733167931437492, 0.164790137903765, 0.249713592929766,
0.836560999508947, 0.692190924426541, 0.570303184678778, 0.984758368227631,
0.328286996111274, 0.917134278221056, 0.612915431382135, 0.242480249842629,
0.52148644789122, 0.660938010783866, 0.448608730221167, 0.691765318857506,
0.341251585632563, 0.304311832180247, 0.0227175708860159, 0.149960547685623,
0.17484575859271, 0.518439768115059, 0.639089794829488, 0.860548116965219,
0.705641856649891, 0.673352806130424, 0.100467606680468, 0.570258957333863,
0.706326137995347, 0.70817924500443, 0.453099914360791, 0.70276203728281,
0.252404118655249, 0.598284177947789, 0.741140816593543, 0.806638022884727,
0.720691111404449, 0.68312481045723, 0.973559587728232, 0.917065860237926,
0.874346838332713, 0.143847105791792, 0.453837911831215, 0.569988392293453,
0.692350469296798, 0.806599335046485, 0.456668916391209, 0.386287276633084,
0.908880397910252, 0.793711767531931, 0.754064193693921, 0.533496839459985,
0.615311589790508, 0.331158632645383, 0.11861584126018, 0.408010076498613,
0.222291024168953), V9 = c(0.125224817777053, 0.914628089172766,
0.377369399182498, 0.751476752338931, 0.0853602013085037, 0.846715328749269,
0.561333758756518, 0.0429282214026898, 0.284403225872666, 0.285037216497585,
0.604412623681128, 0.319880489027128, 0.91246378608048, 0.455940068000928,
0.867939765099436, 0.425015810178593, 0.470896222162992, 0.224302776856348,
0.36109847901389, 0.594786533853039, 0.642710687592626, 0.207564177224413,
0.284501956542954, 0.732748613227159, 0.774612213019282, 0.11284792679362,
0.956883376231417, 0.821359134977683, 0.846147159812972, 0.411850540433079,
0.770936917047948, 0.496579384431243, 0.619339638156816, 0.177370799472556,
0.239690004847944, 0.224350028671324, 0.0038589695468545, 0.179824534803629,
0.719484498025849, 0.798693082761019, 0.842053556581959, 0.294362392276525,
0.221064954530448, 0.898532920051366, 0.723221117863432, 0.863187251612544,
0.0994312795810401, 0.361342326970771, 0.783163121435791, 0.502090321853757,
0.13186371861957, 0.659234342398122, 0.390342378057539, 0.328385862288997,
0.446345729287714, 0.507106031524017, 0.152638432569802, 0.445379047188908,
0.304815052077174, 0.0760775990784168, 0.924889291869476, 0.0672700984869152,
0.199834800790995, 0.0825146653223783, 0.0361869866028428, 0.543002652470022,
0.972956772428006, 0.190790882101282, 0.376062947325408, 0.351586842909455,
0.689982761396095, 0.241374674253166, 0.838397472165525, 0.934896029997617,
0.468875553458929, 0.628862793324515, 0.431529769208282, 0.935133612947538,
0.352268195245415, 0.830344556132331, 0.666667894460261, 0.905744371935725,
0.7915226041805, 0.895240528741851, 0.584306629607454, 0.864065658533946,
0.0452204917091876, 0.55241105845198, 0.101811402710155, 0.749172917334363,
0.27551553142257, 0.83245463995263, 0.274472723482177, 0.257237296551466,
0.959463604493067, 0.109898063354194, 0.971408878453076, 0.0386630205903202,
0.975892154034227, 0.358955943025649, 0.399594583781436, 0.98816026840359,
0.964564403751865, 0.385808907682076, 0.42681196751073, 0.346436691703275,
0.723757655825466, 0.839694014983252, 0.847320703556761, 0.127265167189762,
0.192201391793787, 0.113802545471117, 0.373555006459355, 0.84189727040939,
0.515876488294452, 0.0691332728601992, 0.834895022911951, 0.19373700208962,
0.921317730797455)), .Names = c("startmonth", "V1", "V2", "V3",
"V4", "V5", "V6", "V7", "V8", "V9"), class = "data.frame", row.names = c(NA,
119L))
Here's what I get with a merge that is then restricted to only the results of the first two columns of each df:
> merge(Benret, Allret, by=1)[c(1,2,3,11,12)]
startmonth V1.x V2.x V1.y V2.y
1 Jan 2012 0.005769477 0.27738310 0.02842473 0.029462133
2 Feb 2012 0.205415812 0.05403291 0.00100873 0.003227936
3 Mar 2012 0.503058286 0.33427439 0.00221697 -0.012495625
4 Apr 2012 0.421452875 0.76272515 -0.01468005 -0.013674570
5 May 2012 0.171408428 0.62441319 0.00102046 0.006939153
6 Jun 2012 0.927357386 0.66497853 0.00591266 0.025957210
7 Jul 2012 0.736080843 0.87750471 0.01992673 0.062507000
8 Aug 2012 0.063494948 0.13381045 0.04174836 0.045244532
9 Sep 2012 0.628076423 0.30312729 0.00335606 0.006712116
10 Oct 2012 0.825837342 0.40190400 0.00334483 0.031480766
Calculate the differences and omit the inputs:
res <- merge(Benret, Allret, by=1)[c(1,2,3,11,12)]
> res$diffV1 <- with(res, V1.x - V1.y)
> res$diffV2 <- with(res, V2.x - V2.y)
> res[ -(2:5) ]
startmonth diffV1 diffV2
1 Jan 2012 -0.02265525 0.24792096
2 Feb 2012 0.20440708 0.05080497
3 Mar 2012 0.50084132 0.34677002
4 Apr 2012 0.43613292 0.77639972
5 May 2012 0.17038797 0.61747404
6 Jun 2012 0.92144473 0.63902132
7 Jul 2012 0.71615411 0.81499771
8 Aug 2012 0.02174659 0.08856592
9 Sep 2012 0.62472036 0.29641517
10 Oct 2012 0.82249251 0.37042323
So, I had to change the name of the rows to calculate the Anova. Now I'd like to bring back the previous names of the rows. Is there any way to do it ? Data:
head(tbl_reo)
id Sequence variable value pv
75 AAAAGAAAVANQGKK BiotinControl1 3893050.50 0.8523174
192 AAAAGAAAVANQGKK BiotinControl2 900604.61 0.8523174
3770 AAFTKLDQVWGSE BiotinControl3 90008.14 0.8523174
934 AAIKFIKFINPKINDGE BiotinTreatment1 656399.20 0.8523174
1012 AAIKFIKFINPKINDGE BiotinTreatment2 359937.85 0.8523174
1970 AAIKFIKFINPKINDGE BiotinTreatment3 447089.35 0.8523174
Code:
levels(tbl_reo$variable) <- c("Control","Control","Control","Treat","Treat","Treat")
New data:
head(tbl_reo)
id Sequence variable value pv
75 AAAAGAAAVANQGKK Control 893050.50 0.8523174
192 AAAAGAAAVANQGKK Control 900604.61 0.8523174
3770 AAFTKLDQVWGSE Control 90008.14 0.8523174
934 AAIKFIKFINPKINDGE Treat 656399.20 0.8523174
1012 AAIKFIKFINPKINDGE Treat 359937.85 0.8523174
1970 AAIKFIKFINPKINDGE Treat 447089.35 0.8523174
No idea if it's necessary to use dput on such big data, but here you go.
> dput(head(tbl_reo,4))
structure(list(id = c(75L, 192L, 3770L, 934L), Sequence = structure(c(1L,
1L, 3L, 5L), .Label = c("AAAAGAAAVANQGKK", "AAANILFSDHAD", "AAFTKLDQVWGSE",
"AAIELRE", "AAIKFIKFINPKINDGE", "AAIYKLLKSHFRNE", "AAKKFEE",
"AAKYFRE", "AANVKKTLVE", "AARAGELLKE", "AARDSKSPIILQTSNGGAAYFAGKGISNE",
"AAVKARVASIDE", "AAVKASAPGSVILLE", "AEKLKAE", "AEQVKKE", "AFAKRQGKE",
"AFGSGTAAVVSPIKE", "AGIPAGVVNILPGSGRVVGE", "AGISKDGQTRE", "AGSHGKDAWGVFRVHAFE",
"AHSFSSSLASAE", "AIDALNGMLLNGQE", "AITYSWTLLTE", "AKAKKAAKKAANAKQQQE",
"AKALVAQGVKFIAE", "AKAQGVAVQLKRQPAQPRE", "AKIGKKVE", "AKKVVKKAFEE",
"AKNKKYALDIIKE", "ALANAGVTDIVLAVNYRPE", "ALDFFARGLVKSPIKVVGLSTLPE",
"ALFHPSVLGLE", "ALKIWKRFRE", "ALKPDSQKSYAE", "ALLLKVNQIGTLSE",
"ALRIGSE", "ALRLVE", "AMRIGSE", "ANIIKLVE", "ANNKWNKLFSE", "ANVTGLRDAMFKGE",
"AQKVRDQITSLKE", "AQRVQAKNQLE", "ARINFGIEKAKE", "ASAGIQLSPKFVKLVSWYDNE",
"ASFVLPTWAAKE", "ASKLAAYNKKSTISARE", "ATSSIVMFSLNE", "ATTNAPWYKGWEKE",
"ATVGAGLPIISFLRE", "AVDATHLAE", "AVIRLSKE", "AVSKVYARSVYDSRGNPTVE",
"AVVTVPAYFNDSQRQATKD", "AVWYGPPKAANLGGVAVSGLE", "AWFPYQKE", "AYSPFGSANAPLLKE",
"DGWVPGKNLFTVE", "DHSFINYKQNVE", "DILYKIYKGVSE", "DLYQTFARTVE",
"DNVKAHPRIE", "DQVRNYVMQYIQE", "DSVYLAKLAE", "DVHNLFTYAKE", "DVKISAGTLLKGTE",
"DVQKFRHE", "EAVAKAAE", "EKPLFLVFHGGSGSTVQE", "EQRPGTPLFTVKAYLPVNE",
"ERIAALE", "FAKSLPRLDGLSVQE", "FALKHFPDYE", "FAPYGTITSAKVMRTE",
"FASLFPLKQAPAFLGPKGLKLTE", "FEQKNAGGVGE", "FFIDRPE", "FFIRDTKNLE",
"FGINTTVQE", "FGNRPMILHQIE", "FGVTKARWPAKSVKFVQGLLQNAAANAE",
"FIFSKYQILE", "FKMAKPE", "FKSMIKNNE", "FKTKYNSDPYQLE", "FLAADKE",
"FLKKILVDNKLE", "FLLSYGKE", "FMIAPTGAKTFAE", "FMISQLGLQKKNIKIHGF",
"FNRGLALE", "FNSQPDYRE", "FQKDAKRFE", "FQKTGPFE", "FRNGMINVSPIGRNASTEE",
"FRTASAQLE", "FSGSLLRSKFNTSNE", "FSRWFNHIASKADE", "FSTNVKDE",
"FTADGKVKPE", "FTPETPSPLIGILE", "FVADGLLQGFGNE", "FVAYPIQLLVTKE",
"FVAYPIQLVVTKE", "FVDIINSRN", "FVSSKHGSVSAE", "FYKSLKLPYRIVGIVSGE",
"GAAPGGFPGGAPPAPE", "GAIWNINSFDQWGVE", "GARYGPSLMPGGSEE", "GEKTMGAKITKSLE",
"GGLGPINIPLLADTNHSLSRD", "GGTLIGTARSME", "GGVAPNIQTAEE", "GGVIIITHSAE",
"GIFSGTLSYIFNE", "GIKVPIHIVQE", "GILQRWTKGFDIPNIE", "GILQRWTKGFDIPNVE",
"GILVASPFVELE", "GIPAGWQGLDNGPE", "GIPPDQQRLIFAGKQLE", "GISKKTRIAGVVYHPSNNE",
"GKKIATYQE", "GKLKGVLGYTE", "GLDKIFQE", "GLGWGNNASAAIQRVGLGE",
"GLNHLPVQKQIVE", "GMMANGIAAGGQLTTTTE", "GMRWAGNANE", "GNDIHALAAKLLQE",
"GQALKKKKGKSAE", "GRGINSSYAVWNDVTGKAHE", "GRHSGYGKRKGTRE", "GSDVLNTWKKE",
"GSRKVDGQKVKASKE", "GTAVSQADVTVFKAFQSAYPE", "GTERSAYYTWAAE",
"GTRAVTKYSSSTQA", "GVVLGVEKRATSPLLE", "GWGTPHIKPYGNLSLDPSACVFHYAFE",
"GYSGVEVRVTPTKTE", "GYVGRTQPKKGGEGE", "GYWNGKITVDE", "HDHAGVPLKISPPVVAYRE",
"HIIQAKADAQRE", "HIVPGTYIQE", "HKFAIPAINVTSSSTAVAALE", "HKVVNGNSADAAYAQVE",
"HLQPHFTTNDGLE", "HNLRVISE", "HQIRAKFVE", "HQKYTREQVGCKE", "IAKLNIPTGIPLVFE",
"IAQHIFTKLNSNGYLEE", "IDASKPVKLTLLISE", "IDPINALQAAME", "IDVAWLKANAKE",
"IEAISSPLMGIFE", "IFDANVQRLKE", "IFFKIKKTTPLRRLME", "IFGGYLYFAQAPSAAE",
"IFGGYLYFAQAPSAAEFHTE", "IFGPVVTVAKFKTLEE", "IFGPVVTVTKFKSADE",
"IFLHSLPVKE", "IFSTYADNQPGVLIQVFE", "IGATTSVFPFNKSMIEYLE", "IGNSGMFRPEMLE",
"IGSLTRQGAKGGLLE", "IHPLLPQDE", "IIAGNVATRE", "IIGGTTRIPTLKQSISE",
"IIKKNFDLRPGVLVKE", "IISNAKGQVVGSIIALDRQE", "IKGFFSFATQKLDE",
"IKHLPFNVVNKDGKPAVE", "IKKVVKAAAE", "IKLFNKWSFEE", "IKLVKLIQE",
"ILIVTKDGVRKE", "ILKIVPGRVSTE", "ILKKLDE", "ILLARNLI", "ILNKLALRILKNE",
"ILNSQAPEE", "ILSKLASAGATE", "ILWRQKE", "IMLPVFDAPQNLVE", "IMLQLGIPREE",
"IMLWNLAAKKAMYTLSAQDE", "INLPFITADASGPKHINMKFSRAQFE", "INLWFKKEE",
"IQLNGGSISE", "IQPRANIQLDFPE", "IRVFQGE", "ISAMVLTKMKE", "ISDTMWPGQAMTLKVE",
"ISDVFAKWNNGVLD", "ISLAAFGRKE", "ISTMIRWGLKPYLFVLNNDGYTIE", "ITDIKQFLE",
"ITHPTIVDGWFRE", "ITLGKYLFE", "ITQLMSLIINTVYSNKE", "ITVDAPKAHE",
"ITYLNNRAAAEYE", "IVERLNTASIPE", "IVFASPPRKSDGKVSPE", "IVFNMSREE",
"IVNLTLPDGTVRQGQVLE", "IVPKLME", "IVRVNKNAPAVE", "IVSRFKTE",
"IYEKME", "IYFSLNNGVRIPALGLGTANPHEKLAE", "IYNHIQLRE", "IYNHIQLREE",
"IYPYGSRGPKGLKE", "KAAKVTKAAQKAAKK", "KAALQTYLPKIKE", "KAEQVFQNVKNILAE",
"KAITYSHAAKE", "KAKIGQKE", "KAKKANNPQHSITKDE", "KALTYTPPKKQKAEKPKAE",
"KAQRAPLEE", "KAVAYWKTLKTDE", "KDFATFAKE", "KDYASVVPPRHNE", "KEKLTQQKKE",
"KFAAGTKALLDE", "KFINSEIE", "KFNAVPIFLSDE", "KFSERLLNE", "KGAQTYVKKVIGSFKDWE",
"KGCDIIVKNHE", "KGLLSATQFTQPALTLME", "KGLSKKGGNIPMIPGWVMDFPTGKE",
"KGQIVGRYVVDTSK", "KGTGKSNKITITNDKGRLSKE", "KGVFRSIVPSGASTGVHE",
"KGVFRSIVPSGASTGVHEALE", "KGVITTQAE", "KGWVPVTKLGRLVKAGKITTIEE",
"KIFTVPTETLQAVTKHFISE", "KIGDLNFLIGPKLYE", "KKAADALLLKVNQIGTLSE",
"KKAGKTLDVVE", "KKAPAAAPAASKKKEE", "KKIDLASIRE", "KKPKTKKVKEE",
"KKQAIRSKQE", "KKTISLKTPLKTLVILHSDE", "KKVYPDVLYTSKLSRAIQTANIALE",
"KKYIATIGVE", "KLAAQRKAE", "KLIHGPKAQYNE", "KLIHNDSSYNHE", "KLKAAYSYMFDSLRE",
"KLKAERE", "KLKKVLSANTNAPFSVE", "KLSGIANNRDE", "KLVAVLK", "KLYVGRAQKKNE",
"KMIKPKEGRIE", "KMSSGYYLGE", "KMSTRAAPFEQRLPE", "KNLFTGWVDVKLSAKGQQE",
"KNLKPLAIIKGWGE", "KNMDAATALVGSGPAFVLLMLE", "KNVPLYQHLAD", "KQLYSFDLE",
"KQNLLAVE", "KRVSGMFLGE", "KSPFLDALKAKNFE", "KTFVNVAKPFHKE",
"KTKLQGMLNKAEE", "KTMGAKITKSLE", "KTQAFKKHVDE", "KTVRKHLERQGWIE",
"KVFNAYPAARE", "KVGLRLSPYGVFNSMSGGAE", "KVILKMKE", "KVIRFQEE",
"KVSVLKALGAE", "KYFPGLLRATNE", "KYSLAPVAKE", "LAKKQKKLE", "LDALIAQGRE",
"LDASNKAAKE", "LDLARPIYLPTASYGHFTNQE", "LDTAQKHID", "LDTIRNME",
"LEARIPE", "LFAKFGPIVSASLE", "LFAKHGE", "LFDTHRVIVNE", "LFLKQLVVGGLDRVYE",
"LFRSIGGE", "LGAVALKGALAKVPE", "LGAVNLHE", "LGDKAVYAGE", "LGDNAVFAGE",
"LGGEKKQKGQALQF", "LGGKSPNIVFADAE", "LGGTVVSLSDSKGCIISE", "LGILPRDILTKE",
"LGIYSHE", "LGKGSFKYAWVLD", "LGKVLAKVIGKE", "LGMLAGADRVE", "LGTPFGVTIDFE",
"LHGNAKKAAEE", "LICSYRSKIETE", "LIGGSADLTPSNLTRWKE", "LIHQNEVPLVLLSSGVGVTPLLAMLE",
"LIHSKVGGE", "LIKLTGKLIE", "LILNGGFSPLTGFLNE", "LINNLGTIAKSGTKAFME",
"LKADLRPLQIKSIRE", "LKAHNINVVDK", "LKAQSKPHVGDE", "LKHLVGPEKAAE",
"LKIQRPLHE", "LKKKGILFVGSGVSGGEE", "LKKQKTAE", "LKKQQKEAE", "LKKSANVGKKKE",
"LKPTVILKLLKE", "LKRAIDLNKE", "LKRIMINCFNE", "LKRIPE", "LKTNGKFQVDKSIIE",
"LKYFGKALE", "LLAAFKAYLE", "LLAVPVVKGRKTEKE", "LLAYQFASPVRWIE",
"LLESPNNKVVPATNQIE", "LLKTNRISNE", "LLKTNRISNEQYE", "LLLLDVAPLSLGIE",
"LLLQQRATE", "LLRLVLLE", "LMAKLANDSPLAIE", "LNAATGSKTSLYDKKE",
"LNAAYAADGYARIKGMSCIITTFGVGE", "LNKGLTKKGGNIPMIPGWVMEFPTGKE",
"LNNAAAKKYDLE", "LNQLLKLVSE", "LPAWSKLQKIYE", "LPGVAFLSEKK",
"LPKNSTSSKDKKKKN", "LQAQLDYLNAE", "LQDIANPIMSKLYQAGGAPGGAAGGAPGGFPGGAPPAPE",
"LQGKLTGMAFRVPTVD", "LQGKLTGMAFRVPTVDVSVVD", "LRPLQFKSIRE", "LRQATGGQAFPQMVFD",
"LSALNGIAGSYAE", "LSFALIHKE", "LSFPFEWAFAGTNE", "LSGIPPAPRGVPQIE",
"LSKISGVDSSLIFPALE", "LSKSKTSPYVLPVPFLNVLNGGSHAGGALALQE", "LSLLFNVIPDRYTKYD",
"LSNVINMFFE", "LSQGTWLNKPKSVFQE", "LSRIGEFE", "LTQFPAFVTPMGKGSIDE",
"LTQLKKE", "LTSLTDYVTRMPE", "LTVKLNKE", "LVARVQGGEE", "LVFGQTFTDHMLTIPWSAKE",
"LVGRINGQFGTVE", "LVILGDSNVVRE", "LVKVGHDNLVGE", "LVRDNKLIGNFTLAGIPPAPKGVPQIE",
"LVRFSYGQDLE", "LVSLPNIILTPHIGGSTEE", "LVTAPLDGTILE", "LYGNIVMSGGTTMFPGIAE",
"LYGTAIGARSQGAKTYLE", "LYPGRPLAIALD", "LYPGRPLAIALDTKGPE", "LYSTSSSSTSSSATKN",
"MATLYDGLE", "MAVAIKKE", "MAVKFGKNAF", "MFFVSKVMNNE", "MFNDFLNAGAE",
"MGQKKAPDGRIVE", "MIGSVVGIYNGKAFNQVE", "MIGSVVGIYNGKAFNQVEIRPE",
"MIISSPSDGQVKE", "MIVLHLPSPVTAQAYRAE", "MKPTSIEKE", "MLAQDKLVAE",
"MLGHYLGE", "MLGNWSFGDYFKKE", "MLNRISPKILDE", "MLRAAGKE", "MNNIIAASRRH",
"MPPQIPNE", "MQTAAQSSSIE", "MQTEIKE", "NAHLKKNFKPQGSIE", "NALKLAKWLE",
"NAVDRLLVE", "NFAPFLKTILPE", "NFIQSIME", "NGFINNPIVISPTTTVGE",
"NGHSPSEAFNE", "NGITFPKDAPSPFIFKTLE", "NGIVTYAAKE", "NGLLSFVQYVIAPIQE",
"NGLTAYRLGKE", "NGSFLLYIDRHLVHE", "NGTVTAANASPINDGAAAVILVSE",
"NGVAKGSFKD", "NISNNYMIGAINAE", "NKERILTE", "NKIIRMPLVE", "NKIVRKPLME",
"NKKVIITGAPAAFSPTCTVSHIPGYINYLDE", "NKVQYIAGARPWTHVQKVDIALPCATQNE",
"NKVSPADAAKKAL", "NLHILE", "NLKDTSSGGVTHANE", "NLKKGSTVTSNGLNWE",
"NLKPSKPSYYLDPE", "NLLPLQSGIGNIANAVIE", "NLTVNPNVNE", "NMKLKLIEPLRE",
"NNIFNMVVE", "NNIFNMVVEIPRWTNAKLE", "NNMLHGNTMTVTGDTLAE", "NNYKKRLLGTAFKSSPFGGSSHAKGIVLE",
"NQGLKTATYPYWANE", "NQQGVNNFDE", "NRGIKVANPVWSTIPFISE", "NRKISAVSTYFE",
"NRNVGKTLVNKSTGLKNASDALKGRVVE", "NRNVGKTLVNKSTGLKSASDALKGRVVE",
"NRVKSDMLKE", "NSLLSSGTSQDSLRE", "NTQLMNRDNIE", "NTSLSPKFVYQGE",
"NVKTDSRNMPVKE", "NVQLYNGFVFHTGSLEE", "NVQNGEVRFE", "NVRDKGNSALLE",
"NVSSNIVKNAE", "NYNLLDTSGVAKVIE", "NYTVERVNE", "PKLVLVRHGQSE",
"PKLVLVRHGQSEWNE", "QAKLTAATNAKQ", "QAVGGIYSVLNKKRGQVVSE", "QDKKKKSNH",
"QDPSKSKVVVFE", "QDVPLIVPVVNPE", "QGAQVLFASTDSE", "QGKPYISLPKGKGIKLSIAEE",
"QGLKIDGPGFE", "QHPRYGGVYVGTLSKPE", "QHPRYGGVYVGTLSKPEVKE", "QISSMVLGKMKE",
"QIVYALKLPGIIHIDAAE", "QLADYLIGVQY", "QLASQLVDLKKE", "QLGIHYE",
"QLYKNDSNNNNNNNGNNAE", "QQGKNATVSVEE", "QRPGTPLFTVKAYLPVNE",
"QRTKVEHVRE", "QSPRPGQQAFE", "QSPYVSWVSYPGLASHSHHE", "QSTKFGTEIITE",
"QTYKKPTGGIVTVRSE", "QVADISSAKVNFKSLE", "QVAKWVNLAQKE", "RAVPHGKWQDE",
"RDLIKKKTKNNE", "RFLNVFPKLVEE", "RGEIKKGE", "RGLPVKLTYTDNKTSKE",
"RHIQVGDGLE", "RLAIVGVE", "RLLWYLARE", "RLNTASIPE", "RLNTASIPENVE",
"RLTSLNVVAGSDLRRTSIIGTIGPKTNNPE", "RLVAIAEQD", "RMHVLKKQYE",
"RNAAYQKE", "RNAAYQKEYE", "RNKIHQYLFQE", "RQLNAVLE", "RQLNAVLEE",
"RTSLFLNLANDPTIE", "RTVLFPIKYHE", "RVKANAAKRAE", "RYKLVPGIFADVKNLKE",
"RYKYVDPNVLPE", "RYVFLLDPMLATGGSAIMATE", "SALSYAALILADSE", "SCKMPQTVE",
"SDKSKWLTGVE", "SEYWPRLTKE", "SFAANWGVMVSHRSGETE", "SFQAVYNKLTGKQIVFE",
"SGSKLKFYVLHGDGVE", "SGTLGTKGNTQVIIPRLTE", "SGTVFDSSYSRGSPIAFE",
"SHIGVAKKLE", "SHLIPSATTGE", "SIAVLSAMKME", "SIAYSLKNTISE", "SIQKTKE",
"SIRVGEYNAE", "SISDSLSSKPHNFE", "SKFAKAYAQGVSKKE", "SKFKGDITWLPVRRKAYWE",
"SKLPTYTAKDSAVATRKLSE", "SLDQAWSLLRIYPKE", "SLGSPSGATKARIVVAKSGEL",
"SLIDAKTGLPKE", "SLKTGNAGPRPACGVIGLTN", "SLQNKGQEVKE", "SLRGLGGILLNPITGRRFVNE",
"SLSSKLSVQD", "SNLKKITLE", "SNPTAKLNLIWPATPIHIKKYE", "SQGKTLSVKQE",
"SQLTNLFNEE", "SQRIASKNQLE", "SRFQEIVKE", "SRKLFAATVAKAKTIVWNGPPGVFE",
"SSANRKIVKE", "STKNGTETKRSLE", "STPDTPSPLIAVNE", "STVAGFLVGSE",
"STVLQFKE", "SVKLTKAVAE", "SVLAHL", "SYQKSQKSQKE", "SYVDKHPVVTFNQE",
"TAGLTVGDPVLRTGKPLSVE", "TAIGPSGQLKGYLRPE", "TAKGNYPINAVTTMAE",
"TAQGQFLNFNKLLE", "TARSTATGPSEAVWYGPPKAANLGGVAVSGLE", "TDPFKLSGAQVVD",
"TEPLGTAGPLKLAE", "TFKLFNKYRPE", "TFRTTNTL", "TGGPISVPVGRE",
"TGIVAQYAYVLGELE", "TGVIKPGMVVTFAPAGVTTE", "THINLKVSDGSSE", "TIGNPKYNVPDFE",
"TIQDKSSAGAVVVANAKSLE", "TKAGVVKGKTLLE", "TKKSCCSGK", "TKNFTPEQISSMVLGKMKE",
"TKTFSPQE", "TKYLSGIAKRLNKE", "TLKKFGEE", "TLKQLNASLADKSYIE",
"TLKVVDPE", "TLRKVVKHFIDE", "TMKAVVIE", "TMWDTKKE", "TNPGTDVTVSSVE",
"TQALLLAPIAPHFAE", "TQKGVIFYE", "TRLSLTHMVE", "TSGSTVNDPLANYE",
"TSNFIKKVGYNPKTVPFVPISGWNGDNMIE", "TSSKKPKHPLE", "TTAMITDIATAPAE",
"TTHIKFASDPGCAFTKSIGFE", "TTTGVHHLYRMVKE", "TVAASAVAAVFE", "TVAKSKNAE",
"TVFAGQKP", "TVLRQALGE", "TVQVNLPVSLE", "TVSKMAYLDKTGE", "VAAIMQDPVMQSILQQAQQNPAALQE",
"VAALAAENK", "VAKVVGANPAAIKQAIAANA", "VANPIMSKLYQAGGAPE", "VARVILTQVGSGPQETNE",
"VATRAKQGELLE", "VDASGFRIPQQE", "VDSVLKHMKE", "VEGMRWAGNANE",
"VFIDFTKE", "VFKFFGFTPE", "VFKFFGFTPEGVAE", "VFLNEHPE", "VFVATKE",
"VGAKGVQLLSSYITEE", "VGASMTRSIQTLE", "VGKSVAVDSSE", "VGLKRVVTKAMSSR",
"VGTMNVFFVFLNKVTGKKE", "VGVHVLPKILDE", "VHFVSNIDGTHIAE", "VIAGNVVTRE",
"VIDTILALVKD", "VKILNGFAFVE", "VKVVLFNHSQRD", "VLANYLTSALSE",
"VLAVQKKLIAKSNLAGKPVICATQMLE", "VLGPFAFRIGSVKE", "VLIKRGVKPE",
"VLLLDVTPLSLGIE", "VLLVAVGRRPYIAGLGAE", "VLNPQIIKDVLE", "VLPRAVGSLTFDE",
"VLVLRGTLE", "VLVVRGSKKGQE", "VLYTSKLSRAIQTANIALE", "VMSIGRTFEE",
"VMTKNPVTGIKGITLKE", "VQNSTLAQLTSKLIPE", "VQSAVLGFPRIGPNRE",
"VRVTPTKTE", "VSAAYKNSLE", "VSDVAQKAE", "VSGLQYIDLLLIHSPLE",
"VSNPNIIFFRGTYAAVSPE", "VSNPPAYGAKIVAKLLETPE", "VSNRSTPSVVGFGPKNRYLGE",
"VSSMPTLIFYKGGKE", "VSVVDLTVKLNKE", "VTAALRVTD", "VTPLVQNVTGDKE",
"VTRVVGANPAAIKQAIASNV", "VTSSFFYRGNYTDFE", "VVALNDPFITND", "VVKILRE",
"VVLPVDFIIAD", "VVLTHPADE", "VVNIGIGGSDLGPVMVTE", "VVTSASAPPENASE",
"VVVITKNVSGQDVAAALE", "VYHNLKSLTKKRYGASAGNVGDE", "VYKIGGIGTVPVGRVE",
"WAGLGKIDRE", "WFDIKAPSTFE", "WFRIYKIPDGKPE", "WGKFLTCVEVE",
"WLFKLREE", "WMWNQLGNFLQE", "WSKLISE", "WSLTAPLGKALHDLS", "WTKIFKAIHE",
"YAQLLAKRLSE", "YAQLLAKRLSERKAE", "YFATHITE", "YFATHITERE", "YFQLLNSALTE",
"YGYSARVVDLIE", "YGYSTRVVD", "YGYSTRVVDLVE", "YIALVDQLIAKYSS",
"YIASKKGSISAE", "YKDIPVPKPKANE", "YKPLFTE", "YLLKAIKE", "YLNLPEHIVPGTYIQE",
"YNKIGDILSGRLKLRAE", "YQSKIIVFPRDGKAPE", "YQTQVLKNAKALE", "YQYLRSVVNE",
"YRKVLPIVSVPE", "YRLSVLPDNVPIMSVE", "YSLLAWTNIPRKE", "YVAKVSARE",
"YVNKGQDFHIAGE", "YVPSVIEPSFGIGRIIYSVFE", "YVQAAKNSIAAGADGVE"
), class = "factor"), variable = structure(c(1L, 1L, 1L, 1L), .Label = c("Control",
"Treat"), class = "factor"), value = c(3893050.501, 900604.6087,
90008.14298, 656399.2021), pv = c(0.852317422616546, 0.852317422616546,
0.852317422616546, 0.852317422616546)), .Names = c("id", "Sequence",
"variable", "value", "pv"), row.names = c(NA, 4L), class = "data.frame")
Just save the original variable before changing it.
var.save <- tbl_reo$variable
levels(tbl_reo$variable) <- c("Control","Control","Control","Treat","Treat","Treat")
# do whatever you need to do with changed variables
# ...
# ...
tbl_reo$variable <- var.save # change back
Having an object of class zoo we can select observations for a range of dates of interest using the function:
window(z, start = as.Date("2006-01-05"), end = as.Date("2006-01-08"))
In running this function the following warning message occurs:
Warning messages:
1: In which(in.index & all.indexes >= start & all.indexes <= end) :
Metodi incompatibili ("Ops.POSIXt", "Ops.Date") per ">="
2: In which(in.index & all.indexes >= start & all.indexes <= end) :
Metodi incompatibili ("Ops.POSIXt", "Ops.Date") per "<="
I have checked that the object is of class zoo and that Dates are included in the time series.
How is that possible?
Below the str(z) as requested:
‘zoo’ series from 2006-01-03 to 2013-01-24
Data: num [1:1795, 1:40] 3.65 3.68 3.69 3.72 3.7 ...
- attr(*, "dimnames")=List of 2
..$ : chr [1:1795] "1" "2" "3" "4" ...
..$ : chr [1:40] "EURARS" "EURAUD" "EURBRO" "EURCAD" ...
Index: POSIXct[1:1795], format: "2006-01-03" "2006-01-04" "2006-01-05" "2006-01-06" ...
Below the dput(head(z)) as requested:
structure(c(3.6511, 3.6833, 3.6931, 3.7152, 3.7027, 3.6897, 1.62349,
1.62257, 1.62011, 1.6115, 1.60243, 1.61108, 2.802, 2.7692, 2.7727,
2.7741, 2.7238, 2.729, 1.38937, 1.39109, 1.40716, 1.41627, 1.41196,
1.40666, 1.55055, 1.5472, 1.5448, 1.54335, 1.54215, 1.545, 623.73,
624.16, 628.43, 638.11, 632.27, 630.7, 9.6988, 9.7803, 9.7689,
9.802, 9.7492, 9.7354, 2742.03, 2765.68, 2758.65, 2769.27, 2753.3,
2747.31, 29.047, 28.972, 28.9, 28.88, 28.764, 28.792, 7.4616,
7.4601, 7.4612, 7.458, 7.46, 7.4589, 6.8983, 6.9551, 6.9594,
6.9838, 6.9374, 6.9253, 0.6882, 0.68905, 0.68961, 0.6863, 0.68473,
0.68358, 9.3178, 9.3963, 9.3889, 9.4207, 9.3702, 9.3516, 251.72,
250.28, 250.66, 250.39, 249.89, 250.86, 11657.46, 11677.26, 11612.05,
11603.59, 11433.35, 11403.84, 5.5244, 5.5808, 5.5799, 5.6134,
5.5858, 5.5957, 53.5288, 54.0151, 54.0323, 54.0105, 53.5591,
53.6189, 74.96, 74.88, 74.41, 73.94, 73.79, 73.84, 139.6, 140.71,
140.39, 139.09, 138.39, 137.93, 1208.3493, 1210.0214, 1195.7746,
1200.3966, 1181.3457, 1184.9635, 160.65, 162.15, 162.02, 162.53,
161.7, 161.43, 10.9802, 10.997, 10.9947, 10.9635, 10.9909, 10.9874,
12.7724, 12.8255, 12.8746, 12.8338, 12.7859, 12.8273, 4.5416,
4.5702, 4.5623, 4.5597, 4.5308, 4.5229, 7.9654, 7.9248, 7.9254,
7.914, 7.9574, 8.0011, 1.7571, 1.7626, 1.7622, 1.7574, 1.7411,
1.7391, 4.1276, 4.1608, 4.1665, 4.1818, 4.1606, 4.1534, 63.2627,
63.4733, 63.6725, 63.7198, 63.3608, 63.3412, 3.8295, 3.8116,
3.807, 3.805, 3.7609, 3.7799, 3.6732, 3.6815, 3.6834, 3.6842,
3.6644, 3.6492, 34.5475, 34.8363, 34.7254, 34.8369, 34.68, 34.39,
9.3648, 9.3279, 9.3321, 9.3152, 9.3389, 9.3603, 1.9844, 1.9932,
1.9938, 1.9902, 1.9766, 1.9716, 48.9853, 48.9426, 48.6762, 48.3184,
47.9995, 48.0187, 1.6149, 1.6195, 1.6193, 1.6204, 1.6155, 1.6148,
1.6129, 1.6175, 1.6184, 1.6201, 1.6182, 1.6221, 39.2261, 39.1868,
38.7569, 39.1189, 38.6148, 38.6309, 6.0673, 6.114, 6.1208, 6.1484,
6.1095, 6.1082, 1.2019, 1.2119, 1.211, 1.2151, 1.2088, 1.2065,
7.4834, 7.4559, 7.4658, 7.3872, 7.3206, 7.3497), .Dim = c(6L,
40L), .Dimnames = list(c("1", "2", "3", "4", "5", "6"), c("EURARS",
"EURAUD", "EURBRO", "EURCAD", "EURCHF", "EURCLP", "EURCNO", "EURCOP",
"EURCZK", "EURDKK", "EUREGP", "EURGBP", "EURHKD", "EURHUF", "EURIDO",
"EURILS", "EURINO", "EURISK", "EURJPY", "EURKRO", "EURKZT", "EURMAD",
"EURMXN", "EURMYO", "EURNOK", "EURNZD", "EURPEN", "EURPHO", "EURPLN",
"EURRON", "EURRUB", "EURSEK", "EURSGO", "EURTHO", "EURTND", "EURTRY",
"EURTWO", "EURUAH", "EURUSD", "EURZAR")), index = structure(c(1136242800,
1136329200, 1136415600, 1136502000, 1136761200, 1136847600), class = c("POSIXct",
"POSIXt"), tzone = ""), class = "zoo")
You shouldn't be comparing Date with POSIXct
Try this:
window(z, start = as.POSIXct("2006-01-05"), end = as.POSIXct("2006-01-08"))
Alternatively, as #JoshuaUlrich points out in a comment, your data is daily frequency, so you'd be better off using a Date index class for your data to avoid timezone weirdness.
index(z) <- as.Date(index(z))
window(z, start = as.Date("2006-01-05"), end = as.Date("2006-01-08"))
It looks like window doesn't convert between Date and POSIXct automatically. You have to specify the times in the same class as in the data:
window(z, start = as.POSIXct("2006-01-05"), end = as.POSIXct("2006-01-08"))
EURARS EURAUD EURBRO EURCAD EURCHF EURCLP EURCNO EURCOP
2006-01-05 23:00:00 3.7152 1.6115 2.7741 1.41627 1.54335 638.11 9.802 2769.27
EURCZK EURDKK EUREGP EURGBP EURHKD EURHUF EURIDO EURILS
2006-01-05 23:00:00 28.88 7.458 6.9838 0.6863 9.4207 250.39 11603.59 5.6134
EURINO EURISK EURJPY EURKRO EURKZT EURMAD EURMXN
2006-01-05 23:00:00 54.0105 73.94 139.09 1200.397 162.53 10.9635 12.8338
EURMYO EURNOK EURNZD EURPEN EURPHO EURPLN EURRON EURRUB
2006-01-05 23:00:00 4.5597 7.914 1.7574 4.1818 63.7198 3.805 3.6842 34.8369
EURSEK EURSGO EURTHO EURTND EURTRY EURTWO EURUAH EURUSD
2006-01-05 23:00:00 9.3152 1.9902 48.3184 1.6204 1.6201 39.1189 6.1484 1.2151
EURZAR
2006-01-05 23:00:00 7.3872