Я знаю, что у него есть другие вопросы с этой ошибкой, но, поскольку я новичок в R, они не помогли мне.
Я работаю в этих data.frame:
structure(list(Dates = structure(c(16071, 16071, 16071, 16071,
16071, 16071, 16071, 16071, 16071, 16071, 16071, 16071, 16071,
16071, 16071, 16071, 16071, 16071, 16071, 16071, 16071, 16071,
16071, 16071, 16071, 16071, 16071, 16071, 16071, 16071, 16071,
16071, 16071, 16071, 16071, 16071, 16071, 16071, 16071, 16071,
16071, 16071, 16071, 16071, 16071, 16071, 16071, 16071, 16071,
16071, 16071, 16071, 16071, 16071, 16071, 16071, 16071, 16071,
16071, 16071, 16071, 16071, 16071, 16071, 16071, 16071, 16071,
16071, 16071, 16071, 16071, 16071, 16071, 16071, 16071, 16071,
16071, 16071, 16071, 16071, 16071, 16071, 16071, 16071, 16071,
16071, 16071, 16071, 16071, 16071, 16072, 16072, 16072, 16072,
16072, 16072, 16072, 16072, 16072, 16072), class = "Date"), Identification = structure(c(1L,
2L, 3L, 4L, 5L, 6L, 8L, 9L, 10L, 11L, 12L, 13L, 14L, 15L, 16L,
17L, 18L, 19L, 20L, 21L, 22L, 23L, 24L, 26L, 27L, 28L, 29L, 30L,
31L, 32L, 33L, 34L, 35L, 36L, 37L, 38L, 39L, 40L, 41L, 42L, 43L,
44L, 45L, 46L, 47L, 48L, 49L, 50L, 51L, 53L, 54L, 55L, 56L, 57L,
58L, 59L, 60L, 61L, 62L, 63L, 65L, 67L, 68L, 69L, 71L, 72L, 73L,
74L, 75L, 76L, 77L, 78L, 79L, 80L, 81L, 82L, 83L, 84L, 85L, 86L,
87L, 89L, 90L, 91L, 92L, 93L, 94L, 95L, 96L, 97L, 1L, 2L, 3L,
4L, 5L, 6L, 8L, 9L, 10L, 11L), .Label = c("ALTURAS MUNICIPAL AIRPORT ",
"ARCATA AIRPORT ", "AUBURN MUNICIPAL AIRPORT ",
"BARSTOW-DAGGETT AIRPORT ", "BEALE AIR FORCE BASE ",
"BICYCLE LAKE FORT IRWIN AAF ", "BIG BEAR CITY AIRPORT ",
"BISHOP AIRPORT ", "BLUE CANYON - NYACK ARPT ",
"BRACKETT FIELD AIRPORT ", "BROWN FIELD MUNICIPAL ARPT ",
"BUCHANAN FIELD AIRPORT ", "BURBANK-GLENDALE-PASA ARPT ",
"CAMARILLO AIRPORT ", "CATALINA AIRPORT ",
"CHICO MUNICIPAL AIRPORT ", "CHINO AIRPORT ",
"COLUMBIA AIRPORT ", "CORONA MUNICIPAL AIRPORT ",
"DESERT RESORTS RGNL ARPT ", "FRESNO YOSEMITE INTERNATIONAL ",
"FULLERTON MUNICIPAL ARPT ", "GENERAL WILLIAM J. FOX AIRFIEL",
"GILLESPIE FIELD AIRPORT ", "GNOSS FIELD AIRPORT ",
"HALF MOON BAY AIRPORT ", "HANFORD MUNICIPAL AIRPORT ",
"HAYWARD EXECUTIVE AIRPORT ", "HOLLISTER MUNICIPAL AIRPORT ",
"IMPERIAL COUNTY AIRPORT ", "J. WAYNE APT-ORANGE CO APT ",
"JACK MCNAMARA FIELD ARPT ", "LAKE TAHOE AIRPORT ",
"LIVERMORE MUNICIPAL ARPT ", "LOMPOC AIRPORT ",
"LONG BEACH / DAUGHERTY FIELD /", "LONNIE POOL FIELD WEAVERVILLE ",
"LOS ALAMITOS AAF AIRPORT ", "LOS ANGELES INTERNATIONAL AIRP",
"MADERA MUNICIPAL AIRPORT ", "MAMMOTH YOSEMITE AIRPORT ",
"MCCLELLAN-PALOMAR AIRPORT ", "MDSTO CTY-CO H SHAM FD APT ",
"MEADOWS FIELD AIRPORT ", "METRO OAKLAND INTL AIRPORT ",
"MOFFETT FEDERAL AIRFLD APT ", "MOJAVE AIRPORT ",
"MONTEREY PENINSULA AIRPORT ", "MONTGOMERY FIELD AIRPORT ",
"MRCD MUNI/MACREADY FLD APT ", "N Y. MINETA SN JO INTL APT ",
"NALF/F. SHERMAN FLD ARPT ", "NAPA COUNTY AIRPORT ",
"NEEDLES AIRPORT ", "NEVADA COUNTY AIR PARK ",
"NRTHORP FLD/HATHRN MUNI AP ", "NUT TREE AIRPORT ",
"OCEANSIDE MUNICIPAL ARPT ", "ONTARIO INTERNATIONAL ARPT ",
"OROVILLE MUNICIPAL AIRPORT ", "OXNARD AIRPORT ",
"PALM SPRINGS INTL AIRPORT ", "PASO ROBLES MUNICIPAL ARPT ",
"PETALUMA MUNICIPAL AIRPORT ", "PLACERVILLE AIRPORT ",
"PORTERVILLE MUNICIPAL ARPT ", "RED BLUFF MUNICIPAL ARPT ",
"REDDING MUNICIPAL AIRPORT ", "RIVERSIDE MUNICIPAL ARPT ",
"ROHNERVILLE AIRPORT ", "SACRAMENTO EXECUTIVE AIRPORT ",
"SACRAMENTO INTL AIRPORT ", "SACRAMENTO MATHER AIRPORT ",
"SALINAS MUNICIPAL AIRPORT ", "SAN CARLOS AIRPORT ",
"SAN DIEGO INTERNATIONAL AIRPOR", "SAN FRANCISCO INTERNATIONAL AI",
"SAN LUIS CO REGIONAL ARPT ", "SANTA BARBARA MUNICIPAL AIRPOR",
"SANTA MARIA PUB/CAPT G A HANCO", "SANTA MONICA MUNI AIRPORT ",
"SANTA YNEZ AIRPORT ", "SISKIYOU COUNTY AIRPORT ",
"SOUTH COUNTY AIRPORT OF SANTA ", "SOUTHERN CALIFORNIA LOGISTICS ",
"STOCKTON METROPOLITAN AIRPORT ", "TEHACHAPI MUNICIPAL AIRPORT ",
"TRINITY CENTER AIRPORT ", "TRUCKEE-TAHOE AIRPORT ",
"UKIAH MUNICIPAL AIRPORT ", "UNIVERSITY AIRPORT ",
"VAN NUYS AIRPORT ", "VISALIA MUNICIPAL AIRPORT ",
"WATSONVILLE MUNICIPAL ARPT ", "WHITMAN AIRPORT ",
"YUBA COUNTRY AIRPORT ", "ZAMPERINI FIELD AIRPORT "
), class = "factor"), THI_Mean = c(79.3844194129895, 47.9523117470258,
69.6464461167776, 82.0675696692795, 63.0849350901548, 74.6708232862015,
100.850537076312, 74.8533760212847, 34.5523583142598, 54.1172284209428,
61.8235131140842, 49.6201577656199, 56.4969309931864, 53.8597392147438,
74.4432074783713, 46.837589950195, 64.6122605560279, NaN, 64.2706094582438,
64.337653411354, 52.2603460785691, 94.02576280203, 45.5601654090819,
55.0929418268904, 53.9012313226694, 68.3311804947574, 62.2675014979279,
63.654948977002, 56.0543459752155, 46.7621692338974, 81.9324399493885,
65.7709029415296, 54.8868962009991, 54.8363909689265, 56.8956616154744,
45.3120728155095, 54.5572278050216, 61.4196381163456, 92.7612688311182,
60.6665978683592, 53.7854832395117, 69.5295260343158, 61.3910104784722,
61.392754698528, 91.558905373712, 59.0707460260094, 51.5927530219057,
61.4922843624376, 60.0626408263436, 53.497914887709, 65.3054900336712,
67.0548958858132, 53.2591683242605, 65.8223643726204, 50.8970399973389,
51.2142000658336, 64.1480393838502, 61.4519885276234, 58.1897392633416,
71.5578784068541, 67.8204129142394, 71.601978756553, 58.5528424362473,
54.4067136577575, 56.5068448039998, 64.4002989004288, 59.0557207401386,
56.0689934073247, 55.1991038843164, 55.9231865358109, 61.5743119963135,
50.5540085434313, 53.5401150699814, 52.1935345223173, 61.3208716258652,
46.2035536331972, 80.6273562750492, 52.0374237358233, 87.4361185734358,
58.1078617762253, 91.5909088460021, 75.4622134300844, 62.8279716343081,
65.9364128813625, 50.0901362392351, 67.4125318829913, 56.3157177087556,
31.7479782943601, 60.6632690633163, 56.8037497381875, 67.61491611169,
51.1968002832636, 62.6089733566731, 74.4704531918518, 59.8603592342609,
66.0972847423599, 89.0724055109542, 70.6066530222754, 40.3039349206599,
53.3665353315505), THI_Min = c(52.4557405397089, 33.98, 45.3457093495371,
25.0225954802818, 51.1785334697297, 29.1511758012436, 28.2510119393172,
58.2138874079664, 0.502087905669441, 32.7983334230263, 46.6597894959148,
-11.1629416694279, 42.8124578159009, 27.47094073175, 54.9810099920306,
-7.16578623928804, 43.2890497467431, Inf, -3.70691599739962,
37.4694017813562, 16.7264957827268, 27.4482183124485, -13.5975437364969,
35.6903112057144, 30.2000000000001, 59.4304095606214, 8.7735370709971,
-2.80405898293145, 35.7697087097596, 33.98, 68.3855026732981,
25.3133696920698, 28.556343667171, 33.5599938044734, 37.8359070192615,
30.9228228628589, 46.3999999999999, 40.8603031344418, 64.2502080724629,
52.1493991780154, 36.1767259310531, 42.10373886285, 53.6749547440094,
53.5952658714694, 33.7332666403973, 45.3594438239759, 10.7322063083914,
37.9467060766937, 47.289117456444, 33.0925630861752, 1.90029175361333,
44.2510827881006, 40.9199096994557, 42.3990582610788, 37.4, -10.139479142008,
46.6597894959148, 53.6438500522523, 3.61411799940322, 8.37516788708545,
49.2873234865599, 27.5338335958916, 10.1664311833085, -10.7072039412661,
42.8465883329053, 52.4497314578047, 45.3457093495371, 17.259303962137,
46.1362993029647, 37.7879989441233, 55.2183148794405, 10.129529761709,
46.1362993029647, 10.129529761709, 49.8796278533633, -7.00765625628347,
60.9804990557964, 30.500203826155, 33.7332666403973, 41.793006926767,
40.7748936693757, 49.0989499378863, 32.0322286350935, 32, -8.27786646700535,
47.2702704566806, 19.2826923623029, 1.44285466185288, 44.3337809084761,
46.3999999999999, 35.3138652463278, 41.7170489863234, 33.7351586623816,
30.4879219330503, 43.9911273883178, 44.7926129254152, 32.7467168398368,
58.2191654136387, -12.498112226052, 22.9228778307818), THI_Max = c(109.591786682477,
60.2333967252995, 82.8768296235874, 107.612764911568, 77.0245053208592,
102.168364973031, 139.996109979679, 82.6768921516963, 68.6026287228502,
78.3456284931154, 72.0533015440348, 84.0059030769878, 67.0179354243032,
70.0490739062396, 82.2054720080854, 68.2291538760102, 79.8358713003518,
-Inf, 99.2586259600651, 75.7556525031902, 67.0179354243032, 117.52943488297,
73.5830893311879, 68.65063426021, 69.4335307001204, 77.3076678339767,
76.9529740661977, 100.379670258497, 62.5970088690744, 64.009438459283,
98.9476881025318, 78.3456284931154, 66.3151890399475, 60.9878582718125,
82.2054720080854, 58.3647202290469, 62.1651540978214, 73.5529769693358,
112.070201509673, 72.6703589319581, 71.8716497050957, 81.3559067801607,
68.7415770349638, 66.0921490462558, 122.093934473517, 76.1709284501561,
70.6341779681281, 73.5529769693358, 72.0533015440348, 72.0533015440348,
92.7299425243837, 84.9072999765804, 61.9348463512754, 84.4328754193385,
68.6207402936826, 77.8815726214716, 71.8716497050957, 79.7562300236269,
95.3037028368627, 96.2490043020463, 87.3763071289113, 85.9105487050643,
76.4751770487016, 86.4975370577019, 73.877504195677, 77.521576933463,
75.1948807718176, 78.6082511434511, 66.3151890399475, 71.9551337820645,
67.028237117519, 63.6006960092746, 59.249786941551, 63.8719613651346,
70.6341779681281, 73.1787331352403, 101.186226933919, 73.5830893311879,
115.648856045281, 72.6703589319581, 121.651090979094, 99.5352109263185,
81.5737450432125, 85.3244678886861, 76.5427815713664, 79.1421115384682,
71.3180616536051, 82.0833455621802, 81.9382809187063, 71.0046150736307,
89.7955910777776, 60.4042257053666, 73.1787331352403, 94.3704734177277,
71.4307310204537, 85.4942371284446, 120.376615578723, 79.9743241715952,
80.7722101628022, 73.5571882695921), THIThreshold_1 = c(68, 68,
68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68,
68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68,
68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68,
68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68,
68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68,
68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68, 68,
68, 68), THIThreshold_2 = c(72, 72, 72, 72, 72, 72, 72, 72, 72,
72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72,
72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72,
72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72,
72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72,
72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72,
72, 72, 72, 72, 72, 72, 72, 72, 72, 72, 72)), class = c("grouped_df",
"tbl_df", "tbl", "data.frame"), row.names = c(NA, -100L), groups = structure(list(
Dates = structure(c(16071, 16072), class = "Date"), .rows = list(
1:90, 91:100)), row.names = c(NA, -2L), class = c("tbl_df",
"tbl", "data.frame"), .drop = TRUE))
И мне нужно применить этот код:
period <- 24
amplitude <- (Dados2$THI_Max - Dados2$THI_Min)/2
x1 <- asin ((Dados2$THIThreshold_1 - Dados2$THI_Mean)/amplitude)
x2 <- pi - x1
x11 <- asin ((Dados2$THI_Mean - Dados2$THIThreshold_1)/amplitude)
x22 <- pi - x11
ad <- pi + asin ((Dados2$THI_Mean - Dados2$THIThreshold_1)/amplitude)
ad2 <- pi + asin((Dados2$THIThreshold_1 - Dados2$THI_Mean)/amplitude)
if (Dados2$THI_Mean >= Dados2$THI_Max) {
THI_Load <- 0
} else if (Dados2$THIThreshold_1 <= Dados2$THI_Min){
THI_Load = period * (Dados2$THI_Mean - Dados2$THIThreshold_1)
} else if (Dados2$THIThreshold_1 >= Dados2$THI_Mean){
THI_Load = (cos(x1) - cos(x2)) * amplitude * period/2/pi - (x2 - x1) * period/2/pi*(Dados2$THIThreshold_1-Dados2$THI_Mean)
} else {
THI_Load = amplitude * period/pi + (Dados2$THI_Mean - Dados2$THIThreshold_1) * period/2 + (Dados2$THI_Mean - Dados2$THIThreshold_1) * ((ad-pi) * period/pi) - (cos(ad) - cos(pi)) * amplitude * period/pi
}
Но я получаю это сообщение:
Warning messages:
1: In if (Dados2$THI_Mean >= Dados2$THI_Max) { :
the condition has length > 1 and only the first element will be used
2: In if (Dados2$THIThreshold_1 <= Dados2$THI_Min) { :
the condition has length > 1 and only the first element will be used
3: In if (Dados2$THIThreshold_1 >= Dados2$THI_Mean) { :
the condition has length > 1 and only the first element will be used
Как я могу это исправить?
И еще вопрос, возможно ли создать столбец в data.frame с результатом?