[R-es] pawer law

[Carlos Ortega]

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 using 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 using r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>


-- 
Saludos,
Carlos Ortega
www.qualityexcellence.es

	[[alternative HTML version deleted]]