search for: muggeo

...R] B-spline/smooth.basis derivative matrices > To: "Vassily Shvets" <shv736 at yahoo.com> > Received: Monday, March 12, 2012, 5:15 PM > > > > -------- ???????????? ????????? > -------- > ?? ????: aleksandr shfets <a_shfets at mail.ru> > ????: vito.muggeo <vito.muggeo at unipa.it> > ????: 13 ????? 2012, 03:09 > ????: Re[2]: [R] B-spline/smooth.basis derivative > matrices > > > Hello again,(my message was truncated the first time!) > > Thank you, vito for these functions: I've been able to look > at the smoothi...

...ot;a.xx"),value=TRUE) [1] "a.x" "b.x" "a.xx" Instead, I would like to obtain only "a.x" "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

Hello, R users, I applied segmented regression method contributed by Muggeo and got different slope estimates depending on the initial break points. The results are listed below and I'd like to know what is a reasonable approach handling this kinds of problem. I think applying various initial break points is certainly not a efficient approach. Is there any other method...

...formula(lm(exp(y)~., data=d)) #--> exp(y) ~ age variable 'exp' (maybe indicating "experience") is not included in the model. The same happens with 'log' (and other function names, I suppose..) best, vito -- ============================================== Vito M.R. Muggeo Dip.to Sc Econom, Az e Statistiche Universit? di Palermo viale delle Scienze, edificio 13 90128 Palermo - ITALY tel: 091 23895240 fax: 091 485726 http://dssm.unipa.it/vmuggeo Associate Editor, Statistical Modelling Chair, Statistical Modelling Society

...,574,582,574,471,595, 557,557,584,599,517,649,584,463,591,488,563,553,549) library(MASS) library(Design) summary(polr(factor(grade)~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

Dear all, I would like to write a function like: myfun<-function(x,fn) {xx<-exp(x); x*fn(xx)} where fn is a symbolic description of any function with its argument to be specified. Therefore myfun(5,"2+0.3*y^2") should return 5*(2+0.3*exp(5)^2), myfun(5,"log(y)") should return 5*log(exp(5)) and so on. I tried with "expression" and others, but without success.

...e names.. Here a simple example all.vars(as.formula(y~poly(x,k)+z)) returns [1] "y" "x" "k" "z" and I would like to obtain "y" "x" "z" Where is the trick? 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 23895240 fax: 091 485726/485612 http://dssm.unipa.it/vmuggeo

Functions are first class objects, so some kind of collision is bound to happen if you do this... so don't. -- Sent from my phone. Please excuse my brevity. On January 30, 2018 3:11:56 AM PST, "Vito M. R. Muggeo" <vito.muggeo at unipa.it> wrote: >dear all, >Is the following intentional? Am I missing anything in documentation? > >d<-data.frame(y=rnorm(10,5,.5),exp=rnorm(10), age=rnorm(10)) >formula(lm(exp(y)~exp+age, data=d)) >#--> exp(y) ~ exp + age > >formula(lm(ex...

...) > fo . ~ . + a.log(z) > update(o,formula=fo) Error in eval(expr, envir, enclos) : could not find function "a.log" > 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

dear all, Given a matrix A, say, I would like to apply a bivariate function to each combination of its colums. That is if myfun<-function(x,y)cor(x,y) #computes simple correlation of two vectors x and y then the results should be something similar to cor(A). I tried with mapply, outer,...but without success Can anybody help me? many thanks in advance, vito

...ff Newmiller <jdnewmil at dcn.davis.ca.us> wrote: > > Functions are first class objects, so some kind of collision is bound to happen if you do this... so don't. > -- > Sent from my phone. Please excuse my brevity. > > On January 30, 2018 3:11:56 AM PST, "Vito M. R. Muggeo" <vito.muggeo at unipa.it> wrote: > >dear all, > >Is the following intentional? Am I missing anything in documentation? > > > >d<-data.frame(y=rnorm(10,5,.5),exp=rnorm(10), age=rnorm(10)) > >formula(lm(exp(y)~exp+age, data=d)) > >#--> exp(y) ~ ex...

...uot; not found namely myfun() does not work correctly. In particular the non-working line inside myfun() is update(eval(x$call$obj), 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

...-------------------------------------------------------------------------- INDICES: 1 V1 V2 -0.23682173 -0.01225147 > 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

...I would like to remove all entries within parentheses.. Namely, I aim to obtain respectively "z:f()" or "f()*z" I played with sub() and gsub() but without success.. Thank you very much for your time, best, vito -- ============================================== Vito M.R. Muggeo Dip.to Sc Econom, Az e Statistiche Universit? di Palermo viale delle Scienze, edificio 13 90128 Palermo - ITALY tel: 091 23895240 fax: 091 485726 http://dssm.unipa.it/vmuggeo Associate Editor, Statistical Modelling Chair, Statistical Modelling Society

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,

...nterested in), as call() *does* evaluate the arguments: > my.cc$random $Subject ~age > fm1$call$random list(Subject = ~age) How is it possible to get the right call (similar to the one from fm1$call) by means of call()? 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 23895240 fax: 091 485726/485612 http://dssm.unipa.it/vmuggeo

Dear all, The help file for the generic function vcov states "Classes with methods for this function include: 'lm', 'glm', 'nls', 'lme', 'gls', 'coxph' and 'survreg' (the last two in package 'survival')." Since, I am not able to use vcov.coxph(), I am wondering whether I am missing something (as I suspect..) regards, vito

Hi all, I am using the package ?Segmented? to estimate logistic regression models with unknown breakpoints (see Muggeo 2003 Statistics in Medicine 22:3055-3071). In the documentation it suggests that it might be possible to include several variables with breakpoints in the same model: ?Z = a vector or a matrix meaning the (continuous) explanatory variable(s) having segmented relationships with the response?. Howeve...

...oes not Error in eval(expr, envir, enclos) : object "y" not found > Could anyone to explain such difference? I believed 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

I appreciate to Roger Koenker, Achim Zeileis and Vito Muggeo for their informative answers. Listed below is unedited replies I got followed by the question I posted. Kyong 1. Roger Koenker: You might try rqss() in the quantreg package. It gives piecewise linear fits for a nonparametric form of median regression using total variation of the derivative...