我正在寻找R中的指标函数,即如果向量中元素的值大于0则返回1的函数,如果向量中元素的值小于0则返回0.
我需要在向量中的所有元素上使用此函数,返回一个只有零和一的新向量.
谢谢.
我有一个大小为TxR的矩阵,我正在寻找一个命令来做行的产品(返回产品的1 x R向量).之后我想总结一下列,即总结R项.
在Matlab中,这将完成类似于这个总和(prod(A,1),2),但我不知道R中的代码.
我希望它有意义.
谢谢
我有一个矩阵,在R中包含三列.矩阵看起来像这样:
A=matrix(c(1,2,3,4,0.5,1,7,1.2,3,4,2,1),nrow=4, ncol=3)
我想创建一个基于A的矩阵,在A的每一行中,该行返回1表示该行中的最高值,否则返回0.所以在上面的具体情况中,我需要一个如下所示的矩阵:
B=matrix(c(0,0,0,1,0,0,1,0,1,1,0,0),nrow=4,ncol=3)
我试图搜索论坛,但找不到合适的答案.
谢谢.
我在R中使用了一些代码.我试图在同一个图中创建6个不同的直方图.它工作正常,但我需要在6个直方图的每一个中放置1条垂直线.我正在使用的代码可能如下所示:
require(ggplot2)
require(reshape2)
require(gdata)
MC_beta=rbind(rnorm(1000,-2,0.1),rnorm(1000,-1,0.1),rnorm(1000,0,0.1),rnorm(1000,0.5,0.1),rnorm(1000,1,0.1),rnorm(1000,2,0.1))
df <- data.frame(MC_beta[1,], MC_beta[2,], MC_beta[3,], MC_beta[4,],MC_beta[5,],MC_beta[6,])
names(df)[1:6]<-c("1", "2", "3", "4","5","6")
df2 = melt(df)
z=c(-2,-1,0,0.5,1,2)
ggplot(df2, aes(x=value, fill = variable)) + geom_vline(xintercept = z, colour="black") +geom_histogram(binwidth=0.03,colour = "black") + scale_fill_manual(name = "",values = c('red','blue',"red","blue","red","blue")) +
facet_wrap(~variable,nrow=6, ncol=1) + scale_x_continuous(breaks=seq(-2.5,2.5,0.5)) + guides(fill=FALSE) +
theme_bw() + theme(strip.background = element_blank(),axis.text=element_text(size=14.5),strip.text.x = element_text(size = 14.5)) + stat_function(fun = dnorm)
问题在于声明geom_vline(xintercept = z, colour = "black").显然,不是在每个直方图中放置一条垂直线,而是在每个直方图中放置所有6条线.所以相反,我希望z中的第一个元素在第一个直方图中形成一条垂直线,z中的第二个元素在第二个直方图中形成一条垂直线,所以第四个.
谢谢
我正在尝试最大化一个函数,在我的RI代码中使用以下函数来确保,我得到的值为A,介于-1和1之间:
A=(exp(x)-1)/(exp(x)+1)
问题是,当我例如在x = 3000或某个其他极值处评估函数时,它返回NaN而不是-1,它应该是.
有没有办法写这个以确保,这不会发生？
谢谢.
我想创建一个离散变量,它将变量x除以1000的间隔.我所拥有的是以下内容:
DATA have;
INPUT x;
DATALINES;
1200
3200
5300
49
6500
;
RUN;
我想要的是这样的
data want;
input x y $5-14;
DATALINES;
1200 1000-2000
3200 3000-4000
5300 5000-6000
0049 0000-1000
6500 6000-7000
;
RUN;
我当然可以使用if else语句编写数据步骤来定义间隔,但是我有100个这样的间隔,所以这将是非常耗时的.
谢谢.