我正在处理金融/经济数据,以防你想知道下面一些系数的大尺寸…我的一般问题与从R中的线性随机效应模型输出的参数系数的模拟有关.我是尝试使用来自R中相同模型的模型系数和方差 – 协方差(VCOV)矩阵生成β系数的随机样本.我…
标签:r
Shiny / ggvis对子集图数据的反应性
在看了GGVIS的演讲之后,我一直在努力创建我的第一个Shiny / ggvis应用程序.我的情节在R中有效,但当我尝试将其迁移到Shiny应用程序以便在网络上显示时,我什么也得不到.单选按钮显示,我可以告诉它似乎工作(…
R中的时间序列 – 使用不同的时间戳对齐数据
我有一组二元数据.因此,对话中的每个人都可以在5分钟(300秒)内自由切换2个任务之间的次数,并记录每个人在任务之间切换的时间. Participant A Participant B Time Task Time Ta…
在ggplot2中强制密度图的颜色
我有两个重叠密度函数,我创建如下: require(ggplot2) set.seed(2) a =rnorm(100) b = rnorm(100,.5,1.2) c = rnorm(100,.3,1.2) data …
在knitr中更改注释字体
在官方 knitr manual中,评论是绿色的,它们是正直的而不是斜体.但每当我使用knitr编译PDF时,我的评论都是紫色和斜体.我怎么能改变这个? 最佳答案 回答我自己的问题: 感谢mnel的小费.只需使用knit…
R使用data.table进行多个条件连接
我有一个大型数据集和一个查找表.我需要为数据集中的每一行返回查找中满足条件的行的最小值. 鉴于我的数据集的大小,我不愿意通过交叉连接来共同破解一个iffy解决方案,因为这将创建数百万条记录.我希望有人可以建议一个(理想情…
使用时间序列参数来解决R中的ODE
我试图使用deSolve在R中解决一个简单的ODE:dQ / dt = f(Q)*(P – E).整个事情是Q的时间序列.诀窍是P和E是固定的时间序列数据本身,因此diff eq仅在Q中有效. 用固定的时间步长迭代地解决…
R加速矢量化为方阵
有人能帮我加速一些代码: n = seq_len(ncol(mat)) # seq 1 to ncol(mat) sym.pr<-outer(n,n,Vectorize(function(a,b) { return…
R绘制二进制时间序列
我有一个数据集,我在时间轴上有一个二进制值. 例如: Date, Event January-29-2014, 1 January-29-2014, 0 January-29-2014, 1 January-29-201…
r – tm.package:findAssocs vs Cosine
我是新来的,我的问题是数学而不是编程本质,我希望得到关于我的方法是否有意义的第二意见. 我试图使用tm包中的函数findAssocs在我的语料库中找到单词之间的关联.尽管它似乎在通过该软件包提供的数据上表现得相当不错,例…
在R中更新zeroinfl模型
我有一个示例数据框如下. data.df = data.frame(Y=c(0,0,0,0,0,0,0,0,0,0,1,2,3),X1=c(3,5,3,2,5,6,3,5,1,3,1,7,8),X2=c(6,2,1,6,…
lm()或gls()的正系数
我想知道是否有人知道如何限制lm或gls中的某些系数为正? 例如,我想要lm(y~x1 x2 … xn)或gls(y~x1 x2 … xn),我希望x1和x2的系数是非负的. 最佳答案 一种选择是使用nls(…),它允许…