在SCI文章中交互效应表格通常是表五能为文章锦上添花增加文章的信服力增加结果的可信程度还能进行数据挖掘。什么是亚组通常就是特殊类 型人群比如男女种族等就是说你的数据放入特殊人群中结果还可靠吗如果在各个特殊人群中你的结果很稳定说明你的结论很可靠。如果亚组的结论和你的数据数据结论相反你可以拿来做个新论题。还可以比较不同亚组之间有无区别比如做了心脏支架和没做支架的区别可以发现很多新思路易于数据挖掘。其实交互作用就是一个变量通过对另一个变量的影响导致结局变量改变,最常见的就是我们药物之间的交互作用比如β受体阻滞剂对血糖有影响同时和降糖药使用可能会加强或者减弱降糖药的降糖效果这就是交互作用。既往我已经介绍过多篇交互效应文章比如R语言logistic回归交互项交互作用的可视化分析2这个文章介绍的是随着年龄增大流产2次以上的比没有流产的女性不孕症的概率明显增大因此流产和年龄对不孕症的结局存在交互影响。但是对于连续变量这样的图形不好表示交互效应根据粉丝要求编写了ggsci3d函数能够一步到位生成交互效应数据并绘图下面我来演示一下先导入数据和R包library(foreign)library(mgcv)library(rms)library(ggscitable)library(scitable)be-read.spss(E:/r/test/ozone.sav,use.value.labelsF,to.data.frameT)names(be)数据中有七个变量ozon每日臭氧水平为结局变量Inversion base heightibh反转基准高度Pressure gradient (mm Hg) 压力梯度mm HgVisibility (miles) 能见度英里Temperature (degrees F) 温度华氏度Day of the year日期vh我也不知道是什么反正就是一参数这里所有的变量都是连续的假设我想知道ibh和dpg对结局臭氧ozon的影响out-ggsci3d(xc(ibh,dpg),yozon,databe,usernameusername,tokentoken)这里m1和m2表示ibh和dpg预测值就是臭氧就得到一个随着ibh和dpg变化臭氧改变的3D图可以看到并不是ibh和dpg最大就是臭氧浓度最高而是一个曲线的变化有一个峰值。还可以进一步假如协变量的影响假设这里协变量是temp,“vis”out-ggsci3d(xc(ibh,dpg),yozon,covc(temp,vis),databe,usernameusername,tokentoken)可以自定义颜色有topo, “heat”, “cm”, “terrain”, “gray” ,bw这些颜色参数可以选择out-ggsci3d(xc(ibh,dpg),yozon,covc(temp,vis),colorcm,databe,usernameusername,tokentoken)咱们还可以增加区域网格数这样可以更加详细的显示曲线波动out-ggsci3d(xc(ibh,dpg),yozon,covc(temp,vis),colorcm,n.grid60,databe,usernameusername,tokentoken)给x轴y轴z轴都打上标签out-ggsci3d(xc(ibh,dpg),yozon,covc(temp,vis),colorcm,n.grid60,databe,xlab高度,ylab压力,main3D图,zlab臭氧,usernameusername,tokentoken)还可以更换其他计算3D图的算法out-ggsci3d(xc(ibh,dpg),yozon,databe,types,usernameusername,tokentoken)还可以使用等高线来表示out-ggsci3d(xc(ibh,dpg),yozon,databe,plot.typecontour,usernameusername,tokentoken)这个等高线画得明显不如咱们的ggconinteraction函数绘制效果如果想绘制等高线的可以看文章R语言基于ggscitable包复现一篇3.5分的文章的连续变量交互效应交互作用的可视化图ggsci3d函数还生成有绘图数据你可以借助绘图数据使用其他R包绘图也是可以的先把数据提取出来library(mgcv)library(rgl)library(plotly)x-out[[plotdt]][[x]]y-out[[plotdt]][[y]]z-out[[plotdt]][[z]]使用plot_ly函数绘图 plot_ly(xx,yy,zz,typesurface)%%layout(scenelist(xaxislist(titlex1),yaxislist(titlex2),zaxislist(titlefit)))这个图是能旋转的有需要旋转的可以使用这个函数绘图使用persp3d函数绘图persp3d(xx,yy,zz,xlabx1,ylabx2,zlabfit)文字表达有限下面视频更加精彩