## With a couple of variables
carb(flag=8, var1=8.2, var2=0.00234, S=35, T=25, P=0, Pt=0, Sit=0)
## Using vectors as arguments
flag <- c(8, 2, 8)
var1 <- c(8.2, 7.477544e-06, 8.2)
var2 <- c(0.002343955, 0.001649802, 2400e-6)
S <- c(35, 35, 30)
T <- c(25, 25, 30)
P <- c(0, 0, 0)
Pt <- c(0, 0, 0)
Sit <- c(0, 0, 0)
carb(flag=flag, var1=var1, var2=var2, S=S, T=T, P=P, Pt=Pt, Sit=Sit)
## Test with all flags
flag <- c(1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 21, 22, 23, 24, 25)
var1 <- c(8.200000, 7.477544e-06, 7.477544e-06, 7.477544e-06, 7.477544e-06, 8.2,
8.2, 8.2, 8.2, 0.001685024, 0.001685024, 0.001685024, 0.0002888382, 0.0002888382,
0.002391252, 264.2008, 264.2008, 264.2008, 264.2008, 264.2008)
var2 <- c(7.477544e-06, 0.001685024, 0.0002888382, 0.002391252, 0.001981340,
0.001685024, 0.0002888382, 0.002391252, 0.001981340, 0.0002888382, 0.002391252,
0.001981340, 0.002391252, 0.001981340, 0.001981340, 8.2, 0.001685024, 0.0002888382,
0.002391252, 0.001981340)
carb(flag=flag, var1=var1, var2=var2)
## Test using a data frame
data(seacarb_test)
tab <- seacarb_test
## method 1 using the column numbers
carb(flag=tab[[1]], var1=tab[[2]], var2=tab[[3]], S=tab[[4]], T=tab[[5]], P=tab[[6]],
Sit=tab[[7]], Pt=tab[[8]])
## method 2 using the column names
carb(flag=tab$flag, var1=tab$var1, var2=tab$var2, S=tab$S, T=tab$T, P=tab$P, Sit=tab$Sit,
Pt=tab$Pt)
Run the code above in your browser using DataLab