是否有一个R函数来进行多重关联，而不是重写代码块
目前，我必须以多种方式对数据进行子集，以提取数据，我需要进行关联，并确定这些组的相同金属之间是否存在关系是否有一个R函数来进行多重关联，而不是重写代码块,r,correlation,R,Correlation,目前，我必须以多种方式对数据进行子集，以提取数据，我需要进行关联，并确定这些组的相同金属之间是否存在关系 Group1 <- subset(Data_Set, subset = Data_Set$Sample == "1") Group1A <- subset(Group1, subset = Group1$Sample_Type == "A") GroupX <- subset(Data_Set, subset = Data_Set$Sample == "X") GroupX
Group1 <- subset(Data_Set, subset = Data_Set$Sample == "1")
Group1A <- subset(Group1, subset = Group1$Sample_Type == "A")
GroupX <- subset(Data_Set, subset = Data_Set$Sample == "X")
GroupX<- subset(GroupX, Sample_ID %in% Group1A$Sample_ID )
cor.test(Group1A$Pb,GroupX$Pb, method = "kendall")
structure(list(Sample = c("2", "2", "2", "2", "X", "2", "2",
"2", "2", "2", "2", "X", "2", "2", "5", "5", "5", "5", "5", "X",
"5", "5", "3", "3", "X", "3", "3", "X", "4", "4"), Sample_ID = c("DC001",
"DC001", "DC001", "DC001", "DC001", "DC001", "DC001", "DC002",
"DC002", "DC002", "DC002", "DC002", "DC002", "DC002", "DC003",
"DC003", "DC003", "DC003", "DC003", "DC003", "DC003", "DC003",
"DC004", "DC004", "DC004", "DC005", "DC005", "DC005", "DC006",
"DC006"), Sample_Type = c("A", "D", "E", "F", "X", "I", "J",
"A", "D", "E", "F", "X", "I", "J", "A", "B", "D", "E", "F", "X",
"I", "J", "C", "F", "X", "C", "F", "X", "A", "D"), Co = c(0,
0.204473214269861, 0, 0.50977856054987, 0.262230521160956, 0,
0, 0, 0, 0, 0.465855303428853, 0.229502158969648, 0.214970121592712,
0.588126362402572, 0, 0.0906122639531158, 0.229838105464066,
0, 0.240533898070871, 4.77802122014029, 0.47537095149254, 0.384495379166814,
0.00135414270258444, 0.458235177876183, 0.412977043885698, 0.187579567424379,
0.317854941692133, 0.0271598068567071, 0, 0.293328743450483),
Ni = c(2.32894078024542, 0, 2.75976812547636, 2.35251746719724,
0.351631195258774, 1.25476391714642, 0.0586626807902249,
0, 2.31716731851309, 0, 4.03426936736104, 0.414520597983989,
2.69897385721456, 0.781651988488391, 1.48260693680732, 1.59083944326126,
0.944038748319438, 3.06889126279262, 1.69552165261712, 0.849220149877567,
1.75387912556474, 0, 0.333762199305291, 1.66187141150986,
0.735834552887327, 3.72419677755011, 1.27862769479216, 0.264762516047524,
1.84288031704096, 1.8828793053893), Cu = c(16.6696573471153,
21.377014252538, 16.4581203986139, 6.49438237470201, 1.57054125960644,
5.67180974109468, 23.5835333332964, 38.6483288663375, 15.2589198442198,
21.9746392829346, 7.09307693625389, 0.967127488045321, 6.32542891436958,
16.1173426649179, 11.2222721930992, 8.42093833910001, 11.1332246071585,
16.7442343774396, 10.8140656299147, 14.2632807636599, 5.35502290473828,
7.29141216675894, 2.53789491234011, 16.5791995430022, 1.00648647764661,
26.6313784234462, 0.0413060789264422, 0.656674377606213,
3.98095036332964, 6.17760205144632), Zn = c(76.5281110975817,
2652.50181007495, 1007.00556337852, 206.99812727191, 640.15733114957,
484.221162531697, 3718.61286231799, 131.574098527507, 9826.49966864988,
1827.75831773692, 557.015412652748, 850.519284594127, 955.085171501707,
3039.23169926716, 117.947177178762, 65.7886442827721, 78.1092625035093,
253.691311074245, 980.544294923672, 506.400193234096, 1110.92409209043,
902.659801267825, 284.143460051779, 991.762202132739, 899.71040333897,
1686.99915717559, 27.0835877755038, 956.364728487396, 142.167067778216,
1012.61495002819), As = c(0, 0, 1.91185052013389, 1.32808264279786,
0.141039242323703, 1.74872331719823, 0.1065340816859, 0.812367854870543,
0, 0.797230094696634, 2.38925992872935, 0.305621793073037,
0.664951374730799, 0, 0, 0, 2.52051964809224, 0, 0, 0.392178178336116,
0, 3.08334159340895, 2.32108729394528, 1.62081021652742,
0.171200134084414, 6.19125023716284, 4.43213876523911, 0.289386770990403,
0.313331113399545, 6.41607755268465), Cd = c(8.22465741493669,
22.6126042664945, 34.0150873273517, 13.5844058876617, 5.22665850051452,
24.0465414683255, 109.478598702669, 15.1992477278811, 169.517190223851,
75.2983940524065, 34.5230481628261, 3.75297525105592, 45.6178498733986,
247.435132822196, 2.10793502840313, 1.47647473271431, 0.0848090794945706,
2.98717760781629, 3.13384011407655, 5.31936421369202, 3.73593799828465,
5.36310372449921, 0.298562637256625, 1.82673831232711, 3.78462211601718,
8.0628550389363, 0.138799690323038, 1.32275598609847, 0.285061500560821,
0.635235209786838), Pb = c(0.922803462498185, 5.13959353157866,
1.9525414480789, 0, 2.5902978681043, 1.21865949505257, 7.09067896476338,
0, 3.89524247237658, 0.354938950934777, 2.64634863087263,
0.356658949506862, 1.25701617111933, 4.18799241835111, 0,
0.807369345092201, 0.0263264119388502, 0, 3.32333444396018,
76.7555925603143, 0.613522400825461, 0, 1.72315815094652,
3.21414903849599, 1.03802696495681, 1.73176109371547, 0.72736174943572,
0.23309888503164, 12.8688959655249, 33.2486209089115)), row.names = c(NA,
-30L), class = c("tbl_df", "tbl", "data.frame"))
这里有一个解决方案，让你得到你想要的。它看起来有点长，让人望而生畏，但我试着让它更容易理解，更容易扩展，更容易防止遗漏的观察结果，等等。。我还有一个蛮力解决方案，可以**所有*相关性，然后删除所有不需要的相关性 图书馆弹琴 图书馆长 图书馆咕噜声 图书馆扫帚 组成的数据集与您的数据集相似，但缺少行和NAs 设定2020年种子 数据集% 突变组=X%>% 重命名ECO_X=Co， Ni_X=Ni， Cu_X=Cu， Zn_X=Zn， As_X=As， Cd_X=Cd， Pb_X=Pb%>% 选择组、-Sample、-Sample\u类型 没有XX的df AllNotX% 过滤器组！=XX 按组列出DF的列表 ListofGroupDFs% 分开。$组 格列美脲
组Dfs列表您可以发布示例数据吗？请使用dputData_集的输出编辑问题。或者，如果dputheadData_集的输出太大，则为30。@RuiBarradas done！请注意：在您的示例数据中，没有Group1等。试图复制您的示例数据的人需要弄清楚您的数据是什么样子，并相应地调整您的代码，这并不理想。如果能以问题中的格式化方式设置headData_的输出也会很好，这样可以更容易地看到您试图实现的目标。感谢您让我知道@coffeinjunky我已经编辑了它，所以它应该可以工作。但如果没有，请告诉我。我一直在看数据，看不到获得相关性的方法。我可以很容易地找到一种方法来比较任何一种金属的平均水平或中间水平。无法从您的样本数据判断发生了什么。
structure(list(Sample = c("2", "2", "2", "2", "X", "2", "2",
"2", "2", "2", "2", "X", "2", "2", "5", "5", "5", "5", "5", "X",
"5", "5", "3", "3", "X", "3", "3", "X", "4", "4"), Sample_ID = c("DC001",
"DC001", "DC001", "DC001", "DC001", "DC001", "DC001", "DC002",
"DC002", "DC002", "DC002", "DC002", "DC002", "DC002", "DC003",
"DC003", "DC003", "DC003", "DC003", "DC003", "DC003", "DC003",
"DC004", "DC004", "DC004", "DC005", "DC005", "DC005", "DC006",
"DC006"), Sample_Type = c("A", "D", "E", "F", "X", "I", "J",
"A", "D", "E", "F", "X", "I", "J", "A", "B", "D", "E", "F", "X",
"I", "J", "C", "F", "X", "C", "F", "X", "A", "D"), Co = c(0,
0.204473214269861, 0, 0.50977856054987, 0.262230521160956, 0,
0, 0, 0, 0, 0.465855303428853, 0.229502158969648, 0.214970121592712,
0.588126362402572, 0, 0.0906122639531158, 0.229838105464066,
0, 0.240533898070871, 4.77802122014029, 0.47537095149254, 0.384495379166814,
0.00135414270258444, 0.458235177876183, 0.412977043885698, 0.187579567424379,
0.317854941692133, 0.0271598068567071, 0, 0.293328743450483),
Ni = c(2.32894078024542, 0, 2.75976812547636, 2.35251746719724,
0.351631195258774, 1.25476391714642, 0.0586626807902249,
0, 2.31716731851309, 0, 4.03426936736104, 0.414520597983989,
2.69897385721456, 0.781651988488391, 1.48260693680732, 1.59083944326126,
0.944038748319438, 3.06889126279262, 1.69552165261712, 0.849220149877567,
1.75387912556474, 0, 0.333762199305291, 1.66187141150986,
0.735834552887327, 3.72419677755011, 1.27862769479216, 0.264762516047524,
1.84288031704096, 1.8828793053893), Cu = c(16.6696573471153,
21.377014252538, 16.4581203986139, 6.49438237470201, 1.57054125960644,
5.67180974109468, 23.5835333332964, 38.6483288663375, 15.2589198442198,
21.9746392829346, 7.09307693625389, 0.967127488045321, 6.32542891436958,
16.1173426649179, 11.2222721930992, 8.42093833910001, 11.1332246071585,
16.7442343774396, 10.8140656299147, 14.2632807636599, 5.35502290473828,
7.29141216675894, 2.53789491234011, 16.5791995430022, 1.00648647764661,
26.6313784234462, 0.0413060789264422, 0.656674377606213,
3.98095036332964, 6.17760205144632), Zn = c(76.5281110975817,
2652.50181007495, 1007.00556337852, 206.99812727191, 640.15733114957,
484.221162531697, 3718.61286231799, 131.574098527507, 9826.49966864988,
1827.75831773692, 557.015412652748, 850.519284594127, 955.085171501707,
3039.23169926716, 117.947177178762, 65.7886442827721, 78.1092625035093,
253.691311074245, 980.544294923672, 506.400193234096, 1110.92409209043,
902.659801267825, 284.143460051779, 991.762202132739, 899.71040333897,
1686.99915717559, 27.0835877755038, 956.364728487396, 142.167067778216,
1012.61495002819), As = c(0, 0, 1.91185052013389, 1.32808264279786,
0.141039242323703, 1.74872331719823, 0.1065340816859, 0.812367854870543,
0, 0.797230094696634, 2.38925992872935, 0.305621793073037,
0.664951374730799, 0, 0, 0, 2.52051964809224, 0, 0, 0.392178178336116,
0, 3.08334159340895, 2.32108729394528, 1.62081021652742,
0.171200134084414, 6.19125023716284, 4.43213876523911, 0.289386770990403,
0.313331113399545, 6.41607755268465), Cd = c(8.22465741493669,
22.6126042664945, 34.0150873273517, 13.5844058876617, 5.22665850051452,
24.0465414683255, 109.478598702669, 15.1992477278811, 169.517190223851,
75.2983940524065, 34.5230481628261, 3.75297525105592, 45.6178498733986,
247.435132822196, 2.10793502840313, 1.47647473271431, 0.0848090794945706,
2.98717760781629, 3.13384011407655, 5.31936421369202, 3.73593799828465,
5.36310372449921, 0.298562637256625, 1.82673831232711, 3.78462211601718,
8.0628550389363, 0.138799690323038, 1.32275598609847, 0.285061500560821,
0.635235209786838), Pb = c(0.922803462498185, 5.13959353157866,
1.9525414480789, 0, 2.5902978681043, 1.21865949505257, 7.09067896476338,
0, 3.89524247237658, 0.354938950934777, 2.64634863087263,
0.356658949506862, 1.25701617111933, 4.18799241835111, 0,
0.807369345092201, 0.0263264119388502, 0, 3.32333444396018,
76.7555925603143, 0.613522400825461, 0, 1.72315815094652,
3.21414903849599, 1.03802696495681, 1.73176109371547, 0.72736174943572,
0.23309888503164, 12.8688959655249, 33.2486209089115)), row.names = c(NA,
-30L), class = c("tbl_df", "tbl", "data.frame"))