小编ami*_*mit的帖子

从R中的lm中提取标准化系数

我为这个愚蠢的问题道歉......但我似乎无法找到一个简单的解决方案

我想从拟合的线性模型中提取标准化系数(在R中)必须有一个简单的方法或函数来做到这一点.你能告诉我它是什么吗?

编辑(以下一些评论):我应该提供有关我的问题的更多上下文信息.我正在为一群心理学家教授一个介绍性的R工作室.对于他们来说,没有能够获得标准化系数的线性模型就好像你根本没有运行模型一样(好吧,这有点夸张,但你明白了).当我们做了一些回归时,这是他们的第一个问题,我(我的坏)我没有预料到(我不是心理学家).当然我可以自己编程,当然我可以找到适合我的软件包.但与此同时,我认为这是线性模型的一个基本和常见的必需功能,在现场,我认为应该有一个基本的功能,无需安装越来越多的包(这被认为是初学者的难点).所以我问(这也是向他们展示如何在需要时获得帮助的机会).

我为那些认为我问过一个愚蠢问题的人道歉,并且非常感谢那些花时间回答问题的人.

beta regression r standardized lm

35
推荐指数
2
解决办法
3万
查看次数

绘制围绕一组点的曲线

我在飞机上有一套点.它们被划分为子集.我想围绕属于同一子集的点绘制闭合曲线,以便属于子集的点将位于曲线内部,而不属于该子集的点将位于曲线内部.因此,简单的圆形或凸包可能不起作用.

对于初学者,假设我只想在一组点附近有一条平滑的曲线(不要求它排除其他点)

任何想法如何在R?

---后来补充---

我最终看到的是这里的图形精神:https://tex.stackexchange.com/questions/1175/drawing-a-hypergraph - 虽然上下文不是超图,而是给定的点集和那些分区.

curve r graphics2d

27
推荐指数
3
解决办法
5464
查看次数

使用循环在rmarkdown中生成文本部分

我需要生成一个由几个部分组成的报告,所有部分看起来都很相似,只有一些数据差异.部分的数量也取决于数据.我最终想拥有的是这样的:

```{r}
  section_names = c("A","B","C")
  section_data = c(13,14,16)
```

# some looping mechanism here with variable i

This is section `r section_names[i]`

This section's data is `r section_data[i]`

#more things go here for the section

#end of loop should go here
Run Code Online (Sandbox Code Playgroud)

结果应该是一个单独的html /文档,所有部分一个接一个.

你能指点我用循环生成这样一个Rmd文件的方法吗?

理想情况下,我希望看到像PHP中的东西:

<$php for(i=0;i<10;i++) { ?>
   ## some html template + code chunks here
<$php } ?>
Run Code Online (Sandbox Code Playgroud)

r knitr r-markdown

13
推荐指数
1
解决办法
6428
查看次数

R/shiny中的多选框 - 添加滚动条

我构建了一个R /闪亮的Web应用程序.我想要一个多选框(我使用checkboxGroupInput(),但我可以选择其他选项).但是,选项列表很长,我希望将它包含在一个相对较小的选项框中(一次显示5-6个选项),滚动条可以滚动整个选项列表.

有没有办法可以做到这一点?最小的例子:

ui.R

library(shiny)
choices = paste("A",1:30,sep="_")

shinyUI(pageWithSidebar(

# Application title
headerPanel("my title"),
sidebarPanel(   
  checkboxGroupInput("inp", "choose any of the following", choices)
),
mainPanel(
   tableOutput("result")

)
))
Run Code Online (Sandbox Code Playgroud)

server.R

library(shiny)
shinyServer(function(input, output) {
myInput <- reactive({
    input$inp
})
output$result <- renderTable({
x = myInput()
if(length(x)==0) {
x = "No Choice Made"
}
matrix(x,ncol=1)
})

})
Run Code Online (Sandbox Code Playgroud)

r shiny

10
推荐指数
1
解决办法
2198
查看次数

在R中使用stargazer而对乳胶一无所知

我即将为社会科学家教授R课程.他们可能对LaTeX一无所知,也不愿意听到它(R对他们来说足够复杂).然而,他们很可能喜欢观星包装所代表他们模特的那种桌子.

是否有任何包装器或其他简单的程序使他们能够使用stargazer直接创建表格(或其他格式)的pdf,然后他们可以作为图像插入到word文档中?

latex r ms-word stargazer

10
推荐指数
2
解决办法
9537
查看次数

连接/合并时间序列(在R中)

我必须使用xts/zoo对象.每个都有不同时间跨度的不同变量的度量.我想创建一个包含所有度量的单个时间序列,以及缺少日期/变量组合的NA.我怎么做?人为的例子:

library(xts)
x<-cbind(a=1:3,b=3:1)
rownames(x) = as.character(Sys.Date()-1:3)

y<-cbind(a=5:7,c=3:1)
rownames(y) = as.character(Sys.Date()-5:7)

xs=as.xts(x)
ys=as.xts(y)

#now what?

#desired outcome looks like:
            a  b   c
2013-03-10  7 NA   1
2013-03-11  6 NA   2
2013-03-12  5 NA   3
2013-03-14  3  1  NA
2013-03-15  2  2  NA
2013-03-16  1  3  NA

# regular merge looks like that (adding an a.1 variable)
merge(xs,ys)
            a  b a.1  c
2013-03-10 NA NA   7  1
2013-03-11 NA NA   6  2
2013-03-12 NA NA   5  3
2013-03-14  3  1  NA NA …
Run Code Online (Sandbox Code Playgroud)

merge r concatenation xts

7
推荐指数
1
解决办法
1万
查看次数

按事件拆分xts对象

我有一个带有少量事件指标的xts对象.我需要通过某些事件将其拆分,以便从给定事件到下一个事件的所有条目都将保存在同一个事件中xts,最终创建一个xts对象列表,每个对象都包含事件作为最后一个条目而没有其他事件相同的类型.

一个例子:

ts = as.Date(Sys.Date()-99:0)
e1 = numeric(100);e1[10*1:10]=1
e2 = numeric(100);e2[15*1:6]=1
y  = 1:100  # just a sample content
xs = as.xts(cbind(e1,e2,y),order.by=ts)

ee = e1*e2==1  # the event in which both e1 and e2 are 1, should happen at 30,60,90

# here should be splitting function that gets xs and ee as parameters
# and should return a list of 4 xts: the first with the entries 1 through 30, 
# the second with entries 31 …
Run Code Online (Sandbox Code Playgroud)

split r xts

6
推荐指数
1
解决办法
272
查看次数

列出给定大小的向量的所有子集

该函数choose(n,k)告诉我们k一组n不同元素存在多少个大小子集.假设我需要实际列出这些子集,我是如何创建它的?换句话说,我正在寻找一个接受矢量x(长度n)和数字的函数,k并返回一个矢量列表,每个矢量的大小k,子集的子集x.当然,列表的长度应该是choose(length(x),k).例如

enum.choose = function(x,k) {
    # implementation should go here
{

enum.choose(1:3,2)

# should be:
# [[1]]
#    1  2
# [[2]]
#    1  3
# [[3]]
#    2  3
Run Code Online (Sandbox Code Playgroud)

r list enumerate subset

6
推荐指数
1
解决办法
5941
查看次数

在R中绑定外部变量

假设我有以下功能:

g = function(x) x+h
Run Code Online (Sandbox Code Playgroud)

现在,如果我在我的环境中有一个名为的对象h,我就不会有任何问题:

h = 4
g(2)

## should be 6
Run Code Online (Sandbox Code Playgroud)

现在,我有另一个功能:

f = function() {
    h = 3
    g(2)
}
Run Code Online (Sandbox Code Playgroud)

我希望:

rm(h)
f()

## should be 5, isn't it?
Run Code Online (Sandbox Code Playgroud)

相反,我得到一个错误

## Error in g(2) : object 'h' not found
Run Code Online (Sandbox Code Playgroud)

我期望g在环境中进行评估f,以便hin f将被绑定到h in g,就像我g在内部执行时一样.GlobalEnv.这显然不会发生(显然).任何解释为什么?如何克服这个问题,以便g使用封闭环境评估函数(例如)中的函数?

environment binding r function

6
推荐指数
1
解决办法
131
查看次数

查询雅典娜时将结构转换为json

我有一个雅典娜表,我没有创建或管理它,但可以查询。字段之一是结构类型。为了这个例子,让我们假设它看起来像这样:

my_field struct<a:string,
                b:string,
                c:struct<d:string,e:string>
                >
Run Code Online (Sandbox Code Playgroud)

现在,我知道如何查询此结构中的特定字段。但是在我的查询之一中,我需要提取完整的结构。所以我只用:

select my_field from my_table
Run Code Online (Sandbox Code Playgroud)

结果看起来像一个字符串:

{a=aaa, b=bbb, c={d=ddd, e=eee}}
Run Code Online (Sandbox Code Playgroud)

我想将结果作为json字符串:

{"a":"aaa", "b":"bbb","c":{"d":"ddd", "e":"eee"}}
Run Code Online (Sandbox Code Playgroud)

然后,该字符串将由另一个应用程序处理,这就是为什么我需要json格式的原因。

我该如何实现?

编辑:更好的是,有没有办法以扁平化的方式查询该结构?所以结果看起来像:

a   |   b   |   c.d  |  c.e   |
-------------------------------
aaa |   bbb |   ddd  |  eee   |
Run Code Online (Sandbox Code Playgroud)

sql json struct amazon-athena

6
推荐指数
1
解决办法
2441
查看次数