search for: vito

...een ordered Fit: aov(formula = PAPositionPercentOfVoweldur ~ UtteranceType * SyllLable * Cycle * Speaker * Label, data =PABTSub) $`Cycle:Speaker:Label` diff lwr upr p adj 3:Andrea:!H*L-1:Lavinia:!H*L 2.6069140 -37.499300 42.71313 1.0000000 1:Vito:!H*L-1:Lavinia:!H*L 8.7764075 -85.090794 102.64361 1.0000000 3:Andrea:H*L-1:Lavinia:!H*L 12.3411960 -18.883688 43.56608 0.9792814 1:Vito:H*L-1:Lavinia:!H*L 15.0416018 -32.746962 62.83017 0.9968844 1:Andrea:H*L-1:Lavinia:!H*L 15.2934976 -17.987036 48.57403 0.9379977 1:Lavinia:H*L-1:La...

...y 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. How can I solve my problem? Many thanks, vito

...llection, after two or three iterations, R crashes when used interactively. I thought that running and exiting R after each iteration would be the only effective way to "clean up" memory. If you have any other solutions, I would be happy to hear about them as well. Laurent Franckx, PhD VITO NV Boeretang 200, 2400 MOL, Belgium Tel. + 32 14 33 58 22 Skype: laurent.franckx laurent.franckx at vito.be Visit our website: www.vito.be/english and http://www.vito.be/transport ________________________________ [http://www.vito.be/e-maildisclaimer/bsds.png] <http://bsds.vito.be> VI...

Dear UseRs, maybe is a silly question: how can I get Empirical CDF values from an object created with ecdf()?? Using print I obtain: Empirical CDF Call: ecdf(t) x[1:57] = 4.1, 4.4, 4.5, ..., 491.3, 671.27 Thanks in advance. Regards, Vito Diventare costruttori di soluzioni Became solutions' constructors "The business of the statistician is to catalyze the scientific learning process." George E. P. Box "Statistical thinking will one day be as necessary for efficient citizenship as the ability to read and writ...

...rently running an R program on a computer with 16 Gb memory (Windows7, 64 bit). When I look at the task manager, I see that only 4 out of the 8 CPUs are being used. Is this due to some missing in the R code, or should I change something to the settings of the computer? Laurent Franckx, PhD Expert VITO NV Boeretang 200, 2400 MOL, Belgium Tel. + 32 14 33 58 22 Skype: laurent.franckx laurent.franckx at vito.be Visit our website: www.vito.be/english and http://www.vito.be/transport http://www.vito.be/e-maildisclaimer

...grep("x",c("a.x" ,"b.x","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

...2]: [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...

...ls\bin; C:\R\Rtools\gcc-4.6.3\bin64; C:\MiKTeX\miktex\miktex\bin; c:\windows; c:\windows\system32; C:\Windows\SysWOW64; C:\Program Files (x86)\Inno Setup 5. I have re-installed today the most recent version of Rtools. Does anyone have an idea what is going wrong here? Laurent Franckx, PhD Expert VITO NV Boeretang 200, 2400 MOL, Belgium Tel. + 32 14 33 58 22 Skype: laurent.franckx laurent.franckx at vito.be Visit our website: www.vito.be/english and http://www.vito.be/transport http://www.vito.be/e-maildisclaimer

...ore<-c(525,533,545,582,581,576,572,609,559,543,576,525,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

...10)) formula(lm(exp(y)~exp+age, data=d)) #--> exp(y) ~ exp + age 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, Statist...

...have googled the error message. Initially I had omitted “--arch x64” when compiling the file, but even after specifying the architecture, the error persists. However, I definitely need to be able to run my code on the 64 version. Kind regards Laurent ________________________________ [http://www.vito.be/e-maildisclaimer/bsds.png] <http://bsds.vito.be> VITO Disclaimer: http://www.vito.be/e-maildisclaimer [[alternative HTML version deleted]]

...iate 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

...sefell for me to import/export data from/to Excel: importa.da.excel<-function(){read.delim2("clipboard", dec=",") ## questa funzione consente di importare dati da Excel in R ## selezionare in Excel le celle che contengono i dati, ## compresi in nomi delle colonne ## Autore: Vito Ricci email:vito_ricci at yahoo.com ## Data di creazione: 09/11/04 } esporta.in.excel<-function(dati){write.table(dati,"clipboard", sep="\t", dec=",", col.names=NA) ## questa funzione consente di esportare dati in Excel da R ## passare come argomento il dataframe,...

...(".~.+",paste(myname, collapse = "+"))) > 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

...e all.vars(), I am interest in extracting only the variable 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

...or ss2<- "f(5, a=3, b=4, c=\"1:4\", d=2)*z" 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, Statist...

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 > >...

...ial, data = dati) Coefficients: (Intercept) x -348.23 11.23 Degrees of Freedom: 12 Total (i.e. Null); 11 Residual Null Deviance: 17.94 Residual Deviance: 7.011e-10 AIC: 4 I don't understand the meaning of warning. Can anyone help me? Many thanks. Cheers, Vito ===== Diventare costruttori di soluzioni Became solutions' constructors "The business of the statistician is to catalyze the scientific learning process." George E. P. Box Top 10 reasons to become a Statistician 1. Deviation is considered normal 2. We feel complete a...

...l, "is.loaded("equil_distC")" would still lead to a TRUE. I would like to understand the structural reason for these problems. Again, the problem is not the C++ code - it does run correctly now, although it had crashed in the previous run. Best regards Laurent Franckx, PhD VITO NV Boeretang 200, 2400 MOL, Belgium Tel. + 32 14 33 58 22 Skype: laurent.franckx laurent.franckx at vito.be ________________________________ [http://www.vito.be/e-maildisclaimer/bsds.png] <http://bsds.vito.be> VITO Disclaimer: http://www.vito.be/e-maildisclaimer

...cept) anno -427017.1740 217.0588 [[2]] (Intercept) anno -39625.82146 21.78025 ..... [[12]] (Intercept) anno 257605.0343 -129.7646 I want create a data frame with two columns (intercept and anno)using data in these list. Any help will be appreciated. Best Vito ===== Diventare costruttori di soluzioni "The business of the statistician is to catalyze the scientific learning process." George E. P. Box Visitate il portale http://www.modugno.it/ e in particolare la sezione su Palese http://www.modugno.it/archivio/cat_palese.shtml