1

我正在尝试解析使用.step函数生成的 LTSpice 数据。这会将生成的数据集拆分为我不知道如何有效处理的部分R。

library(tidyverse)

data <- read_tsv("FileName")
## 
## ── Column specification ────────────────────────────────────────────────────────
## cols(
##   time = col_character(),
##   `V(out)` = col_double(),
##   `V(in)` = col_double()
## )
## Warning: 24 parsing failures.
## row col  expected    actual                           file
##   1  -- 3 columns 1 columns 'FileName'
##   6  -- 3 columns 1 columns 'FileName'
##  11  -- 3 columns 1 columns 'FileName'
##  16  -- 3 columns 1 columns 'FileName'
##  21  -- 3 columns 1 columns 'FileName'
## ... ... ......... ......... ..............................
## See problems(...) for more details.

head(data)
## # A tibble: 6 x 3
##   time                                                          `V(out)` `V(in)`
##   <chr>                                                            <dbl>   <dbl>
## 1 Step Information: Resistor1=10 Resistor2=1M Offset=-3  (Run…  NA       NA     
## 2 0.000000000000000e+00                                         -1.50e-5  0     
## 3 4.913626562387649e-09                                          3.08e-2  0.0309
## 4 4.996634374887648e-09                                          3.13e-2  0.0314
## 5 5.000000000000000e-09                                          3.14e-2  0.0314
## 6 Step Information: Resistor1=10.01K Resistor2=1M Offset=-3  …  NA       NA

如您所见，数据被分成几部分，这些部分由有关步骤的信息分隔。我想更改数据，以便删除包含“步骤信息”的行，并将相应的参数写入匹配的列：

## # A tibble: 6 x 6
##   time                    `V(out)` `V(in)` Resistor1 Resistor2 Offset
##   <chr>                      <dbl>   <dbl> <lgl>     <lgl>     <lgl> 
## 1 0.000000000000000e+00 -0.0000150  0      NA        NA        NA    
## 2 4.913626562387649e-09  0.0308     0.0309 NA        NA        NA    
## 3 4.996634374887648e-09  0.0313     0.0314 NA        NA        NA    
## 4 5.000000000000000e-09  0.0314     0.0314 NA        NA        NA    
## 5 0.000000000000000e+00 -0.00749    0      NA        NA        NA    
## 6 4.908743749887649e-09  0.00789    0.0308 NA        NA        NA

编辑：

部分原始数据：

time    V(out)  V(in)
Step Information: Resistor1=10 Resistor2=1M Offset=-3  (Run: 1/24)
0.000000000000000e+00   -1.498494e-05   0.000000e+00
4.913626562387649e-09   3.082234e-02    3.086832e-02
4.996634374887648e-09   3.134312e-02    3.138962e-02
5.000000000000000e-09   3.136424e-02    3.141076e-02
Step Information: Resistor1=10.01K Resistor2=1M Offset=-3  (Run: 2/24)
0.000000000000000e+00   -7.485015e-03   0.000000e+00
4.908743749887649e-09   7.887614e-03    3.083766e-02
4.996634374887648e-09   8.162769e-03    3.138962e-02
5.000000000000000e-09   8.173305e-03    3.141076e-02
Step Information: Resistor1=20K Resistor2=1M Offset=-3  (Run: 3/24)
...
4

1 回答 1

1

它总是有点难看，但我认为你最好使用原始数据。将行拉入，标记具有 Step 信息的行，将数据分离出来，然后合并：

## read everything in as raw lines
## for this example use this code
##    tmp <- readLines(textConnection(txt))
## for your real data use this code
tmp <- readLines("FileName")

## flag the Step Information rows
stl <- grepl("^Step", tmp)

## get all name=value pairs out of the Step Information rows and join together
nms <- regmatches(tmp[stl], gregexpr("\\S+(?=\\=)",    tmp[stl], perl=TRUE))
val <- regmatches(tmp[stl], gregexpr("(?<=\\=)(\\S+)", tmp[stl], perl=TRUE))
nvp <- do.call(rbind, Map(function(v,n) as.data.frame(as.list(v), col.names=n), val, nms))

## finally put it all together
cbind(
  read.table(text=tmp[!stl], header=TRUE),
  nvp[ cumsum(stl)[!stl][-1], ]
)

#            time        V.out.      V.in. Resistor1 Resistor2 Offset
#1   0.000000e+00 -1.498494e-05 0.00000000        10        1M     -3
#1.1 4.913627e-09  3.082234e-02 0.03086832        10        1M     -3
#1.2 4.996634e-09  3.134312e-02 0.03138962        10        1M     -3
#1.3 5.000000e-09  3.136424e-02 0.03141076        10        1M     -3
#2   0.000000e+00 -7.485015e-03 0.00000000    10.01K        1M     -3
#2.1 4.908744e-09  7.887614e-03 0.03083766    10.01K        1M     -3
#2.2 4.996634e-09  8.162769e-03 0.03138962    10.01K        1M     -3
#2.3 5.000000e-09  8.173305e-03 0.03141076    10.01K        1M     -3

我使用的示例数据是：

txt <- "time\tV(out)\tV(in)
Step Information: Resistor1=10 Resistor2=1M Offset=-3  (Run: 1/24)
0.000000000000000e+00\t-1.498494e-05\t0.000000e+00
4.913626562387649e-09\t3.082234e-02\t3.086832e-02
4.996634374887648e-09\t3.134312e-02\t3.138962e-02
5.000000000000000e-09\t3.136424e-02\t3.141076e-02
Step Information: Resistor1=10.01K Resistor2=1M Offset=-3  (Run: 2/24)
0.000000000000000e+00\t-7.485015e-03\t0.000000e+00
4.908743749887649e-09\t7.887614e-03\t3.083766e-02
4.996634374887648e-09\t8.162769e-03\t3.138962e-02
5.000000000000000e-09\t8.173305e-03\t3.141076e-02"
于 2021-05-19T00:41:45.390 回答