Aiuto con S-PLUS

[mingus]

Qualcuno bravo ad usare il software S-PLUS o in alternativa R è disposto a darmi una mano con un piccolo problema?

Speriamo che ci sia qualcuno...

Grazie comunque
Andrea

[Slater]

...non so quanta gente troverai che smanetta con software di data mining...

[mingus]

Posso immaginarlo, ma magari photo4u mi stupisce

Sto lavorando sulla tesi e S-Plus mi serve perchè è una forza nel massimizzare funzioni di verosimiglianza piuttosto complesse. E' un lavoro interessante ma che dal punto di vista computazionale è da mani nei capelli


Andrea

[Slater]

guarda, non so se posso esserti d'aiuto ma le ho studiate in basi di dati le data mining... se vuoi provare a postare qualcosa vedo se riesco ad aiutarti...ma non garantisco niente

[mingus]

Il problema è nella programmazione in S-Plus.
Ho una routine che presi n valori li "sbatte" dentro una funzione e ne plotta l'andamento. Sullo stesso grafico ho anche una retta parallela all'asse delle ascisse che intercetta la prima curva 2 volte. Io devo trovare ardinata e ascissa dei 2 punti. Per l'ordinata non è un problema, il problema è nel determinare le ordinate. La funzione *non* è invertibile, quindi non posso buttargli dentro l'ordinata sperando che mi restituisca le ascisse.
E' proprio un problema di programmazione, che non è il mio campo. A differenza delle funzioni di verosimiglianza che ormai me le mangio a colazione

Il codice comunque è questo:

[Slater]

non conoscendo il linguaggio non è che possa aiutarti molto.
magari sarò banale, ma hai provato a vedere se non esiste una funzione che dati in input le 2 curva/retta ti resituisca i punti di intersezione? perchè ho provato a pensare qualcosa ma non puoi manipolare 2 variabili matematicamente (nel senso pensando alla matematica non puoi applicare il principio di sostituzione su 2 variabili che hanno assegnate le due curve/retta)
....sorry.... di + non so come aiutarti

[photolais]

v[i] e x[i] rappresentano le coordinate dei punti della curva ?
se si, allora forse puoi fare così (i dettagli dipendono dal linguaggio di programmazione che non so quale sia)
L'idea è che se le due curve si intersecano allora prima (dell'intersezione) una sta sopra, poi sotto.
Allora ad ogni ciclo calcolo la differenza di ordinata tra le due, quando cambia di segno ho il punto di intersezione (v[i],x[i])
Quindi nel dettaglio farei così:
utilizzo due variabili: deltaCorr e deltaPrec, inizializzate a +100000 (o comunque ad un valore sicuramente maggiore della differenza piu grande tra le due)
Nel ciclo che hai riportato, dopo l'ultima operazione (calcolo di v[i]) io calcolerei una variabile 'deltaCorr' così
delta1<-v[i]-yretta (cioè credo ma - 0.5 * qchisq(conf, 1))
if (deltaPrec=10000 ) or (deltaCorr*deltaPrec)>0 then
deltaPrec=deltaCorr
else
hai il punto cercato (v[i],x[i]) e te lo salvi o stampi come ti pare
ovviamente la sintassi può cambiare parecchio, il senso del test della if è il seguente:
deltaPrec=10000 --->siamo al primo giro, non ha senso considerare deltaPrec
(deltaCorr*deltaPrec)>0 i due delta hanno lo stesso segno--> le curve non si sono intersecate.
Spero ti possa servire
Fammi sapere se va, se non va o se servono maggiori chiarimenti

[mingus]

Grazie mille ad entrambi.

photolais, ad una soluzione simile alla tua ci avevo pensato pure io ora devo studiarmi un po' meglio il linguaggio. Non sono mica un informatico, io!

Se esce qualcosa vi farò sapere.


Andrea