小编Sve*_*ein的帖子

是否可以使用该omit选项省略观星者中的互动？通常我会将变量名写入,omit=c('varname')但在交互的情况下,我不知道写什么.有什么提示吗？你如何在其他包中解决这个问题xtable？

\documentclass{article}   
\begin{document}   
%Load dataset and run regression    

 << lm, echo=FALSE >>=   
 load('dataset.RData')  
 library(stargazer)  
 lm1 <- lm(y~ x + factor(v)*z ,data=dataset)  
 @


 << table_texstyle, echo=FALSE, comment=NA, results='asis' >>=  
 stargazer(lm1 ,omit=c('???'), omit.labels=c('Omitted interactions'),
      omit.yes.no = c("Yes", "No")))  
  @  
   \end{document}

我有一个数据框z,我想根据两个旧列的值创建新列z.以下是过程:

>z<-cbind(x=1:10,y=11:20,t=21:30)
> z<-as.data.frame(z)
>z
    x  y  t
1   1 11 21
2   2 12 22
3   3 13 23
4   4 14 24
5   5 15 25
6   6 16 26
7   7 17 27
8   8 18 28
9   9 19 29
10 10 20 30

#生成列q等于列t时间4 的值,如果x=3和其他值x,则等于列的值t.

for (i in 1:nrow(z)){
  z$q[i]=if (z$x[i]==4) 4*z$t[i] else z$t[i]}

但是,我的问题是我想应用多个条件:

例如,我想得到这样的东西:

(If x=2, q=t*2; x=4, q=t*4; …

我有一个以下面的方式创建的数据框.

library(ggplot2)

x <- data.frame(letters[1:10],abs(rnorm(10)),abs(rnorm(10)),type="x")
y <- data.frame(letters[1:10],abs(rnorm(10)),abs(rnorm(10)),type="y")
 # in reality the number of row could be larger than 10 for each x and y

all <- rbind(x,y)
colnames(all) <- c("name","val1","val2","type")

我想要做的是创建一个看起来大致如下的刻面ggplot:

因此,上面的每个方面都是以下相关图:

# Top left facet
subset(all,type=="x")$val1 
subset(all,type=="y")$val1

# Top right facet
subset(all,type=="x")$val1 
subset(all,type=="y")$val2

# ...etc..

但我坚持使用以下代码:

p <- ggplot(all, aes(val1, val2))+ geom_smooth(method = "lm")  + geom_point() +
facet_grid(type ~ ) 
# Calculate correlation for each group
cors <- ddply(all, c(type ~ ), summarise, cor = round(cor(val1, val2), 2)) …

我试图计算一个r在T日期范围内的回报时间序列的滚动总和.但是,在我计算滚动总和的每个日期时,我想要计算滚动总和w中每个数字的权重.

该公式适用于每个日期t:

[Sum from i=1 to m](w(i)*r(t-i-1))

让我们看一个非常简单的例子.我有一个T = 6返回的返回系列r.对于每个日期,t我想计算过去两个日期(m = 2)的滚动总和.我还想把第一次观察的重量增加到第二次.

r <- c(100,110,100,110,100,110)                 
w <- c(1,0.5)

我知道我可以使用过滤功能轻松完成滚动总和:

filter(r, rep(1, 2))

但是我无法将权重因子包含在滚动总和中.以下行给出了错误的结果c(155, 155, 155, 155, 155, NA)

filter(r*w, rep(1, 2))

在哪里我真的想得到结果 c(155, 160, 155, 160, 155, NA)

任何帮助表示赞赏.

我创建了各种参考类来适应一些 arima、garch 过程，并希望在并行计算中使用它们 parSapply

我先做了一些导出

cl <- makeCluster(mc <- getOption("cl.cores", 20))
        clusterExport(cl, c("merge.xts", "index", "coredata", "xts", "lag.xts", "zoo", "LearnerPredict", "arima", "generic_learner", "arma_simple", "logwarn"))
        clusterEvalQ(cl, "arma_simple")
        clusterEvalQ(cl, "generic_learner")


generic_learner <- setRefClass(

        Class = "generic_learner",

        fields = list(

            params = "list"

        ),

        methods = list(

            fitModel = function() {cat("overload function with fitting function \n")},

            fcastModel = function() {cat("overload function with forecast function \n")},

            fmt_params = function() {cat("overload function with formatted parameters \n")},

            fmt_class = function() {cat("overload class\n")},

            fmt_ref = function() {paste(.self$fmt_class(), .self$fmt_params(), …

我有一个Rhtml文件,我从中获取一个R文件.在这个R文件中,我正在做一些情节.

p=ggplot(data)
p+geom_line()

现在,我可以在另一个之后产生一个情节,knit(".Rhtml")然后在做的 时候我接着另一个.

但我想把这些数字并排.(数据的数量因报告而异).

有没有办法在Rhtml文件中设置一个选项,以便数字并排排列(例如两列,三列或四列).

所以,实际上它会像是一样par(mfrow).

我想创建一个矩阵n通过k与k使用循环MVN协变量.相当简单,但到目前为止还没有...这是我的代码:

n=1000
k=5
p=100
mu=0
sigma=1
x=matrix(data=NA, nrow=n, ncol=k)


for (i in 1:k){
        x [[i]]= mvrnorm(n,mu,sigma)
       }

少了什么东西？

如何标记此轴标题的单位.该""代表我需要上标部分:Photosynthetically available radiation (µE m"-2"d"-1").

到目前为止,我使用了公式并且没有运气:

plot(PAR~SST,data=brazilw, pch=15,col="red", main ="Fig. 1. Relationship between photosynthically available radiation\n and sea surface temperature",
ylab=expression("Photosynthetically available radiation (µE m"^-2~d^-1))

我有一个大型数据集,其中所有列标题都是单独的IDS,每个长度为8个字符.我想将这些单独的ID分成2行,其中第一行ID包含前7个字符,第二行仅包含最后一个字符.

当前数据集:

ID1:    Indiv01A    Indiv01B    Indiv02A    Indiv02B    Speci03A    Speci03B

预期数据集:

ID1:    Indiv01 Indiv01 Indiv02 Indiv02 Speci03 Speci03  
ID2:    A   B   A   B   A   B

我查看了有关拆分数据的其他帖子,但它们似乎都有一种独特的方式来分隔列名称(即:有一个逗号分隔2个组件或句点).

这是我认为最好的代码,但我无法弄清楚如何编写"7个字符"作为分割点,而不是逗号:

sapply(strsplit(as.character(d$ID), ",")

任何帮助,将不胜感激.

苦苦挣扎几个小时才能获得这场比赛并在R中取代gsub工作但仍然没有成功.我正在尝试匹配"Reason:"字符串中的模式,并在此模式之后触发所有内容,直到第一次出现一个点(.)例如:

Offer Disposition. MSISDN: 7183067962. Offer: . Disposition: DECLINED. Reason: Not interested. ChannelID: CARE.

会回来的 "Not interested"