530471_it-IT

liliconcom

Oggetto: JMP > JSL > Metodo migliore per calcolare l'area sotto la curva con segno per serie temporali semplici

Regola del trapezio?  

Names default to here(1);
dt = current data table();
x_col = Column(dt, "Time (Week)");
y_col = Column(dt, "Y");
by_col = Column(dt, "Curve");


auc = function({x_col, y_col, by_col}, 
 {DEFAULT LOCAL}, 
 Summarize(dt, by_values=By(by_col));
 areas = [=>];
 for(i=1, i<=nitems(by_values), i++, 
  by_value = by_values[i];
  rows = dt << Get Rows Where(by_col[]==by_value); // get he rows for the curve
  x_values = x_col[rows];
  r = rank(x_values); // make sure it's ordered right
  
  x_values = x_values[r];
  y_values = y_col[rows][r];
  n = nitems(x_values);
  //trapezoid rule
  x_1 = x_values[1::n-1];
  x_2 = x_values[2::n];
  y_1 = y_values[1::n-1];
  y_2 = y_values[2::n];
  
  x_intervals = abs(x_2 - x_1);
  y_avg = (y_2 + y_1)/2;
  area = sum(x_intervals :*</img> y_avg);
  areas[by_value] = area
 );
 return(areas);
);
areas = auc(x_col, y_col, by_col);
show(areas); // areas = ["A" => 10.5, "B" => -7.5, "C" => -16.5, "D" => 6, "E" => 24];

 

Inoltre, quello che stai facendo non è un poligono (nel caso lo usi altrove per cose simili).L'area di un poligono non può mai essere negativa o non esisterebbe.
L'esempio nell'indice di scripting mostra cosa sta facendo.  

 

Names Default To Here( 1 );
area = Polygon Area( {0, 0}, {0, 10}, {10, 10}, {10, 0} );

Chiude il ciclo (in questo caso un quadrato) tra {10, 0} e {0, 0}.Quindi l'area del poligono è 100.Non importa dove lo metti nel piano x/y.L'area di un quadrato 10x10 dovrebbe essere 100.

 

Names Default To Here( 1 );
area = Polygon Area( {-100, -100}, {-100, -90}, {-90, -90}, {-90, -100} );
show(area); //100

JMP > JSL > Metodo migliore per calcolare l'area sotto la curva con segno per serie temporali semplici

Ciao comunità JMP,

 

Ho bisogno di calcolare l'area firmata sotto la curva per semplici dati di serie temporali (vedi i dati fittizi di seguito).

Nello specifico, devo calcolare l'area (come mostrato nei grafici sopra) sopra lo zero come positiva e sotto lo zero come negativa, compreso lo spostamento locale da positivo a negativo e viceversa.

 

Esempi:

• AUC per la curva A = Area 1 positiva + Area 2 negativa + Area 3 positiva
• AUC per la curva C = Area 1 positiva + Area 2 negativa
• AUC per la curva E = Area positiva 1

La funzione Polygon Area () non sembra gestire correttamente le serie temporali più complesse, probabilmente a causa dei segmenti che si incrociano. Inoltre, questa funzione non restituisce il valore con segno per l'area. Conosceresti un componente aggiuntivo o una funzione che potrei utilizzare per raggiungere il mio obiettivo?

Grazie.

Migliore,

ST

Modellazione statistica avanzataQualità e ingegneria di processo