関数を実装するのに手助けが必要ですR

Question

私はこの形式で出力を定期的に取得し、Excelで1つ1つコピー＆ペーストする必要があるため、スクリプトを活用するような方法でコードを作成しようとしています雑用。しかし、私はそれが機能を実装することになると固執しています。

だから、私のデータは次の形式である：

Condition Sample1 Sample2 .... Sample n 
T1  6.99 5.80 .... n_1  
T2  2.05 3.04 .... n_1  
T3  4.50 4.69 .... n_1  
T4  4.71 5.22 .... n_1  
T5  5.66 3.65 .... n_1  
T6  9.76 2.89 .... n_1  

私は、次の式に適用する必要があります。、xは、個々のエントリで、nは係数で、完全な方程式が何かに見えるようこのように：

基本的には、各要素を順番に考え、逐次係数（1：length条件からの奇数）を掛けて、答えをSとする必要があります。私のデータセットのサイズは変更されません。それは常にT1：T6になります。変更されるのはサンプル1 ... nです。理想的にはSの値が列の末尾に追加されるか、属しているサンプルを参照して別のデータセットに保存されることが理想的です。

私は転置を含むいくつかの解決策を試しましたが、私の周りに頭を包んでいるように見えません。

現在、データセットの一部でより単純な関数を実装しようとする私の試みは成功しませんでした。

適切なソリューションには*適用されていると思いますが、正しく使用する方法は完全にわかりません。

EDIT：再現可能な例を追加：

N Condition Sample A Sample B Sample C Sample D 

1 T1   91.323  78.758  70.298  66.765 
3 T2   -3.737  -1.5  -7.744  -9.247 
5 T3   5.205  4.533  2.284  2.178 
7 T4   -0.486  -0.068  -1.386  -0.927 
9 T5   0.337  -0.139  0.087  0.055 
    S  -0.046296296 -0.123654391 0.394039047 0.445258425 

Answer 1

あなたが何をしたいのか、これをしていますか？

per_row <- function(row){ 
    l <- length(row) 
    exp <- 2*(1:l)-1 # all the exponents 
    each <- row*(-1)^exp # compute all of these at once 
    return(sum(each)) # return sum 
} 

#some sample data 
datafr <- data.frame(a = sample(1:6), b = 1:6) 

#apply per column 
apply(datafr, 2, per_row) 

Answer 2

この内容はどうですか？あなたが持っているSampleN列の数には無関心でなければなりません。 6つの奇数乗算器を使用する6条件ケースに対して特別に設計されていることに注意してください。ただし、それは変更されないので、問題はないはずです。

suppressPackageStartupMessages(library(dplyr)) 
suppressPackageStartupMessages(library(tidyr)) 

samples <- tribble(
    ~Sample1, ~Sample2, ~Sample3, 
    6.99, 5.80, 2.5, 
    2.05, 3.04, 3.4,  
    4.50, 4.69, 8.7,  
    4.71, 5.22, 8.6,  
    5.66, 3.65, 3.4,  
    9.76, 2.89, 5.6 
) 

samples 
#> # A tibble: 6 x 3 
#> Sample1 Sample2 Sample3 
#>  <dbl> <dbl> <dbl> 
#> 1 6.99 5.80 2.50 
#> 2 2.05 3.04 3.40 
#> 3 4.50 4.69 8.70 
#> 4 4.71 5.22 8.60 
#> 5 5.66 3.65 3.40 
#> 6 9.76 2.89 5.60 

samples_modified <- samples %>% 
    # Add the multipliers as a column 
    mutate(multiplier = c(1,3,5,7,9,11)) %>% 

    # Gather all the samples. Make it 'tidy' 
    gather(key = "sample", value = "x", -multiplier) %>% 

    # Perform the multiplication on each element, we will sum later 
    mutate(x_modified = x * (-1)^((multiplier - 1)/2)) 

# Now we want to sum the x_modified column for each sample group 
samples_modified %>% 
    group_by(sample) %>% 
    summarise(S = sum(x_modified)) 
#> # A tibble: 3 x 2 
#> sample  S 
#> <chr> <dbl> 
#> 1 Sample1 0.630 
#> 2 Sample2 2.99 
#> 3 Sample3 -3.00 

Answer 3

簡単なようです。 Rはベクトル化されたベクトルの乗算であるためです。そして、sum。

zirconium <- function(x){ 
    n <- 2*seq_along(x) - 1 
    sum(x * (-1)^((n - 1)/2)) 
} 

sapply(dT[-1], zirconium) 
#Sample1 Sample2 
# 0.63 2.99 

データ。

dT <- 
structure(list(Condition = structure(1:6, .Label = c("T1", "T2", 
"T3", "T4", "T5", "T6"), class = "factor"), Sample1 = c(6.99, 
2.05, 4.5, 4.71, 5.66, 9.76), Sample2 = c(5.8, 3.04, 4.69, 5.22, 
3.65, 2.89)), .Names = c("Condition", "Sample1", "Sample2"), class = "data.frame", row.names = c(NA, 
-6L))