小编ECI*_*CII的帖子

我想将数据框的最后一列移到开头(作为第一列).我怎么能在R？

我的data.frame有大约一千列来改变订单.我只想选择一列并"将其移至开头".

我有4个绘图的汇编par(mfrow=c(2,2)).我想为上面的2个图绘制一个共同的标题,并在下面的2个面板中绘制一个共同的标题,这些标题位于2个左右图之间.

这可能吗？

可以在R中的2行或更多行中分割plot()的主标题吗？

此外,我有一个标题,其中包含许多粘贴在一起的论证

plot(1, main=paste("X:",1," ","Y:", 2," ","Z:",3))

如何在主要的Y和Z之间插入换行符？

谢谢EC

似乎虽然grep有一个反转参数,grepl但没有.

我想使用2个过滤器进行子集化

data$ID[grepl("xyx", data$ID) & data$age>60]

如何为年龄> 60而ID和不包含 "xyx"的ID进行子集化？我做的是

data$ID[abs(grepl("xyx", data.frame$ID)-1) & data$age>60]

这看起来很有效,但看起来很糟糕且不直观.有更好的解决方案/参数吗？

谢谢

我有以下模型

x  <- rep(seq(0, 100, by=1), 10)
y  <- 15 + 2*rnorm(1010, 10, 4)*x + rnorm(1010, 20, 100)
id <- NULL
for(i in 1:10){ id <- c(id, rep(i,101)) }
dtfr <- data.frame(x=x,y=y, id=id)
library(nlme)
with(dtfr, summary(     lme(y~x, random=~1+x|id, na.action=na.omit)))
model.mx <- with(dtfr, (lme(y~x, random=~1+x|id, na.action=na.omit)))
pd       <- predict( model.mx, newdata=data.frame(x=0:100), level=0)
with(dtfr, plot(x, y))
lines(0:100, predict(model.mx, newdata=data.frame(x=0:100), level=0), col="darkred", lwd=7)

用predict,level=0我可以绘制平均人口反应.如何从nlme对象中提取和绘制整个群体的95%置信区间/预测带？

我有

library(ISLR)
attach(Wage)

# Polynomial Regression and Step Functions

fit=lm(wage~poly(age,4),data=Wage)
coef(summary(fit))

fit2=lm(wage~poly(age,4,raw=T),data=Wage)
coef(summary(fit2))

plot(age, wage)
lines(20:350, predict(fit, newdata = data.frame(age=20:350)), lwd=3, col="darkred")
lines(20:350, predict(fit2, newdata = data.frame(age=20:350)), lwd=3, col="darkred")

预测线似乎是相同的,但为什么系数如此不同？你如何解读为此他们raw=T和raw=F.

我看到产生的系数与poly(...,raw=T)那些系数相匹配~age+I(age^2)+I(age^3)+I(age^4).

如果我想使用系数来"手动"(不使用predict()函数)获得预测,那么我应该注意什么？我该如何解释正交多项式的系数poly().

我在R中使用lme4来适应混合模型

lmer(value~status+(1|experiment)))

值是连续的,状态(N/D/R)和实验是因素,我得到

Linear mixed model fit by REML 
Formula: value ~ status + (1 | experiment) 
  AIC   BIC logLik deviance REMLdev
 29.1 46.98 -9.548    5.911    19.1
Random effects:
 Groups     Name        Variance Std.Dev.
 experiment (Intercept) 0.065526 0.25598 
 Residual               0.053029 0.23028 
Number of obs: 264, groups: experiment, 10

Fixed effects:
            Estimate Std. Error t value
(Intercept)  2.78004    0.08448   32.91
statusD      0.20493    0.03389    6.05
statusR      0.88690    0.03583   24.76

Correlation of Fixed Effects:
        (Intr) statsD
statusD -0.204       
statusR -0.193  0.476

我想以图形方式表示固定效应评估.但是这些对象似乎没有绘图功能.有什么方法可以用图形方式描绘固定效果吗？

我目前正在使用R进行所有科学计算和绘图,但我想探索Python.我一直使用RStudio作为R的IDE,作为IDE,它满足了我对科学计算,数字运算,数据分析和可视化的100%需求.Python有类似的东西吗？基本上是语法高亮,代码完成,智能缩进,直接从源编辑器执行代码,在IDE和版本控制(git)中绘图.

我觉得挑战是,建立在R.美观的三维曲面我熟悉的解决方案(persp,image,wireframe,lattice,rgl和几个在SO其他问题的其他解决方案),但是结果不是很好.

是否可以像在MATLAB中一样在R中创建3D曲面图？

Here is the MATLAB code
    % Create a grid of x and y points
    points = linspace(-2, 0, 20);
    [X, Y] = meshgrid(points, -points);

    % Define the function Z = f(X,Y)
    Z = 2./exp((X-.5).^2+Y.^2)-2./exp((X+.5).^2+Y.^2);

    % "phong" lighting is good for curved, interpolated surfaces. "gouraud"
    % is also good for curved surfaces
    surf(X, Y, Z); view(30, 30);
    shading interp;
    light;
    lighting phong;
    title('lighting phong', 'FontName', 'Courier', 'FontSize', 14);

情节现代,色彩鲜艳,美观,代码语法非常易读. …

除了base和之外,如何删除所有已安装的软件包recommended？