search for: 6626240

...; "b.x" How is it possible to get this result with grep()? many thanks for your attention, best, vito -- ==================================== Vito M.R. Muggeo Dip.to Sc Statist e Matem `Vianelli' Universit? di Palermo viale delle Scienze, edificio 13 90128 Palermo - ITALY tel: 091 6626240 fax: 091 485726/485612 http://dssm.unipa.it/vmuggeo

...e)~score)) lrm(factor(grade)~score) It seems to me that results from lrm() are right.. Am I wrong? Thanks, vito -- ==================================== Vito M.R. Muggeo Dip.to Sc Statist e Matem `Vianelli' Universit? di Palermo viale delle Scienze, edificio 13 90128 Palermo - ITALY tel: 091 6626240 fax: 091 485726/485612 http://dssm.unipa.it/vmuggeo

...data=d) #d is created in myfun() My question is: why the problem occurs just via source()ing??? many thanks, vito ==================================== Vito M.R. Muggeo Dip.to Sc Statist e Matem `Vianelli' Universit? di Palermo viale delle Scienze, edificio 13 90128 Palermo - ITALY tel: 091 6626240 fax: 091 485726/485612

...47 > by(H[,-1],H[,1],median) Error in median.default(data[x, , drop = FALSE], ...) : need numeric data > -- ==================================== Vito M.R. Muggeo Dip.to Sc Statist e Matem `Vianelli' Universit? di Palermo viale delle Scienze, edificio 13 90128 Palermo - ITALY tel: 091 6626240 fax: 091 485726/485612 http://dssm.unipa.it/vmuggeo

...that "with(data,glm(formula))" and "glm(formula,data)" were equivalent. Many thanks, vito -- ==================================== Vito M.R. Muggeo Dip.to Sc Statist e Matem `Vianelli' Universit? di Palermo viale delle Scienze, edificio 13 90128 Palermo - ITALY tel: 091 6626240 fax: 091 485726/485612

With version 8 of acrobat reader, it is now possible to have 3D in PDf documents. Does it exist already an R package who manage to produce 3D plots which can be saved as interactive 3D graphs in a PDF file? Best Regards Bruno Cavestro ------------------------------------------------------ Leggi GRATIS le tue mail con il telefonino i-mode? di Wind http://i-mode.wind.it/

...> How can fit the model? namely how can I use "a.log(z)" in the example above? Many thanks, vito -- ==================================== Vito M.R. Muggeo Dip.to Sc Statist e Matem `Vianelli' Universit? di Palermo viale delle Scienze, edificio 13 90128 Palermo - ITALY tel: 091 6626240 fax: 091 485726/485612

Hi all, Could anybody recommend some easy-to-understand and example based notes/tutorials on how to use cubic splines to do filtering on multi-dimension data? I am confused by the 1-dimensional case, and more confused by multi-dimensional case. I found all the books suddenly become very abstract when it comes to this subject. They don't provide examples in R or Splus at all. Specifically,

I have estimated a multiple nonparametric regression using the loess command in R. I have also estimated an additive version of the model using the gam function. Is there a way of using the output of these two models to test the restrictions imposed by the additive model?

Hi. I'm a student at SFU in Canada. The basic thing I want to do is calculate means of different strata. I have 2 vectors. One has the values I want to take the means from, the other is the four strata I am interested in. So I essentially want to break up the information vector into the four strata and calculate four means, one for each stratum. How can I do this in a reasonable way? Thanks

Is there a way to fit mixed proportional odds models in R? -- Using Opera's revolutionary e-mail client: http://www.opera.com/m2/

Hi, I am trying to creat a matrix with diagnal blocks, say, I have a matrix X of any dimension (nxm) ,and would like to have: X X X

Is it possible to use the '\ell' (i.e. the log likelihood) in plots? I've been browsing the plotmath documentation unsucesfully. Cheers, Mario dos Reis mdosrei at nimr.mrc.ac.uk +44 (0)20 8816 2300 Division of Mathematical Biology National Institute for Medical Research The Ridgeway Mill Hill London, NW7 1AA, UK

...11 day 26 svn rev 43537 language R version.string R version 2.6.1 (2007-11-26) > -- ==================================== Vito M.R. Muggeo Dip.to Sc Statist e Matem `Vianelli' Universit? di Palermo viale delle Scienze, edificio 13 90128 Palermo - ITALY tel: 091 6626240 fax: 091 485726/485612

I'm using the monoreg function (with weights) from the fdrtool package. How can I calculate the R square for this type of regression? Thanks for your help, Thierry -- View this message in context: http://www.nabble.com/R-square-for-Monotone-regression-tp15580803p15580803.html Sent from the R help mailing list archive at Nabble.com.

One of the assumptions, on which the (General) Linear Modelling is based is that the response variable is measured with some uncertainties (or weighted), but the explanatory variables are fixed. Is it possible to extend the model by assigning the weights to the explanatory variables as well? Is there a package for doing such a model fit? Thanks

Hello, I'm trying to find an algorithm to estimate a switching regression model based on the 1990 Economics Letters paper by Ohtani/Kakimoto/Abe or the earlier version from 1985 (Ohtani/Katayama, Economic Studies Quarterly; assuming as a transition path a polynomial of order 1). I found an idea for using nls here: http://www.biostat.wustl.edu/archives/html/s-news/2000-04/msg00223.html.

Dear All, I have a matrix say, X ( 100 X 40,000) and a vector say, y (100 X 1) . I want to perform linear regression. I have scaled X matrix by using scale () to get mean zero and s.d 1 . But still I get very high values of regression coefficients. If I scale X matrix, then the regression coefficients will bahave as a correlation coefficient and they should not be more than 1.

Hello R users I ran more than 100 logistic regression analyses. Some of the analyses gave me this kind warning below. ########################################################### Warning messages: 1: algorithm did not converge in: glm.fit(x = X, y = Y, weights = weights, start = start, etastart = etastart, ... 2: fitted probabilities numerically 0 or 1 occurred in: glm.fit(x = X, y = Y,

Hi all, I'm having sometrouble with managing the seach path, in a function , I need to attach some data set at the begining and detach them at the end, say, myfunction<- function() { attach(mylist); ............. detach(mylist) } , the problem is, since I am still debugging this code, sometimes it got error and ended before reaching the end, thus the data is left in the