Hola,
Una forma de hacerlo es esta:
#-------------------------> my_pol <- function( n = 100, grad_ini = 1, grad_end = 5) {
+ val_rnd <- rnorm(n)
+ df <- data.frame(c(0))
+ for( i in grad_ini:grad_end) {
+ df_tmp <- as.data.frame(val_rnd^i)
+ df <- cbind(df, df_tmp)
+ }
+ df <- df[ , 2:ncol(df)]
+ names(df) <- paste("V_", grad_ini:grad_end, sep = "")
+ return(df)
+ }>
> df_out <- my_pol(10, 3, 8)
> head(df_out)
V_3 V_4 V_5 V_6 V_7
V_8
1 -2.152566721 2.7793431004 -3.588622e+00 4.633543e+00 -5.982721e+00
7.724748e+00
2 -3.848810314 6.0316390414 -9.452445e+00 1.481334e+01 -2.321463e+01
3.638067e+01
3 0.681408379 0.5996207120 5.276498e-01 4.643174e-01 4.085866e-01
3.595450e-01
4 -0.003075413 0.0004472372 -6.503878e-05 9.458164e-06 -1.375439e-06
2.000211e-07
5 -0.620395129 0.5291244605 -4.512813e-01 3.848901e-01 -3.282662e-01
2.799727e-01
6 0.023210872 0.0066209750 1.888654e-03 5.387446e-04 1.536786e-04
4.383731e-05
#-------------------------
"caret" tiene un par de funciones "SLC14_1()" y
"SLC14_2" con las que
puedes generar dataframes con diferentes variables correladas entre ellas,
te sugiero que las veas por si además de esta función que te incluyo, te
ayuda en otros posibles casos.
Gracias,
Carlos Ortega
www.qualityexcellence.es
El sáb., 23 nov. 2019 a las 12:11, Manuel Mendoza (<
mmendoza en fulbrightmail.org>) escribió:
> Buenos días erreros, ¿sabe alguien cómo generar una df cuyas variables
> sigan una ley de potencias?
> Gracias
>
> [[alternative HTML version deleted]]
>
> _______________________________________________
> R-help-es mailing list
> R-help-es en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
--
Saludos,
Carlos Ortega
www.qualityexcellence.es
[[alternative HTML version deleted]]