小编Naz*_*zer的帖子

我相信我的数据帧是可以的,我的代码还可以.实际上,我已经删除了部分数据框和大部分图形代码,以使事情尽可能基本.但是,我仍然得到:

Error in .Call.graphics(C_palette2, .Call(C_palette2, NULL)) : 
  invalid graphics state

这有什么不对？这是数据:

 date   trt var val
1/8/2008    cc  sw5 0.2684138
1/8/2008    cc  sw15    0.2897586
1/8/2008    cc  sw5 0.2822414
2/8/2008    cc  sw5 0.2494583
2/8/2008    cc  sw5 0.2692917
2/8/2008    cc  sw15    0.2619167
2/8/2008    cc  sw5 0.204375
3/8/2008    cc  sw5 0.2430625
3/8/2008    cc  sw5 0.2654375
3/8/2008    cc  sw5 0.2509583
3/8/2008    cc  sw5 0.2055625
1/8/2008    ccw sw15    0.2212414
1/8/2008    ccw sw5 0.3613448
1/8/2008    ccw sw5 0.2607586
2/8/2008    ccw sw5 0.2087917
2/8/2008    ccw sw15    0.3390417
2/8/2008 …

当标题有两个必要的标题行时,将文件读入R的最佳方法是什么？

这种情况一直发生在我身上,因为人们经常使用一行作为列名,然后在其下面包含另一行作为测量单位.我不想跳过任何东西.我想要名字和单位进行.

以下是具有两个标头的典型文件可能如下所示:

trt   biomass    yield
crop    Mg/ha    bu/ac
C2      17.76   205.92
C2      17.96   207.86
CC      17.72   197.22
CC      18.42   205.20
CCW     18.15   200.51
CCW     17.45   190.59
P       3.09    0.00
P       3.34    0.00
S2      5.13    49.68
S2      5.36    49.72

我有一些数据必须格式化为(%d /%m /%Y).数据不按时间顺序排列,因为它按第一个数字排序,即第一个数字,而不是月份.

我希望我可以指定order或者reorder我希望排序以不同的方式发生.我只是不确定该怎么做.

以下是要订购的日期数据:

date
1/1/2009  
1/1/2010
1/1/2011
5/4/2009
5/4/2011
10/2/2009
10/3/2011
15/9/2010
15/3/2009
31/12/2011
31/7/2009

谢谢你的任何建议.

我想使用相同的ggplot代码生成8个不同的数字,条件是我的数据框中的数字.通常我会使用facet_grid,但在这种情况下,我想得到每个人的数字的pdf.例如,我想在这里每行一个pdf:

df <- read.table(text = "
xvalue     yvalue    location    planting    crop
  1          5          A          early      corn
  2          3          A          late       corn
  6          2          A          early      soy
  7          4          A          late       soy
  4          7          S          early      corn
  2          6          S          late       corn
  3          2          S          early      soy
  5          1          S          late       soy
", sep = "", header = TRUE)

基本ggplot:

library(ggplot2)

ggplot()+
  geom_point(aes(x=xvalue, y=yvalue), data=df)

但是不是facet_grid为了获得x种植x作物组合的位置,我想要一个单独的pdf.

我试图sum在函数中使用,但结果是NA,我认为可能是由于整数溢出.但我使用的数字类是数字.

功能最简单

sum((columnA-columnB)^2)

来自columnA的值是0.1376146来自columnB的值0.272

是小数位的不同长度？我知道如何更改显示的内容,但我不确定这会改变R使用的内容sum.

我正在尝试使用 xml2 读取 Excel xml 文件，但遇到了困难，因为我拥有的文件与 xml2 文档中使用的结构非常不同。我想阅读工作簿中的一张工作表并用作数据框。

该片段包含完整的结构，但只有一个单元格填充了一堆文本，而我想要读取的内容有 50,000 行数据。

<?xml version='1.0'?>
<?mso-application progid='Excel.Sheet'?>
<s:Workbook xmlns:x="urn:schemas-microsoft-com:office:excel" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:s="urn:schemas-microsoft-com:office:spreadsheet">
  <s:Worksheet s:Name="DBCitation">
    <s:Table>
      <s:Row>
        <s:Cell>
          <s:Data s:Type="String">The suggested citation for your download is below. See metadata folder and citationsyntax.xls for more explanation</s:Data>
        </s:Cell>
      </s:Row>
      <s:Row>
        <s:Cell>
          <s:Data s:Type="String" />
        </s:Cell>
      </s:Row>
      <s:Row>
        <s:Cell>
          <s:Data s:Type="String">Acosta-Martinez, Veronica ; Balkcom, Kipling; Caesar-TonThat, Thecan; Franzluebbers, Alan; Gollany, Hero; Jabro, Jalal; Jin, Virginia; Johnson, Jane; Liebig, Mark; Phillips, Rebecca; Sainju, Upendra; Sistani, Karamat; Skinner, R; …

我曾经lmList适合 480 个关系，我想要每个关系的 R2。

这是一个示例数据集和模型，它非常接近实际情况，但我有 480 eu（实验单位）：

eu  mass  day
11  .02    1
11  .03    2
11  .04    3
11  .06    4
12  .01    1
12  .03    2
12  .04    3
12  .05    4

fit<-lmList(mass ~ day | eu, data=df)

打印fit或summary不给我我想要的信息。我最终试图制作一个新的数据框，看起来像：

eu    intercept    slope    R2
11     .01          .95     .98
12     .01          .96     .98

我已经得到了系数coef，现在我需要 R 平方。

我正在拟合一些线条,我觉得我正在告诉R究竟如何适合它们,但我觉得有一些东西(一些因素或影响)我不知道这是在阻止一个良好的契合.

我的实验单位是"情节",如在田野情节中,我很抱歉令人困惑.

可以找到以下数据:https://www.dropbox.com/s/a0tplyvs8lxu1d0/rootmeansv2.csv.同

df$plot.f<-as.factor(df$plot)
dfG<-groupedData(mass ~ year|plot.f, data=df)
dfG30<-dfG[dfG$depth == 30,]

简单地说,随着时间的推移,我的质量很大,我将其与模型的每个实验单元相匹配:

fit <- lme(mass ~ year , random = ~ 1 | plot, data = df)

并且plot (augPred(fit))我得到了每个实验单位的这些拟合("情节"):

我需要做些什么来让实验单位之间的坡度变化更大？从统计角度来看,我对此并不感兴趣,但从预测的角度来看 - 所以模型中的任何东西都可以被操纵以使这些线条移动.

我试图通过基于标识符进行划分来改变数据.例如,我想mass在以下数据中进行变异.如果它是深度10,我想除以2.如果它是深度20,我想除以3.

day    year    depth    mass 
1      2008    10       13
2      2008    10       15
1      2008    20       14
2      2008    20       12
1      2009    10       14
2      2009    10       16
1      2009    20       12
2      2009    20       18   

分部导致:

day    year    depth    mass 
1      2008    10      6.5
1      2008    10      6.5
2      2008    10      7.5
2      2008    10      7.5
1      2008    20      4.6
1      2008    20      4.6
2      2008    20      4
2      2008    20      4
1      2009    10      7 …

我有一个数据框列表,所有数据框都有相同数量的列(和行).我想用它们绑定它们purrr::map_df.

我试着map_df(my_list)去

Error in as_mapper(.f, ...) : argument ".f" is missing, with no default

我不确定我的清单有什么问题.它看起来不错(每个数据帧都有一个唯一的名称):