##Triple Helix Calculation for 4 Variables #Ed Lefley #Shawn Leu #2021 12 02 #2021 12 08 #2021 12 10 #2022 01 19 Added Overseas and changed to using full dataset #2022 03 24 Changed to match Leydesdorff 2021 revision #Based on Leydesdorff's original calculations #Expanded using Hu et al (2021) library(stringr) library(data.table) Ag7_5H_2003 <- subset(Ag7_TH_M_TH_SS_5, Publication_Year == 2003) #With the newly created set, now create a key based on the set patterns: Entropy_2003_5_Ag7 <- cbind.data.frame(Ag7_5H_2003$TH_U, Ag7_5H_2003$TH_I, Ag7_5H_2003$TH_G, Ag7_5H_2003$TH_S, Ag7_5H_2003$TH_F) #write.csv(Entropy_2003, "data.txt", col.names =FALSE) #Set the total number of rows: denom_2003_5_Ag7 <- nrow(Ag7_5H_2003) #Entropy calculation for UI covers #(U,I,G,O) #Using the Ag7_5H_2003$TH_U, Ag7_5H_2003$TH_I, Ag7_5H_2003$TH_G, Ag7_5H_2003$TH_S and Ag7_5H_2003$TH_F columns, to replicate the 2021 revision of the TH4.exe model the following is performed: #Convert character to numeric Ag7_5H_2003$TH_U <- as.numeric(Ag7_5H_2003$TH_U) Ag7_5H_2003$TH_I <- as.numeric(Ag7_5H_2003$TH_I) Ag7_5H_2003$TH_G <- as.numeric(Ag7_5H_2003$TH_G) Ag7_5H_2003$TH_S <- as.numeric(Ag7_5H_2003$TH_S) Ag7_5H_2003$TH_F <- as.numeric(Ag7_5H_2003$TH_F) #U Single PU_2003_5LL_Ag7 <- sum((Ag7_5H_2003$TH_U/denom_2003_5_Ag7)) PNU_2003_5LL_Ag7 <- sum((denom_2003_5_Ag7-sum(Ag7_5H_2003$TH_U))/denom_2003_5_Ag7) #I Single PI_2003_5LL_Ag7 <- sum((Ag7_5H_2003$TH_I/denom_2003_5_Ag7)) PNI_2003_5LL_Ag7 <- sum((denom_2003_5_Ag7-sum(Ag7_5H_2003$TH_I))/denom_2003_5_Ag7) #G Single PG_2003_5LL_Ag7 <- sum((Ag7_5H_2003$TH_G/denom_2003_5_Ag7)) PNG_2003_5LL_Ag7 <- sum((denom_2003_5_Ag7-sum(Ag7_5H_2003$TH_G))/denom_2003_5_Ag7) #S Single PS_2003_5LL_Ag7 <- sum((Ag7_5H_2003$TH_S/denom_2003_5_Ag7)) PNS_2003_5LL_Ag7 <- sum((denom_2003_5_Ag7-sum(Ag7_5H_2003$TH_S))/denom_2003_5_Ag7) #F Single PF_2003_5LL_Ag7 <- sum((Ag7_5H_2003$TH_F/denom_2003_5_Ag7)) PNF_2003_5LL_Ag7 <- sum((denom_2003_5_Ag7-sum(Ag7_5H_2003$TH_F))/denom_2003_5_Ag7) #Pairs (10) #Multiply denom_2003 by 2 for the pairs denom_2003_5_Ag7_2 <- denom_2003_5_Ag7*2 #UI PUI1_2003_5LL_Ag7 <- sum((Ag7_5H_2003$TH_U+Ag7_5H_2003$TH_I)/denom_2003_5_Ag7_2) PUI2_2003_5LL_Ag7 <- sum((denom_2003_5_Ag7_2-sum(Ag7_5H_2003$TH_U+Ag7_5H_2003$TH_I))/denom_2003_5_Ag7_2) #UG Pairs PUG1_2003_5LL_Ag7 <- sum((Ag7_5H_2003$TH_U+Ag7_5H_2003$TH_G)/denom_2003_5_Ag7_2) PUG2_2003_5LL_Ag7 <- sum((denom_2003_5_Ag7_2-sum(Ag7_5H_2003$TH_U+Ag7_5H_2003$TH_G))/denom_2003_5_Ag7_2) #US Pairs PUS1_2003_5LL_Ag7 <- sum((Ag7_5H_2003$TH_U+Ag7_5H_2003$TH_S)/denom_2003_5_Ag7_2) PUS2_2003_5LL_Ag7 <- sum((denom_2003_5_Ag7_2-sum(Ag7_5H_2003$TH_U+Ag7_5H_2003$TH_S))/denom_2003_5_Ag7_2) #UF PUF1_2003_5LL_Ag7 <- sum((Ag7_5H_2003$TH_U+Ag7_5H_2003$TH_F)/denom_2003_5_Ag7_2) PUF2_2003_5LL_Ag7 <- sum((denom_2003_5_Ag7_2-sum(Ag7_5H_2003$TH_U+Ag7_5H_2003$TH_F))/denom_2003_5_Ag7_2) #IG Pairs PIG1_2003_5LL_Ag7 <- sum((Ag7_5H_2003$TH_I+Ag7_5H_2003$TH_G)/denom_2003_5_Ag7_2) PIG2_2003_5LL_Ag7 <- sum((denom_2003_5_Ag7_2-sum(Ag7_5H_2003$TH_I+Ag7_5H_2003$TH_G))/denom_2003_5_Ag7_2) #IS Pairs PIS1_2003_5LL_Ag7 <- sum((Ag7_5H_2003$TH_I+Ag7_5H_2003$TH_S)/denom_2003_5_Ag7_2) PIS2_2003_5LL_Ag7 <- sum((denom_2003_5_Ag7_2-sum(Ag7_5H_2003$TH_I+Ag7_5H_2003$TH_S))/denom_2003_5_Ag7_2) #IF Pairs PIF1_2003_5LL_Ag7 <- sum((Ag7_5H_2003$TH_I+Ag7_5H_2003$TH_F)/denom_2003_5_Ag7_2) PIF2_2003_5LL_Ag7 <- sum((denom_2003_5_Ag7_2-sum(Ag7_5H_2003$TH_I+Ag7_5H_2003$TH_F))/denom_2003_5_Ag7_2) #GS Pairs PGS1_2003_5LL_Ag7 <- sum((Ag7_5H_2003$TH_G+Ag7_5H_2003$TH_S)/denom_2003_5_Ag7_2) PGS2_2003_5LL_Ag7 <- sum((denom_2003_5_Ag7_2-sum(Ag7_5H_2003$TH_G+Ag7_5H_2003$TH_S))/denom_2003_5_Ag7_2) #GF PGF1_2003_5LL_Ag7 <- sum((Ag7_5H_2003$TH_G+Ag7_5H_2003$TH_F)/denom_2003_5_Ag7_2) PGF2_2003_5LL_Ag7 <- sum((denom_2003_5_Ag7_2-sum(Ag7_5H_2003$TH_G+Ag7_5H_2003$TH_F))/denom_2003_5_Ag7_2) #SF PSF1_2003_5LL_Ag7 <- sum((Ag7_5H_2003$TH_S+Ag7_5H_2003$TH_F)/denom_2003_5_Ag7_2) PSF2_2003_5LL_Ag7 <- sum((denom_2003_5_Ag7_2-sum(Ag7_5H_2003$TH_S+Ag7_5H_2003$TH_F))/denom_2003_5_Ag7_2) #Multiply denom_2003 by 3 for the Triples (10) denom_2003_5_Ag7_3 <- denom_2003_5_Ag7*3 #UIG Triples PUIG1_2003_5LL_Ag7<- sum((Ag7_5H_2003$TH_U+Ag7_5H_2003$TH_I+Ag7_5H_2003$TH_G)/denom_2003_5_Ag7_3) PUIG2_2003_5LL_Ag7<- sum((denom_2003_5_Ag7_3-sum(Ag7_5H_2003$TH_U+Ag7_5H_2003$TH_I+Ag7_5H_2003$TH_G))/denom_2003_5_Ag7_3) #UIS Triples PUIS1_2003_5LL_Ag7<- sum((Ag7_5H_2003$TH_U+Ag7_5H_2003$TH_I+Ag7_5H_2003$TH_S)/denom_2003_5_Ag7_3) #P(1,1,1,*) PUIS2_2003_5LL_Ag7<- sum((denom_2003_5_Ag7_3-sum(Ag7_5H_2003$TH_U+Ag7_5H_2003$TH_I+Ag7_5H_2003$TH_S))/denom_2003_5_Ag7_3) #UIF Triples PUIF1_2003_5LL_Ag7<- sum((Ag7_5H_2003$TH_U+Ag7_5H_2003$TH_I+Ag7_5H_2003$TH_F)/denom_2003_5_Ag7_3) PUIF2_2003_5LL_Ag7<- sum((denom_2003_5_Ag7_3-sum(Ag7_5H_2003$TH_U+Ag7_5H_2003$TH_I+Ag7_5H_2003$TH_F))/denom_2003_5_Ag7_3) #USG Triples PUGS1_2003_5LL_Ag7<- sum((Ag7_5H_2003$TH_U+Ag7_5H_2003$TH_S+Ag7_5H_2003$TH_G)/denom_2003_5_Ag7_3) #P(1,1,1,*) PUGS2_2003_5LL_Ag7<- sum((denom_2003_5_Ag7_3-sum(Ag7_5H_2003$TH_U+Ag7_5H_2003$TH_S+Ag7_5H_2003$TH_G))/denom_2003_5_Ag7_3) #UGF PUGF1_2003_5LL_Ag7<- sum((Ag7_5H_2003$TH_U+Ag7_5H_2003$TH_F+Ag7_5H_2003$TH_G)/denom_2003_5_Ag7_3) PUGF2_2003_5LL_Ag7<- sum((denom_2003_5_Ag7_3-sum(Ag7_5H_2003$TH_U+Ag7_5H_2003$TH_F+Ag7_5H_2003$TH_G))/denom_2003_5_Ag7_3) #USF PUSF1_2003_5LL_Ag7<- sum((Ag7_5H_2003$TH_U+Ag7_5H_2003$TH_S+Ag7_5H_2003$TH_F)/denom_2003_5_Ag7_3) PUSF2_2003_5LL_Ag7<- sum((denom_2003_5_Ag7_3-sum(Ag7_5H_2003$TH_U+Ag7_5H_2003$TH_S+Ag7_5H_2003$TH_F))/denom_2003_5_Ag7_3) #IGS Triples PIGS1_2003_5LL_Ag7<- sum((Ag7_5H_2003$TH_I+Ag7_5H_2003$TH_G+Ag7_5H_2003$TH_S)/denom_2003_5_Ag7_3) #P(1,1,1,*) PIGS2_2003_5LL_Ag7<- sum((denom_2003_5_Ag7_3-sum(Ag7_5H_2003$TH_I+Ag7_5H_2003$TH_G+Ag7_5H_2003$TH_S))/denom_2003_5_Ag7_3) #IGF PIGF1_2003_5LL_Ag7<- sum((Ag7_5H_2003$TH_I+Ag7_5H_2003$TH_G+Ag7_5H_2003$TH_F)/denom_2003_5_Ag7_3) PIGF2_2003_5LL_Ag7<- sum((denom_2003_5_Ag7_3-sum(Ag7_5H_2003$TH_I+Ag7_5H_2003$TH_G+Ag7_5H_2003$TH_F))/denom_2003_5_Ag7_3) #ISF PISF1_2003_5LL_Ag7<- sum((Ag7_5H_2003$TH_I+Ag7_5H_2003$TH_S+Ag7_5H_2003$TH_F)/denom_2003_5_Ag7_3) PISF2_2003_5LL_Ag7<- sum((denom_2003_5_Ag7_3-sum(Ag7_5H_2003$TH_I+Ag7_5H_2003$TH_S+Ag7_5H_2003$TH_F))/denom_2003_5_Ag7_3) #GSF PGSF1_2003_5LL_Ag7<- sum((Ag7_5H_2003$TH_G+Ag7_5H_2003$TH_S+Ag7_5H_2003$TH_F)/denom_2003_5_Ag7_3) PGSF2_2003_5LL_Ag7<- sum((denom_2003_5_Ag7_3-sum(Ag7_5H_2003$TH_G+Ag7_5H_2003$TH_S+Ag7_5H_2003$TH_F))/denom_2003_5_Ag7_3) #Multiply denom_2003 by 4 for the quads (5) denom_2003_5_Ag7_4 <- denom_2003_5_Ag7*4 #UIGS PUIGS1_2003_5LL_Ag7 <- sum((Ag7_5H_2003$TH_U+Ag7_5H_2003$TH_I+Ag7_5H_2003$TH_G+Ag7_5H_2003$TH_S)/denom_2003_5_Ag7_4) PUIGS2_2003_5LL_Ag7 <- sum((denom_2003_5_Ag7_4-sum(Ag7_5H_2003$TH_U+Ag7_5H_2003$TH_I+Ag7_5H_2003$TH_G+Ag7_5H_2003$TH_S))/denom_2003_5_Ag7_4) #UIGF PUIGF1_2003_5LL_Ag7 <- sum((Ag7_5H_2003$TH_U+Ag7_5H_2003$TH_I+Ag7_5H_2003$TH_G+Ag7_5H_2003$TH_F)/denom_2003_5_Ag7_4) PUIGF2_2003_5LL_Ag7 <- sum((denom_2003_5_Ag7_4-sum(Ag7_5H_2003$TH_U+Ag7_5H_2003$TH_I+Ag7_5H_2003$TH_G+Ag7_5H_2003$TH_F))/denom_2003_5_Ag7_4) #UISF PUISF1_2003_5LL_Ag7 <- sum((Ag7_5H_2003$TH_U+Ag7_5H_2003$TH_I+Ag7_5H_2003$TH_S+Ag7_5H_2003$TH_F)/denom_2003_5_Ag7_4) PUISF2_2003_5LL_Ag7 <- sum((denom_2003_5_Ag7_4-sum(Ag7_5H_2003$TH_U+Ag7_5H_2003$TH_I+Ag7_5H_2003$TH_S+Ag7_5H_2003$TH_F))/denom_2003_5_Ag7_4) #UGSF PUGSF1_2003_5LL_Ag7 <- sum((Ag7_5H_2003$TH_U+Ag7_5H_2003$TH_G+Ag7_5H_2003$TH_S+Ag7_5H_2003$TH_F)/denom_2003_5_Ag7_4) PUGSF2_2003_5LL_Ag7 <- sum((denom_2003_5_Ag7_4-sum(Ag7_5H_2003$TH_U+Ag7_5H_2003$TH_G+Ag7_5H_2003$TH_S+Ag7_5H_2003$TH_F))/denom_2003_5_Ag7_4) #IGSF PIGSF1_2003_5LL_Ag7 <- sum((Ag7_5H_2003$TH_I+Ag7_5H_2003$TH_G+Ag7_5H_2003$TH_S+Ag7_5H_2003$TH_F)/denom_2003_5_Ag7_4) PIGSF2_2003_5LL_Ag7 <- sum((denom_2003_5_Ag7_4-sum(Ag7_5H_2003$TH_I+Ag7_5H_2003$TH_G+Ag7_5H_2003$TH_S+Ag7_5H_2003$TH_F))/denom_2003_5_Ag7_4) #Quintuple (1) denom_2003_5_Ag7_5 <- denom_2003_5_Ag7*5 PUIGSF1_2003_5LL_Ag7 <- sum((Ag7_5H_2003$TH_U+Ag7_5H_2003$TH_I+Ag7_5H_2003$TH_G+Ag7_5H_2003$TH_S+Ag7_5H_2003$TH_F)/denom_2003_5_Ag7_5) PUIGSF2_2003_5LL_Ag7 <- sum((denom_2003_5_Ag7_5-sum(Ag7_5H_2003$TH_U+Ag7_5H_2003$TH_I+Ag7_5H_2003$TH_G+Ag7_5H_2003$TH_S+Ag7_5H_2003$TH_F))/denom_2003_5_Ag7_5) #Shannon's Entropy calculations: HU_2003_5LL_Ag7 <- (-PU_2003_5LL_Ag7*(log2(PU_2003_5LL_Ag7))) - (PNU_2003_5LL_Ag7*(log2(PNU_2003_5LL_Ag7))) HI_2003_5LL_Ag7 <- (-PI_2003_5LL_Ag7*(log2(PI_2003_5LL_Ag7))) - (PNI_2003_5LL_Ag7*(log2(PNI_2003_5LL_Ag7))) HG_2003_5LL_Ag7 <- (-PG_2003_5LL_Ag7*(log2(PG_2003_5LL_Ag7))) - (PNG_2003_5LL_Ag7*(log2(PNG_2003_5LL_Ag7))) HS_2003_5LL_Ag7 <- (-PS_2003_5LL_Ag7*(log2(PS_2003_5LL_Ag7))) - (PNS_2003_5LL_Ag7*(log2(PNS_2003_5LL_Ag7))) HF_2003_5LL_Ag7 <- (-PF_2003_5LL_Ag7*(log2(PF_2003_5LL_Ag7))) - (PNF_2003_5LL_Ag7*(log2(PNF_2003_5LL_Ag7))) #Pairs HUI_2003_5LL_Ag7 <- -(PUI1_2003_5LL_Ag7*(log2(PUI1_2003_5LL_Ag7)))-(PUI2_2003_5LL_Ag7*(log2(PUI2_2003_5LL_Ag7))) HUG_2003_5LL_Ag7 <- -(PUG1_2003_5LL_Ag7*(log2(PUG1_2003_5LL_Ag7)))-(PUG2_2003_5LL_Ag7*(log2(PUG2_2003_5LL_Ag7))) HUS_2003_5LL_Ag7 <- -(PUS1_2003_5LL_Ag7*(log2(PUS1_2003_5LL_Ag7)))-(PUS2_2003_5LL_Ag7*(log2(PUS2_2003_5LL_Ag7))) HUF_2003_5LL_Ag7 <- -(PUF1_2003_5LL_Ag7*(log2(PUF1_2003_5LL_Ag7)))-(PUF2_2003_5LL_Ag7*(log2(PUF2_2003_5LL_Ag7))) HIG_2003_5LL_Ag7 <- -(PIG1_2003_5LL_Ag7*(log2(PIG1_2003_5LL_Ag7)))-(PIG2_2003_5LL_Ag7*(log2(PIG2_2003_5LL_Ag7))) HIS_2003_5LL_Ag7 <- -(PIS1_2003_5LL_Ag7*(log2(PIS1_2003_5LL_Ag7)))-(PIS2_2003_5LL_Ag7*(log2(PIS2_2003_5LL_Ag7))) HIF_2003_5LL_Ag7 <- -(PIF1_2003_5LL_Ag7*(log2(PIF1_2003_5LL_Ag7)))-(PIF2_2003_5LL_Ag7*(log2(PIF2_2003_5LL_Ag7))) HGS_2003_5LL_Ag7 <- -(PGS1_2003_5LL_Ag7*(log2(PGS1_2003_5LL_Ag7)))-(PGS2_2003_5LL_Ag7*(log2(PGS2_2003_5LL_Ag7))) HGF_2003_5LL_Ag7 <- -(PGF1_2003_5LL_Ag7*(log2(PGF1_2003_5LL_Ag7)))-(PGF2_2003_5LL_Ag7*(log2(PGF2_2003_5LL_Ag7))) HSF_2003_5LL_Ag7 <- -(PSF1_2003_5LL_Ag7*(log2(PSF1_2003_5LL_Ag7)))-(PSF2_2003_5LL_Ag7*(log2(PSF2_2003_5LL_Ag7))) #Triples HUIG_2003_5LL_Ag7 <- (-(PUIG1_2003_5LL_Ag7)*(log2(PUIG1_2003_5LL_Ag7))) - ((PUIG2_2003_5LL_Ag7)*(log2(PUIG2_2003_5LL_Ag7))) HUIS_2003_5LL_Ag7 <- (-(PUIS1_2003_5LL_Ag7)*(log2(PUIS1_2003_5LL_Ag7))) - ((PUIS2_2003_5LL_Ag7)*(log2(PUIS2_2003_5LL_Ag7))) HUIF_2003_5LL_Ag7 <- (-(PUIF1_2003_5LL_Ag7)*(log2(PUIF1_2003_5LL_Ag7))) - ((PUIF2_2003_5LL_Ag7)*(log2(PUIF2_2003_5LL_Ag7))) HUGS_2003_5LL_Ag7 <- (-(PUGS1_2003_5LL_Ag7)*(log2(PUGS1_2003_5LL_Ag7))) - ((PUGS2_2003_5LL_Ag7)*(log2(PUGS2_2003_5LL_Ag7))) HUGF_2003_5LL_Ag7 <- (-(PUGF1_2003_5LL_Ag7)*(log2(PUGF1_2003_5LL_Ag7))) - ((PUGF2_2003_5LL_Ag7)*(log2(PUGF2_2003_5LL_Ag7))) HUSF_2003_5LL_Ag7 <- (-(PUSF1_2003_5LL_Ag7)*(log2(PUSF1_2003_5LL_Ag7))) - ((PUSF2_2003_5LL_Ag7)*(log2(PUSF2_2003_5LL_Ag7))) HIGS_2003_5LL_Ag7 <- (-(PIGS1_2003_5LL_Ag7)*(log2(PIGS1_2003_5LL_Ag7))) - ((PIGS2_2003_5LL_Ag7)*(log2(PIGS2_2003_5LL_Ag7))) HIGF_2003_5LL_Ag7 <- (-(PIGF1_2003_5LL_Ag7)*(log2(PIGF1_2003_5LL_Ag7))) - ((PIGF2_2003_5LL_Ag7)*(log2(PIGF2_2003_5LL_Ag7))) HISF_2003_5LL_Ag7 <- (-(PISF1_2003_5LL_Ag7)*(log2(PISF1_2003_5LL_Ag7))) - ((PISF2_2003_5LL_Ag7)*(log2(PISF2_2003_5LL_Ag7))) HGSF_2003_5LL_Ag7 <- (-(PGSF1_2003_5LL_Ag7)*(log2(PGSF1_2003_5LL_Ag7))) - ((PGSF2_2003_5LL_Ag7)*(log2(PGSF2_2003_5LL_Ag7))) #Quads HUIGS_2003_5LL_Ag7 <- (-(PUIGS1_2003_5LL_Ag7)*(log2(PUIGS1_2003_5LL_Ag7))) - ((PUIGS2_2003_5LL_Ag7)*(log2(PUIGS2_2003_5LL_Ag7))) HUIGF_2003_5LL_Ag7 <- (-(PUIGF1_2003_5LL_Ag7)*(log2(PUIGF1_2003_5LL_Ag7))) - ((PUIGF2_2003_5LL_Ag7)*(log2(PUIGF2_2003_5LL_Ag7))) HUISF_2003_5LL_Ag7 <- (-(PUISF1_2003_5LL_Ag7)*(log2(PUISF1_2003_5LL_Ag7))) - ((PUISF2_2003_5LL_Ag7)*(log2(PUISF2_2003_5LL_Ag7))) HUGSF_2003_5LL_Ag7 <- (-(PUGSF1_2003_5LL_Ag7)*(log2(PUGSF1_2003_5LL_Ag7))) - ((PUGSF2_2003_5LL_Ag7)*(log2(PUGSF2_2003_5LL_Ag7))) HIGSF_2003_5LL_Ag7 <- (-(PIGSF1_2003_5LL_Ag7)*(log2(PIGSF1_2003_5LL_Ag7))) - ((PIGSF2_2003_5LL_Ag7)*(log2(PIGSF2_2003_5LL_Ag7))) #QUintuple HUIGSF_2003_5LL_Ag7<- (-(PUIGSF1_2003_5LL_Ag7)*(log2(PUIGSF1_2003_5LL_Ag7))) - ((PUIGSF2_2003_5LL_Ag7)*(log2(PUIGSF2_2003_5LL_Ag7))) #Conversion to Transmission (in millibits) #BASED ON "A Routine for Measuring Synergy in University-Industry-Government Relations: Mutual Information as a Triple-Helix and Quadruple-Helix Indicator" Loet Leydesdorff, Han Woo Park, & Balazs Lengyelc #PAIRS TUI_2003_5LL_Ag7 <- (HU_2003_5LL_Ag7 + HI_2003_5LL_Ag7 - HUI_2003_5LL_Ag7)#*1000 TUG_2003_5LL_Ag7 <- (HU_2003_5LL_Ag7 + HG_2003_5LL_Ag7 - HUG_2003_5LL_Ag7)#*1000 TUS_2003_5LL_Ag7 <- (HU_2003_5LL_Ag7 + HS_2003_5LL_Ag7 - HUS_2003_5LL_Ag7)#*1000 TUF_2003_5LL_Ag7 <- (HU_2003_5LL_Ag7 + HF_2003_5LL_Ag7 - HUF_2003_5LL_Ag7)#*1000 TIG_2003_5LL_Ag7 <- (HI_2003_5LL_Ag7 + HG_2003_5LL_Ag7 - HIG_2003_5LL_Ag7)#*1000 TIS_2003_5LL_Ag7 <- (HI_2003_5LL_Ag7 + HS_2003_5LL_Ag7 - HIS_2003_5LL_Ag7)#*1000 TIF_2003_5LL_Ag7 <- (HI_2003_5LL_Ag7 + HF_2003_5LL_Ag7 - HIF_2003_5LL_Ag7)#*1000 TGS_2003_5LL_Ag7 <- (HG_2003_5LL_Ag7 + HS_2003_5LL_Ag7 - HGS_2003_5LL_Ag7)#*1000 TGF_2003_5LL_Ag7 <- (HG_2003_5LL_Ag7 + HF_2003_5LL_Ag7 - HGF_2003_5LL_Ag7)#*1000 TSF_2003_5LL_Ag7 <- (HS_2003_5LL_Ag7 + HF_2003_5LL_Ag7 - HSF_2003_5LL_Ag7)#*1000 #TRIPLE TUIG_2003_5LL_Ag7 <- (HU_2003_5LL_Ag7 + HI_2003_5LL_Ag7 + HG_2003_5LL_Ag7 - HUI_2003_5LL_Ag7 - HUG_2003_5LL_Ag7 - HIG_2003_5LL_Ag7 + HUIG_2003_5LL_Ag7)#*1000 TUIS_2003_5LL_Ag7 <- (HU_2003_5LL_Ag7 + HI_2003_5LL_Ag7 + HS_2003_5LL_Ag7 - HUI_2003_5LL_Ag7 - HUS_2003_5LL_Ag7 - HIS_2003_5LL_Ag7 + HUIS_2003_5LL_Ag7)#*1000 TUIF_2003_5LL_Ag7 <- (HU_2003_5LL_Ag7 + HI_2003_5LL_Ag7 + HF_2003_5LL_Ag7 - HUI_2003_5LL_Ag7 - HUF_2003_5LL_Ag7 - HIF_2003_5LL_Ag7 + HUIF_2003_5LL_Ag7)#*1000 TUGS_2003_5LL_Ag7 <- (HU_2003_5LL_Ag7 + HG_2003_5LL_Ag7 + HS_2003_5LL_Ag7 - HUG_2003_5LL_Ag7 - HUS_2003_5LL_Ag7 - HGS_2003_5LL_Ag7 + HUGS_2003_5LL_Ag7)#*1000 TUGF_2003_5LL_Ag7 <- (HU_2003_5LL_Ag7 + HG_2003_5LL_Ag7 + HF_2003_5LL_Ag7 - HUG_2003_5LL_Ag7 - HUF_2003_5LL_Ag7 - HGF_2003_5LL_Ag7 + HUGF_2003_5LL_Ag7)#*1000 TUSF_2003_5LL_Ag7 <- (HU_2003_5LL_Ag7 + HS_2003_5LL_Ag7 + HF_2003_5LL_Ag7 - HUS_2003_5LL_Ag7 - HUF_2003_5LL_Ag7 - HSF_2003_5LL_Ag7 + HUSF_2003_5LL_Ag7)#*1000 TIGS_2003_5LL_Ag7 <- (HI_2003_5LL_Ag7 + HG_2003_5LL_Ag7 + HS_2003_5LL_Ag7 - HIG_2003_5LL_Ag7 - HIS_2003_5LL_Ag7 - HGS_2003_5LL_Ag7 + HIGS_2003_5LL_Ag7)#*1000 TIGF_2003_5LL_Ag7 <- (HI_2003_5LL_Ag7 + HG_2003_5LL_Ag7 + HF_2003_5LL_Ag7 - HIG_2003_5LL_Ag7 - HIF_2003_5LL_Ag7 - HGF_2003_5LL_Ag7 + HIGF_2003_5LL_Ag7)#*1000 TISF_2003_5LL_Ag7 <- (HI_2003_5LL_Ag7 + HS_2003_5LL_Ag7 + HF_2003_5LL_Ag7 - HIS_2003_5LL_Ag7 - HIF_2003_5LL_Ag7 - HSF_2003_5LL_Ag7 + HISF_2003_5LL_Ag7)#*1000 TGSF_2003_5LL_Ag7 <- (HG_2003_5LL_Ag7 + HG_2003_5LL_Ag7 + HF_2003_5LL_Ag7 - HGS_2003_5LL_Ag7 - HGF_2003_5LL_Ag7 - HGF_2003_5LL_Ag7 + HGSF_2003_5LL_Ag7)#*1000 #QUAD TUIGS_2003_5LL_Ag7 <- (HU_2003_5LL_Ag7 + HI_2003_5LL_Ag7 + HG_2003_5LL_Ag7 + HS_2003_5LL_Ag7 - HUI_2003_5LL_Ag7 - HUG_2003_5LL_Ag7 - HUS_2003_5LL_Ag7 - HIG_2003_5LL_Ag7 - HIS_2003_5LL_Ag7 - HGS_2003_5LL_Ag7 + HUIG_2003_5LL_Ag7 + HUIS_2003_5LL_Ag7 + HUGS_2003_5LL_Ag7 + HIGS_2003_5LL_Ag7 - HUIGS_2003_5LL_Ag7)#*1000 TUIGF_2003_5LL_Ag7 <- (HU_2003_5LL_Ag7 + HI_2003_5LL_Ag7 + HG_2003_5LL_Ag7 + HF_2003_5LL_Ag7 - HUI_2003_5LL_Ag7 - HUG_2003_5LL_Ag7 - HUF_2003_5LL_Ag7 - HIG_2003_5LL_Ag7 - HIF_2003_5LL_Ag7 - HGF_2003_5LL_Ag7 + HUIG_2003_5LL_Ag7 + HUIF_2003_5LL_Ag7 + HUGF_2003_5LL_Ag7 + HIGF_2003_5LL_Ag7 - HUIGF_2003_5LL_Ag7)#*1000 TUISF_2003_5LL_Ag7 <- (HU_2003_5LL_Ag7 + HI_2003_5LL_Ag7 + HS_2003_5LL_Ag7 + HF_2003_5LL_Ag7 - HUI_2003_5LL_Ag7 - HUS_2003_5LL_Ag7 - HUF_2003_5LL_Ag7 - HIS_2003_5LL_Ag7 - HIF_2003_5LL_Ag7 - HSF_2003_5LL_Ag7 + HUIS_2003_5LL_Ag7 + HUIF_2003_5LL_Ag7 + HUSF_2003_5LL_Ag7 + HISF_2003_5LL_Ag7 - HUISF_2003_5LL_Ag7)#*1000 TUGSF_2003_5LL_Ag7 <- (HU_2003_5LL_Ag7 + HG_2003_5LL_Ag7 + HS_2003_5LL_Ag7 + HF_2003_5LL_Ag7 - HUG_2003_5LL_Ag7 - HUS_2003_5LL_Ag7 - HUF_2003_5LL_Ag7 - HGS_2003_5LL_Ag7 - HGF_2003_5LL_Ag7 - HSF_2003_5LL_Ag7 + HUGS_2003_5LL_Ag7 + HUGF_2003_5LL_Ag7 + HUSF_2003_5LL_Ag7 + HGSF_2003_5LL_Ag7 - HUGSF_2003_5LL_Ag7)#*1000 TIGSF_2003_5LL_Ag7 <- (HI_2003_5LL_Ag7 + HG_2003_5LL_Ag7 + HS_2003_5LL_Ag7 + HF_2003_5LL_Ag7 - HIG_2003_5LL_Ag7 - HIS_2003_5LL_Ag7 - HIF_2003_5LL_Ag7 - HGS_2003_5LL_Ag7 - HGF_2003_5LL_Ag7 - HSF_2003_5LL_Ag7 + HIGS_2003_5LL_Ag7 + HIGF_2003_5LL_Ag7 + HISF_2003_5LL_Ag7 + HGSF_2003_5LL_Ag7 - HIGSF_2003_5LL_Ag7)#*1000 #QUINTUPLE TUIGSF_2003_5LL_Ag7 <- (HU_2003_5LL_Ag7 + HI_2003_5LL_Ag7 + HG_2003_5LL_Ag7 + HS_2003_5LL_Ag7 + HF_2003_5LL_Ag7 - HUI_2003_5LL_Ag7 - HUG_2003_5LL_Ag7 - HUS_2003_5LL_Ag7 - HUF_2003_5LL_Ag7 - HIG_2003_5LL_Ag7 - HIS_2003_5LL_Ag7 - HIF_2003_5LL_Ag7 - HGS_2003_5LL_Ag7- HGF_2003_5LL_Ag7 - HSF_2003_5LL_Ag7 + HUIG_2003_5LL_Ag7 + HUIS_2003_5LL_Ag7 + HUIF_2003_5LL_Ag7 + HUGS_2003_5LL_Ag7 + HUGF_2003_5LL_Ag7 + HUSF_2003_5LL_Ag7 + HIGS_2003_5LL_Ag7 + HIGF_2003_5LL_Ag7 + HISF_2003_5LL_Ag7 + HGSF_2003_5LL_Ag7 - HUIGS_2003_5LL_Ag7 - HUIGF_2003_5LL_Ag7 - HUISF_2003_5LL_Ag7 - HUGSF_2003_5LL_Ag7 - HIGSF_2003_5LL_Ag7 + HUIGSF_2003_5LL_Ag7)#*1000 TH4_2003_5LL_Ag7 <- do.call("cbind", list(HU_2003_5LL_Ag7, HI_2003_5LL_Ag7, HG_2003_5LL_Ag7, HS_2003_5LL_Ag7, HF_2003_5LL_Ag7, HUI_2003_5LL_Ag7, HUG_2003_5LL_Ag7, HUS_2003_5LL_Ag7, HUF_2003_5LL_Ag7, HIG_2003_5LL_Ag7, HIS_2003_5LL_Ag7, HIF_2003_5LL_Ag7, HGS_2003_5LL_Ag7, HGF_2003_5LL_Ag7, HSF_2003_5LL_Ag7, HUIG_2003_5LL_Ag7, HUIS_2003_5LL_Ag7, HUIF_2003_5LL_Ag7, HUGS_2003_5LL_Ag7, HUGF_2003_5LL_Ag7, HUSF_2003_5LL_Ag7, HIGS_2003_5LL_Ag7, HIGF_2003_5LL_Ag7, HISF_2003_5LL_Ag7, HGSF_2003_5LL_Ag7, HUIGS_2003_5LL_Ag7, HUIGF_2003_5LL_Ag7, HUISF_2003_5LL_Ag7, HUGSF_2003_5LL_Ag7, HIGSF_2003_5LL_Ag7, HUIGSF_2003_5LL_Ag7, TUI_2003_5LL_Ag7, TUG_2003_5LL_Ag7, TUS_2003_5LL_Ag7, TUF_2003_5LL_Ag7, TIG_2003_5LL_Ag7, TIS_2003_5LL_Ag7, TIF_2003_5LL_Ag7, TGS_2003_5LL_Ag7, TGF_2003_5LL_Ag7, TSF_2003_5LL_Ag7, TUIG_2003_5LL_Ag7, TUIS_2003_5LL_Ag7, TUIF_2003_5LL_Ag7, TUGS_2003_5LL_Ag7, TUGF_2003_5LL_Ag7, TUSF_2003_5LL_Ag7, TIGS_2003_5LL_Ag7, TIGF_2003_5LL_Ag7, TISF_2003_5LL_Ag7, TGSF_2003_5LL_Ag7, TUIGS_2003_5LL_Ag7, TUIGF_2003_5LL_Ag7, TUISF_2003_5LL_Ag7, TUGSF_2003_5LL_Ag7, TIGSF_2003_5LL_Ag7, TUIGSF_2003_5LL_Ag7, denom_2003_5_Ag7)) colnames(TH4_2003_5LL_Ag7) <- c("HU_2003_5LL_Ag7", "HI_2003_5LL_Ag7", "HG_2003_5LL_Ag7", "HS_2003_5LL_Ag7", "HF_2003_5LL_Ag7", "HUI_2003_5LL_Ag7", "HUG_2003_5LL_Ag7", "HUS_2003_5LL_Ag7", "HUF_2003_5LL_Ag7", "HIG_2003_5LL_Ag7", "HIS_2003_5LL_Ag7", "HIF_2003_5LL_Ag7", "HGS_2003_5LL_Ag7", "HGF_2003_5LL_Ag7", "HSF_2003_5LL_Ag7", "HUIG_2003_5LL_Ag7", "HUIS_2003_5LL_Ag7", "HUIF_2003_5LL_Ag7", "HUGS_2003_5LL_Ag7", "HUGF_2003_5LL_Ag7", "HUSF_2003_5LL_Ag7", "HIGS_2003_5LL_Ag7", "HIGF_2003_5LL_Ag7", "HISF_2003_5LL_Ag7", "HGSF_2003_5LL_Ag7", "HUIGS_2003_5LL_Ag7", "HUIGF_2003_5LL_Ag7", "HUISF_2003_5LL_Ag7", "HUGSF_2003_5LL_Ag7", "HIGSF_2003_5LL_Ag7", "HUIGSF_2003_5LL_Ag7", "TUI_2003_5LL_Ag7", "TUG_2003_5LL_Ag7", "TUS_2003_5LL_Ag7", "TUF_2003_5LL_Ag7", "TIG_2003_5LL_Ag7", "TIS_2003_5LL_Ag7", "TIF_2003_5LL_Ag7", "TGS_2003_5LL_Ag7", "TGF_2003_5LL_Ag7", "TSF_2003_5LL_Ag7", "TUIG_2003_5LL_Ag7", "TUIS_2003_5LL_Ag7", "TUIF_2003_5LL_Ag7", "TUGS_2003_5LL_Ag7", "TUGF_2003_5LL_Ag7", "TUSF_2003_5LL_Ag7", "TIGS_2003_5LL_Ag7", "TIGF_2003_5LL_Ag7", "TISF_2003_5LL_Ag7", "TGSF_2003_5LL_Ag7", "TUIGS_2003_5LL_Ag7", "TUIGF_2003_5LL_Ag7", "TUISF_2003_5LL_Ag7", "TUGSF_2003_5LL_Ag7", "TIGSF_2003_5LL_Ag7", "TUIGSF_2003_5LL_Ag7", "N") write.csv(TH4_2003_5LL_Ag7, "TH4_2003_5LL_Ag7.csv")