nav[aria-label="Primary Navigation"] { padding: 0; & ul { list-style: none; width: 100%; display: flex; flex-direction: row; justify-content: start; align-items: start; gap: 30px; padding: 0; & li { margin: 0; } & ul li { list-style: none; } } }

274941_it-IT

liliconcom

Individuazione del vertice del grafico bivariato

Ciao a tutti,

Sto cercando di ottenere un punto di vertice dal grafico bivariato

finora, ottengo il parametro extract con script e lo calcolo per foglio Excel.

Sono abbastanza sicuro che la funzione di vertice sia inclusa in JMP ma non sono riuscita a trovarla.

C'è un modo per ottenere direttamente il vertice?

Grazie

Ri: Trovare il vertice del grafico bivariato

JMP non ha un risolutore in quanto tale. Ma stai adattando una funzione di interpolazione polinomiale quadratica, quindi devi solo impostare la prima derivata rispetto a Offset_Z uguale a zero e risolvere per Offset_Z. Suppongo che per "vertice" intendi il minimo o il massimo della funzione.

Names Default To Here( 1 );

// example of data
dt = Current Data Table();

// fit quadratic polynomial interpolating function
biv = dt << Bivariate(
 Y( :contrast_100cm ),
 X( :Offset_Z ),
 Fit Polynomial( 2, {Line Color( {212, 73, 88} )} ),
 Where( :site == "Para#4" )
);

// get reference to report layer
rbiv = biv << Report;

// get parameter estimates
parms = rbiv[NumberColBox( 1 )] << Get As Matrix;

// compute extremum
extremum = -parms[2] / (2*parms[3]);

Show( parms, extremum );

Ri: Trovare il vertice del grafico bivariato

In alternativa, un altro approccio consiste nell'utilizzare il modello di adattamento, utilizzare la macro della superficie di risposta in modo che JMP imposti l'attributo della superficie di risposta. Quindi otterrai un output della superficie di risposta che risolve i valori critici.

Ho fatto creare a JMP questo script con alcuni dati che ho creato.

Fit Model(
 Y( :contrast_100cm ),
 Effects( :Offset_Z & RS, :Offset_Z * :Offset_Z ),
 Personality( "Standard Least Squares" ),
 Emphasis( "Effect Leverage" ),
 Run(
  :contrast_100cm << {Summary of Fit( 0 ), Analysis of Variance( 0 ),
  Parameter Estimates( 0 ), Lack of Fit( 0 ), Scaled Estimates( 0 ),
  Plot Actual by Predicted( 0 ), Plot Residual by Predicted( 0 ),
  Plot Studentized Residuals( 0 ), Plot Effect Leverage( 0 ),
  Plot Residual by Normal Quantiles( 0 ), Box Cox Y Transformation( 0 )}
 ),
 SendToReport(
  Dispatch(
   {"Response contrast_100cm", "Whole Model"},
   "Summary of Fit",
   OutlineBox,
   {Close( 1 )}
  ),
  Dispatch(
   {"Response contrast_100cm", "Whole Model"},
   "Analysis of Variance",
   OutlineBox,
   {Close( 1 )}
  ),
  Dispatch(
   {"Response contrast_100cm", "Whole Model"},
   "Parameter Estimates",
   OutlineBox,
   {Close( 1 )}
  ),
  Dispatch(
   {"Response contrast_100cm"},
   "Response Surface",
   OutlineBox,
   {Close( 0 )}
  )
 )
)

Ecco la parte dell'output che trova il valore critico.

Ri: Trovare il vertice del grafico bivariato

Grazie per la risposta.

Ho testato la sceneggiatura e ho scoperto che era esattamente quello che stavo cercando.

(Anche se non ho capito l'intero codice)

ma ho scoperto che il valore Offset_Z che viene calcolato dalla superficie di risposta è leggermente diverso dal mio calcolo manuale (per differenziale)

la differenza era inferiore a 0,0003, quindi non era un grosso problema

Grazie ancora per il tuo aiuto

Ri: Trovare il vertice del grafico bivariato

Grazie per la cortese risposta.

Poiché hai spiegato ogni codice per osservazione, è stato facile da capire.

Stavo cercando di modificare il tuo codice per scoprire il x intercetta in funzione lineare.

ma non è venuto fuori nulla.

In realtà non so davvero come funzioni "show function".

La prego di aiutarmi a risolvere questo problema?

Grazie

Ri: Trovare il vertice del grafico bivariato

Il Mostrare() La funzione scrive il testo in base ai suoi argomenti nella finestra Registro, che non è aperta per impostazione predefinita. Selezionare Visualizza> Registro se si utilizza Windows o Finestra> Registro se stai usando Macintosh.

Ho modificato lo script per stimare l'intercetta x, supponendo che esista sul piano reale. È solo la soluzione quadratica per le radici.

Names Default To Here( 1 );

// example of data
x = -5::5;
dt = New Table( "Non-Linear Data",
 New Column( "X", Numeric, Continuous, Values( x ) ),
 New Column( "Y", Numeric, Continuous, Values( x^2 - x - 12 ) )
);

// fit quadratic polynomial interpolating function
biv = dt << Bivariate(
 Y( :Y ),
 X( :X ),
 Fit Polynomial( 2 )
);

// get reference to report layer
rbiv = biv << Report;

// get parameter estimates
parm = rbiv["Parameter Estimates"][NumberColBox( 1 )] << Get As Matrix;

// compute x intercept
ans = parm[2]^2 - 4*parm[3]*parm[1];
x intercepts = If( ans > 0,
 (-parm[2] + { -1, 1 } * Sqrt( ans )) / (2*parm[3]);
);

Show( parm, x intercepts );

Ri: Trovare il vertice del grafico bivariato

Ciao Makrbailey,

Posso chiedere più cose?

Ho leggermente modificato lo script per risolvere la funzione lineare e ha funzionato correttamente.

tuttavia, mi sono trovato che controllare il risultato attraverso il registro non è conveniente tanto quanto controllare attraverso la casella ad albero.

per quanto ne so, la funzione di invio è quella giusta per creare una casella ad albero da segnalare.

così ho cambiato lo script come segue. tuttavia, non funziona.

e voglio ottenere la risposta per ogni sito (ci sono 4 siti). quindi ho provato a farlo, ma non ha funzionato.

Potresti darmi la tua mano per risolverlo, per favore

Grazie mille

-> Il mio scopo

1. aggiungi ans di ogni sito.

2. aggiungi ans per segnalare.

-> La mia sceneggiatura

Nomi predefiniti qui (1);

dt = Current Data Table ();

// adatta la funzione di interpolazione polinomiale quadratica
biv = dt << Bivariato (
Y (: CALU_CALU_intr_cy),
X (: OC offset Y um),
Adatta linea ({Colore linea ({212, 73, 88})}),
by (: sito),
);

// ottiene riferimento al livello di report
rbiv = biv << Rapporto;

// ottiene le stime dei parametri
parm = rbiv ["Stime dei parametri"] [Casella numerica (1)] << Ottieni come matrice;

// calcola x intercetta
ans = -parm [1] / parm [2];

rbiv << spedizione (
{ "ans"},
"ans",
framebox,
{dimensione marker (3)}
);

Grazie

Ri: Trovare il vertice del grafico bivariato

Come ho pubblicato molte volte, il Invia al rapporto () e Spedizione() gli argomenti sono destinati a JMP per salvare la personalizzazione. Non sono destinati agli utenti. Certo, puoi usarli. Ma ci sono modi migliori.

Sembra che tu voglia aggiungere la risposta alla trama. Questo esempio illustra come farlo.

Names Default To Here( 1 );

dt = Open( "$SAMPLE_DATA/Big Class.jmp" );

// fit quadratic polynomial interpolating function
biv = dt << Bivariate(
 Y( :weight ),
 X( :height ),
 Fit Line
);

// get reference to report layer
rbiv = biv << Report;

// get parameter estimates
parm = rbiv["Parameter Estimates"][Number Col Box( 1 )] << Get As Matrix;

// compute x intercept
ans = -parm[1] / parm[2];

rbiv[FrameBox(1)] << Add Graphics Script(
 Text( { 52.5, 160 }, "X intercept is " || Char( ans ) )
);