bucles – como hacer/ modificar este loop en R?

Tengo una base de datos con 121 filas y como 10 columnas. Una de estas columnas corresponde las estaciones, otra a profundidad y el resto a variables químicas (temperatura, salinidad, etc.). Quiero calcular el valor integrado de estas propiedades químicas por estación, usando la función de oce:: integrationTrapezoid. Es la primera vez que hago un loop, así que no sé cómo. ¿me podrían ayudar/orientar en cómo hacerlo? Gracias.

dA<-matrix(data=NA, nrow=121, ncol=3) for (Station in unique(datos$Station)) {dA(Station, cd) <- integrateTrapezoid(cd, Profundidad..m., "cA") }

Station Depth temp
1 10 28
1 50 25
1 100 15
1 150 10
2 9 27
2 45 24
2 98 14
2 152 11
3 11 28.7
3 48 23
3 102 14
3 148 9