[R] Matching Values issue

2019-10-25 Thread Roberto Martinez Caballero
__ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible cod

Re: [R] Different stack barplots - same color legends

2018-11-04 Thread P. Roberto Bakker
Dear Richard, I ran your syntax and the plots look great and cool. Again, thank you so much for helping me out a lot and for all the programming. And again, Jim en David, thank you also so much. Now, I can continue with my work and I learned a lot from you. Best, Roberto Op za 3 nov. 2018 om 14

Re: [R] Different stack barplots - same color legends

2018-11-03 Thread P. Roberto Bakker
Thank you Op do 1 nov. 2018 om 00:25 schreef Jim Lemon : > Hi Roberto, > Here is a snippet of code that translates the text responses of the > BIS-11 into numeric values. Note the reversal of the order in the > second item: > > BIS$Q1<-as.numeric(factor(BIS$Q1, >

Re: [R] Different stack barplots - same color legends

2018-11-03 Thread P. Roberto Bakker
Dear Richard, Thank you so much for all your work and time you punt in it. I will start with your suggestions and let you know how far I come. Also thanks to the others who helpt me. Best Roberto Op do 1 nov. 2018 om 23:42 schreef Richard M. Heiberger : > ## reminder on how the lev

Re: [R] Different stack barplots - same color legends

2018-11-01 Thread P. Roberto Bakker
> > On Thu, Nov 1, 2018 at 9:01 PM P. Roberto Bakker > wrote: > > > > Hi Jim, > > > > Thank you for your quick reply. It is a great procedure. > > The response options in my data.frame are (fortunately) similar in all > columns. > > It would be nice

Re: [R] Different stack barplots - same color legends

2018-11-01 Thread P. Roberto Bakker
names. This means that I need to put each (long) column name after "$: dataname$"very long name". So, is ther a way to do this procedure for all columns? Roberto Op do 1 nov. 2018 om 10:50 schreef Jim Lemon : > Hi Roberto, > What I suggested is a brute force method of

Re: [R] Different stack barplots - same color legends

2018-11-01 Thread P. Roberto Bakker
Hi Jim, Thank you. An additional question: as I have many columns to change in numeric, and the columns are long sentences, what is an efficient way to do this? I checked in StackOverflow but could not find the right answer Best Roberto Op do 1 nov. 2018 om 00:25 schreef Jim Lemon : >

Re: [R] Different stack barplots - same color legends

2018-10-31 Thread P. Roberto Bakker
might want to convert them to factors where you control the order > of the levels. > > Factor.wrong <- factor(c("mm", "cm", "m", "km")) > > levels(Factor.wrong) ## alphabetical order, not meaning order > [1] "cm" "km" &quo

Re: [R] Different stack barplots - same color legends

2018-10-31 Thread P. Roberto Bakker
rengroep", Op wo 31 okt. 2018 om 16:24 schreef Richard M. Heiberger : > part is fine. just be sure that the small part causes the problem. > I will need that to investigate what is happening. > > > On Wed, Oct 31, 2018 at 11:15 AM, P. Roberto Bakker > wrote: > > It is

Re: [R] Different stack barplots - same color legends

2018-10-31 Thread P. Roberto Bakker
It is a very long result text. I can send it to you, or is part of it ok?[ Op wo 31 okt. 2018 om 14:27 schreef Richard M. Heiberger : > Please send me the > dput(teamq) > > > On Wed, Oct 31, 2018 at 03:51 P. Roberto Bakker > wrote: > >> Thank you for you information.

Re: [R] Different stack barplots - same color legends

2018-10-31 Thread P. Roberto Bakker
m(teamq) [1] 4 2 > ncol(teamq) [1] 2 So it should be good. I used 'make.names' , in case the spaces in the variable names would be a problem. Same error. What could I do? Best and thank you in advance. Roberto Op ma 22 okt. 2018 om 20:10 schreef Richard M. Heiberger : > T

Re: [R] Different stack barplots - same color legends

2018-10-22 Thread P. Roberto Bakker
e situatie enigszins van toepassing is voor u of uw supervisorengroep", "de situatie in hoge mate van toepassing is voor u of uw supervisorengroep", "de situatie in zeer hoge mate van toepassing is voor u of uw supervisorengroep"),

[R] Different stack barplots - same color legends

2018-10-22 Thread P. Roberto Bakker
the same likert-item, which I understand - the colors go in sequence along the table. Question: how can I write a syntax that each likert-item has the same legend color? Thank you in advance, Roberto SYNTAX: counts19 <- table(teamq[,19]) counts20 <- table(teamq[,20]) barplot(as.matrix(counts1

Re: [R] Copy text from Script syntax into .txt

2018-04-25 Thread P. Roberto Bakker
Dear Don, Thank you for your information. Bw Roberto 2018-04-25 23:13 GMT+02:00 MacQueen, Don : > Try putting this >options(echo=TRUE) > at the beginning of your script > > See ?source for a clue > > -Don > > -- > Don MacQueen > Lawrence Livermore Nation

Re: [R] Copy text from Script syntax into .txt

2018-04-25 Thread P. Roberto Bakker
Hi Greg, Thank you for all information. Bw Roberto Op di 24 apr. 2018 21:42 schreef Greg Snow <538...@gmail.com>: > Look at the spin and stitch functions in the knitr package if you want > to process an existing script into an output that mixes the code run > with the outp

Re: [R] Copy text from Script syntax into .txt

2018-04-25 Thread P. Roberto Bakker
Hi Thierry, Thank you for your information. Bw Roberto Op di 24 apr. 2018 12:01 schreef Thierry Onkelinx : > Dear Roberto, > > The easiest way IMHO is to convert your script into an R markdown > document. See https://rmarkdown.rstudio.com/ > > Best regards, > &

[R] Copy text from Script syntax into .txt

2018-04-24 Thread P. Roberto Bakker
Hi everybody, How can I get text from RScript (e.g. syntax, reminder) into the result text. Sink() does not do that - I only read the results and therefore I have to 'guess' which syntax was used where - reminders I wrote are lost. Bw and thank you in advance, Roberto [[a

[R] Dataframe is character

2017-11-17 Thread P. Roberto Bakker
t I do not know why it did sometimes. Thank you in advance, Roberto PS I used "guess". The problem is not solve by using "text", "numeric" etc My syntax (I think I cannot send the excel file as binary?) > library(readxl) > library(readxl) > library(metafor) &g

[R] Problem with the port algorithm in nls function

2015-11-16 Thread roberto marrone
Dear all, using the following code I have find a problem with the port algorithm. If I use the nls function without lower bound for my parameter it compute the parameter's value but I want the parameter positive. Running the following code I had this error. Thanks at all. optim <- nls(Prezzo

[R] Problem with nls function

2015-11-16 Thread roberto marrone
Dear all, I have a problem using the R finction nls. I am trying to perform an optimisation of the volatility parameter in the Black and Scholes formula. In the function nls I wrote as a formula the call option price with the only unknown parameter the volatility that I called theta. The code

Re: [R] remove outliers with scatterplotMatrix

2015-03-19 Thread Pagliari, Roberto
Thank you John, I will try that one. Bob, From: John Fox [j...@mcmaster.ca] Sent: Thursday, March 19, 2015 2:08 PM To: Pagliari, Roberto Cc: r-help@r-project.org Subject: RE: [R] remove outliers with scatterplotMatrix Dear Roberto, This is, I assume

[R] remove outliers with scatterplotMatrix

2015-03-19 Thread Pagliari, Roberto
outliers may sometimes make a scatter plot less understandable. when using scatterplotMatrix, is it possible to automatically remove them? [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, se

[R] MICE skipping function

2013-02-05 Thread roberto
Hi Sir, I'm Roberto Palloni, student in economics. My question is about Multiple Imputation with MICE. I would like to know if skipping the imputation of a variable used as predictor could affect the estimation of the dependent variable. In other words, does the imputed values of a predicto

Re: [R] 3D Gaussian with different colors

2012-12-27 Thread Roberto
*cos(phi)), 1+(1*sin(phi)), fn(2, 1, 1, 0.5), pmat = res), col = "black", lwd = 1, lty = 2) polygon(trans3d(-3+(1*cos(phi)), 1+(0.5*sin(phi)), fn(-3, 1, 1, 0.5), pmat = res), col = rgb(0,0,1,0.5), border = NA) lines(trans3d(-3+(1*cos(phi)), 1+(0.5*sin(phi)), fn(-3, 1, 1, 0.5), pmat =

[R] 3D Gaussian with different colors

2012-12-26 Thread Roberto
I wrote, finding pieces of code in the web, a simple script to draw a 3D gaussian plot. For the next step, I need that the 2 gaussian distributions have different colors. Can someone help me to do this? fn <- function(x, y, scale, scale2) dnorm(x,mean=1,sd=scale)*dnorm(y,mean=-1,sd=scale) + dnorm

[R] Training, Validation and Test Sets and LDA with MASS package

2012-11-30 Thread Roberto
Hi all, please, can someone give me an example of LDA with MASS package with training, validation and test sets? I never used it, so I need and example to avoid to make errors. Thank you very much. Best, Roberto -- View this message in context: http://r.789695.n4.nabble.com/Training

Re: [R] Problem to read non-standard csv file

2012-10-12 Thread Roberto
.319 ... > # $ Tempo : num 0 0.002 0.004 0.104 0.204 0.304 > 0.404 0.504 0.604 0.704 ... > # $ Deformazione.in.compressione: num 0 0 0 0.0108 0.0245 ... > A.K. Thank you very much. It work very well! Best, Roberto -- View this message i

[R] Problem to read non-standard csv file

2012-10-12 Thread Roberto
ata_1.csv", header = TRUE, sep = ";", dec = ",", nrows = -1, encoding = "ANSI") The possibility to open the file into a text editor and replace comma with point is a crazy way, because I have around 1000 files to elaborate. The only way for me is to use R. I hope that so

[R] PLS-DA and LV variance with mixOmics

2012-09-27 Thread Roberto
Hi all, I need to obtain the LV variance from my PLS-DA analysis. I tried to read the reference manuale of the package, but I do not found information about that. Someone know a way to do it? Thank you, Roberto -- View this message in context: http://r.789695.n4.nabble.com/PLS-DA-and-LV

Re: [R] Package to remove collinear variables

2012-08-09 Thread Roberto Moscetti
Thank you very much for your help. Really appreciated. Best regards, Roberto Il 07/08/2012 12:18, John C Frain ha scritto: For background have a look at http://en.wikipedia.org/wiki/Multicollinearity. I have also used Regression Diagnostics: Identifying Influential Data and Sources of

Re: [R] Package to remove collinear variables

2012-08-05 Thread Roberto Moscetti
book (or literature) than can help me? Thank you very much for your time and suggestions. Best regards, Roberto Il 05/08/2012 12:55, Jeff Newmiller ha scritto: There is no "magic bullet" (package) for your problem. You must either learn enough statistics to understand how to analyze

Re: [R] Package to remove collinear variables

2012-08-05 Thread Roberto
. I tried to use VIF-Regression to select variables, but rfe function advise me with the same warning messages again. What do you think about that? Thank you very much for your help. Best, Roberto -- View this message in context: http://r.789695.n4.nabble.com/Package-to-remove-collinear

[R] Package to remove collinear variables

2012-08-04 Thread Roberto
Hi, I need to remove collinear variables to my Near-Infrared table of spectra. What package can I use? Something simple, because I am a novice about statistic. Thank you. Best regards, Roberto -- View this message in context: http://r.789695.n4.nabble.com/Package-to-remove-collinear

Re: [R] How to link two R packages together

2012-08-01 Thread Roberto
Yuo can use source function. Looking for ?source for an example. Best, Roberto -- View this message in context: http://r.789695.n4.nabble.com/How-to-link-two-R-packages-together-tp4638765p4638766.html Sent from the R help mailing list archive at Nabble.com

[R] Different results between lda(mass) and spss discriminant analysis

2012-08-01 Thread Roberto
ag(prop.table(ct, 1)) min.err <- 1-sum(diag(prop.table(ct))) fit_cv <- lda(samples ~., data = r, CV =T) ct_cv <- table(r[,1], fit_cv$class) errors_cv <- 1-diag(prop.table(ct_cv, 1)) min.err_cv <- 1-sum(diag(prop.table(ct_cv))) thank you for your help! Best regards, Roberto -- Vi

Re: [R] How to create a floating bar plot

2012-05-31 Thread Roberto Brunelli
rary(HH) > twobar <- data.frame(neg=c(8,0,9,2), pos=c(12,9,0,4)) > likert(twobar) > > Rich > > On Thu, May 31, 2012 at 9:30 AM, Sarah Goslee > wrote: >> >> Hi Roberto, >> >> The R Graph Gallery is an excellent resource for this kind of >> ques

[R] How to create a floating bar plot

2012-05-31 Thread Roberto Brunelli
sitive and negative parts of the bar should have a different color. Is there any function/package supporting this kind of plot? Thanks a lot, Roberto __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the pos

[R] Hyperspec package: need to change spectra names in a stacked plot

2012-04-27 Thread Roberto
regards, Roberto -- View this message in context: http://r.789695.n4.nabble.com/Hyperspec-package-need-to-change-spectra-names-in-a-stacked-plot-tp4593690p4593690.html Sent from the R help mailing list archive at Nabble.com. __ R-help@r-project.org

[R] error loading tcltk2

2012-04-20 Thread Roberto Ugoccioni
Hello, I just installed R-2.15.0 on windows XP and cannot load package tcltk2 (which I just downloaded from CRAN as tcltk2_1.2-1.zip; package install reported no problems): > library(tcltk2) Carico il pacchetto richiesto: tcltk Loading Tcl/Tk interface ... done Error : .onLoad failed in loadNamesp

Re: [R] ANOVA and LSD TEST - contrasting results

2011-11-20 Thread Roberto
I didn't meant to rewrite textbook. I asked only for a clarification to my doubt, because I have started to study statistic recently. -- View this message in context: http://r.789695.n4.nabble.com/ANOVA-and-LSD-TEST-contrasting-results-tp4086090p4088810.html Sent from the R help mailing list arch

Re: [R] ANOVA and LSD TEST - contrasting results

2011-11-19 Thread Roberto
First of all, thank you for the reply to my topic. Anova summary has shown no significative difference into Treatment (Pr 0.374), otherwise LSD Test has shown difference between 2 Treatment. I suppose that's something strange. Roberto -- View this message in context: http://r.7896

[R] ANOVA and LSD TEST - contrasting results

2011-11-19 Thread Roberto
0.02964944 ab D2 0.02913658 ab B2 0.02891764 b C2 0.0284452 how is it possible? my mistake to set anova? Best regards, Roberto -- View this message in context: http://r.789695.n4.nabble.com/ANOVA-and-LSD-TEST-contrasting-results-tp4086090p4086090.html Sent fro

[R] R-GoogleVis: I need to import my chart into powerpoint file

2011-10-08 Thread Roberto
"Simple" question: how can I do that? I think that I need to find a possible way to save my googlevis motion plot into a SWF file, but how? Thanks for the help, Roberto. -- View this message in context: http://r.789695.n4.nabble.com/R-GoogleVis-I-need-to-import-my-chart-into-power

[R] Problems to extract data from anova table

2011-08-06 Thread Roberto
Hi, this is my script > anova <- aov(data ~ Ts*Te*t + Error(R/Ts*Te*t)) > results <- summary(anova$Within) this is results Df Sum Sq Mean Sq F valuePr(>F) Ts 2 1232.2 616.11 53.606 3.965e-10 *** Ts:Te 4 4889.5 1222.37 106.356 4.075e-16 *** Ts:t 4 6472.1 16

[R] How does locpoly (KernSmooth package) estimate densities?

2011-07-16 Thread Roberto Molinari
! Roberto. [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html

[R] Tukey HSD with repeated measure ANOVA

2011-07-13 Thread Roberto
quot;) how can I do? my apologies for the banal questions, but I'm taking the first steps in statistical computing Thanks, Roberto -- View this message in context: http://r.789695.n4.nabble.com/Tukey-HSD-with-repeated-measure-ANOVA-tp3665872p3665872.html S

Re: [R] Computing row differences in new columns

2011-03-21 Thread Roberto Lodeiro Muller
-Original Message- From: Roberto Lodeiro Muller To: roberto.mul...@doctor.com Sent: Mon, Mar 21, 2011 3:37 pm Subject: Re: [R] Computing row differences in new columns Sorry, my data appeared badly formatted to me, so I converted it to plain text: And just to clarify, for each

[R] Computing row differences in new columns

2011-03-21 Thread Roberto Lodeiro Muller
ished that ? I guess questions like that have been already answered a thousand times, so I apologize for asking one more time. Thanks Roberto SUBJECT Date Result DateDiff ResultDiff 10751 22-Jul-03 3.5 10751 13-Feb-04 1.3 10751 20-Aug-04 1.6 10751 08-Mar-05 1.7 10751 30-A

[R] Fwd: How to get all combinations between two character vectors?

2011-03-11 Thread roberto . muller
Hi Rainer Or maybe you are referring to the outer function. I'm a newbie in R, but I recently read something about it in the pdf book named below (pages 91-2). I send to you an excerpt: An introduction to R An introduction to R Longhow Lam 6.2.5 The outer function outer function Th

Re: [R] waldtest and nested models - poolability (parameter stability)

2010-12-06 Thread Roberto Patuelli
onformable arguments Cheers Roberto - Original Message - From: "Achim Zeileis" To: "Patuelli Roberto" Cc: Sent: Monday, December 06, 2010 9:38 PM Subject: Re: [R] waldtest and nested models - poolability (parameter stability) On Mon, 6 Dec 2010, Roberto Patuelli wro

[R] waldtest and nested models - poolability (parameter stability)

2010-12-06 Thread Roberto Patuelli
lts are not reliable because of the overdispersion - basically without correcting for overdispersion almost every variable is highly significant anova(inv.log.leva.base, inv.log.leva, test = "Chisq") Thanks everyone! Best regards, Roberto Patuelli Roberto Pat

Re: [R] TimeStamp

2010-09-14 Thread Roberto Badilla Fuentes
Thanks. Both methods definitely help. -Roberto 2010/9/14 Uwe Ligges > > > On 06.09.2010 03:51, Roberto Badilla Fuentes wrote: > >> Hi, >> >> I have a dataset in .dbf format. It contains Coordinates and Time. >> The TIMESTAMP is as follows: >> >>

[R] TimeStamp

2010-09-05 Thread Roberto Badilla Fuentes
Anyone help me out why R does this and how I can correct it. Thanks -Roberto [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R

[R] anova for plm objects

2010-08-26 Thread Roberto Patuelli
'panelmodel')" waldtest > waldtest(depFE.plm, depLDFE.plm) Errore in modelCompare(objects[[i - 1]], objects[[i]], vfun = vcov.) : models are not nested Did anyone code a plm method for anova? OR Does anyone know why according to waldtest my models are not nested? Thanks Ro

Re: [R] Percentage effects in logistic regression

2009-11-09 Thread Roberto Patuelli
think the same concept is used in other fields (such as medicine) by the name "percentage effects". It should not be difficult to find but somehow I'm having problems in solving this small issue... Anyone has used this in R before? Thanks Roberto - Original Message -

Re: [R] Percentage effects in logistic regression

2009-11-09 Thread Roberto Patuelli
ed in the effect on y of a 1% increase in x (called percentage effects, right?). Could you please clarify? Thanks Roberto - Original Message - From: "Daniel Malter" To: "Patuelli Roberto" ; Sent: Monday, November 09, 2009 7:44 PM Subject: AW: [R] Percentage

[R] Percentage effects in logistic regression

2009-11-09 Thread Roberto Patuelli
asticity"). For example, if my independent variables are in logs, the betas can be directly interpreted as percentage effects both in OLS and Poisson regression. What about the logistic regression? Is there a package (maybe effects?) that can compute these automatically? Thanks and b

[R] Regular expression substitution ...

2009-11-02 Thread Roberto Brunelli
I would like to replace all occurences of /pc with \% using something like sub("/pc", "\%", x) but I am unable to make the '\' pass through. Any hint? Thanks, Roberto __ R-help@r-project.org mailing lis

Re: [R] loop vs. apply(): strange behavior with data frame?

2009-10-22 Thread Roberto Perdisci
Thanks for the suggestion. I found some documentation on why accessing a data.gram using the matrix notation (e.g., [i,j]) is so expensive, which was the cause of the problem. regards, Roberto On Thu, Oct 22, 2009 at 12:05 AM, Jim Holtman wrote: > try running Rprof on the two examples

[R] loop vs. apply(): strange behavior with data frame?

2009-10-21 Thread Roberto Perdisci
"tf = " user system elapsed 0.088 0.000 0.087 [1] "ta = " user system elapsed 0.128 0.000 0.128 [1] "--" [1] "Same experiment on data.frame..." [1] "diff = 0" [1] "tf = " user s

[R] spplot and cex

2009-09-21 Thread Roberto Brunelli
Dear community, I'm using R 2.6.1 on a GNU/Linux Slackware box. I would like to obtain larger text when using spplot (from sp 0.9-43) using spplot(... cex=2) but I get no visible effect. Any help? Thanks, __ R-help@r-project.org mailing list https:/

Re: [R] Winsorized mean and variance

2009-08-27 Thread Roberto Perdisci
This is of great help, thanks! Roberto On Thu, Aug 27, 2009 at 7:20 AM, Jim Lemon wrote: > Roberto Perdisci wrote: >> >> Hello everybody, >>  after searching around for quite some time, I haven't been able to >> find a package that provides a function to co

[R] Winsorized mean and variance

2009-08-26 Thread Roberto Perdisci
Hello everybody, after searching around for quite some time, I haven't been able to find a package that provides a function to compute the Windorized mean and variance. Also I haven't found a function that computes the trimmed variance. Is there any such package around? thank

[R] How do I retrieve column and row names after comparing two matrices?

2009-02-25 Thread Nussenzveig, Roberto H.
Hello, I have two matrices as shown below: Matrix 1 IDAB1BC1CD1... name1 1,1 2,1 0,2 ... name2 2,0 1,2 1,2 ... name3 0,2 1,1 2,0 ... name4 2,0 0,2 0,2

[R] Fitted values and residuals from glmmPQL (MASS package)

2009-01-30 Thread Roberto Patuelli
columns of fitted values. There's probably something about the way mixed models are treated in R that I'm not catching here... Can someone help me? Thanks and best regards, Roberto Patuelli ******** Roberto Patuelli, Ph.D. Post-doc researcher Institute for Economic Res

[R] Joint significance of more regressors in summary

2009-01-21 Thread Roberto Patuelli
ard error and a p-value for each set. But the degrees of freedom of the model would in reality be different... Thanks in advance for any help! Cheers Roberto Patuelli Post-doc researcher Institute for Economic Research (IRE) University of Lugano Email: roberto.patue...@lu.unisi.ch Homepage: h

[R] do.call and plotting functions ...

2008-11-19 Thread Roberto Brunelli
I'm trying to write a simple wrapper for plotting functions to make them print to postscript, something like ploteps <- function(file, plotFunction, ...) { args <- list(bquote(...)) # prepare postscript device do.call(plot, args) # close postscript device } I have inserted the bquote

[R] Marginal effects in negative binomial

2008-10-28 Thread Roberto Patuelli
compute them. Thanks in advance for your help. Roberto Patuelli Roberto Patuelli, Ph.D. Post-doc researcher Institute for Economic Research (IRE) University of Lugano (USI) via Maderno 24, CP 4361 CH-6904 Lugano Switzerland Phone: +41-(0)58-666-4166 Email: [EMAIL PROTECTED

[R] Problem with glm.nb estimation

2008-10-03 Thread Roberto Patuelli
stepwise procedure fails at the first step (I know step() is not made for the glm.nb() function, but it usually works just fine for me). Thanks in advance to everyone who might drop his/her two cents on what the problem is. Sincerely, Roberto **** Roberto Patuelli, Ph.D. Post-doc re

[R] Help please! How to code a mixed-model with 2 within-subject factors using lme or lmer?

2008-09-14 Thread roberto toro
example? Thanks again! roberto ps1. How would you code this with lmer()? ps2. this is part of the output of mod<-lme: > summary(mod) Linear mixed-effects model fit by REML Data: vlt AIC BIClogLik 57528.35 57639.98 -28745.17 Random effects: Formula: ~1 | Subjec

[R] Help please! How to code a mixed-model with 2 within-subject factors using lme or lmer?

2008-09-14 Thread roberto toro
|Subject),data.vslt) but I can't figure out the /(Lobe*Tissue) part... Thank you very much in advance! roberto __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.or

Re: [R] moving from aov() to lmer()

2008-09-13 Thread roberto toro
is straightforward mod2<-lmer(Volume~Sex*Lobe*Tissue+(1|Subject),data.vslt) but I can't figure out the /(Lobe*Tissue) part... Thank you very much in advance! roberto __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/list

[R] moving from aov() to lmer()

2008-09-13 Thread roberto toro
bject),data.vslt) but I can't figure out the /(Lobe*Tissue) part... Thank you very much in advance! roberto [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do r

[R] write dataframes

2008-09-09 Thread Roberto Olivares-Hernández
ave tried to generate objects like ob1<-object1$ORF ob2<-object2$ORF and then use cbind like, HG<-cbind(on1,ob2) but there is an error. Is there any other function I can use? Thanks for the help Roberto __ R-help@r-project.org mailing list https:/

Re: [R] Short courses for beginners?

2008-06-23 Thread Roberto Laforgia
See the following link: http://www.xlsolutions-corp.com/courselist.htm here you can find some courses in the U.S. Regards, R.L. David Hewitt wrote: > > Recently I've had a couple folks ask me whether there are short courses or > tutorials offered to help them get started in R. They're think

Re: [R] Book on R graphics

2008-06-14 Thread Roberto Laforgia
See the following link: http://www.amazon.com/Graphics-Computer-Science-Data-Analysis/dp/158488486X/ref=wl_itt_dp?ie=UTF8&coliid=I3TIYS7LRDPIED&colid=1ZQCB91VGW5UR Regards, R.L. Marcin Kozak wrote: > > Hi all, > > I am looking for a book from which one could learn a wide range of > graphics

[R] Some problem with the OLS

2008-06-11 Thread roberto laforgia
Dear all!! I estimated the following OLS model with R 2.5.6: Output R 2.5.6 Call: lm(formula = UN ~ log(x) + time2, data = dati) Residuals: Min 1Q Median 3QMax -5.649 -2.753 -1.015 1.225 16.199 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercep

[R] R-packages installation problems

2008-03-31 Thread Roberto Cremonini
Indeed R headers (R.h and Rinternals.h) don' t exist on my filesystem. Where can I find them? How can I install them? Thanks for help, Best regards, Roberto -- ARPA Piemonte SC05 - Area Previsione e Monitoraggio Ambientale Roberto Cremonini via Pio VII, 9 10135 TORINO ITALIA Tel. +

[R] approxfun + save: problem as R_approx changed package

2008-03-26 Thread Ugoccioni Roberto
"R_approx" not in DLL for package "stats" It seems to me the two versions of R are looking for the same thing in different places, is it because of the environment they carry? Can anyone suggest a way to work around this situation? Thanks, Roberto Roberto Ugoccioni I

[R] dimnames

2008-02-13 Thread Roberto Olivares Hernandez
#matrix dimnames(XX)<-NULL write.table(XX,"XX.txt",quote=FALSE,sep="\t") Thanks in advance Roberto [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-

[R] filter data from data frame

2008-02-12 Thread Roberto Olivares Hernandez
7C" "YAL012W" "YAL016W" and I want to generate a new data frame exclusively with the elements from the object X and its respective spectra value from the data frame Y. how to do it? I have tried to look into several functions, match, subset, command Y[ X,

[R] filter data from data frame

2008-02-12 Thread Roberto Olivares Hernandez
7C" "YAL012W" "YAL016W" and I want to generate a new data frame exclusively with the elements from the object X and its respective spectra value from the data frame Y. how to do it? I have tried to look into several functions, match, subset, command Y[ X,

[R] time-dependent ROC package

2008-02-07 Thread Roberto Grandi
Dear all, i am newbie to R software and i work with survival functions. Anyone knows a package to build and analyze time-dependent ROC curves? I need some help ;) Best regards Roberto __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman

[R] SQLiteDF SQLITE ERROR after attach

2007-12-07 Thread Roberto Iacopetti
eDF/iris.db", "iris.sdf") : .SQLiteDF/iris.db is not a valid SDF. perhaps anyone know thi problem thanks in advance Roberto Iacopetti -- View this message in context: http://www.nabble.com/SQLiteDF-SQLITE-ERROR-after-attach-tf4962680.html#a14214520 Sent from th

Re: [R] R2HTML how to pair graphic.png and table

2007-12-06 Thread Roberto Iacopetti
Eric, your code work well for my need, i'm not skillful in html environment, but your functions in R2HTML give simple many output requirements thanks Roberto Eric Lecoutre wrote: > > Hi Roberto, > > here is a way that presumes you know some (basic) HTML tags: &g

[R] R2HTML how to pair graphic.png and table

2007-12-06 Thread Roberto Iacopetti
co") browseURL(myfile) Thanks in advance Roberto Iacopetti -- View this message in context: http://www.nabble.com/R2HTML--how-to-pair-graphic.png-and-table-tf4956321.html#a14193218 Sent from the R help mailing list archive at Nabble.com. __ R-help@r

Re: [R] Help with K-means Clustering

2007-11-15 Thread Roberto Perdisci
to different final clusters (and therefore different centers and validity measures). see ?kmeans and look at the parameter 'centers'. Roberto http://roberto.perdisci.googlepages.com On Nov 14, 2007 6:07 PM, Alejandro Rodríguez <[EMAIL PROTECTED]> wrote: > > Hello, I'm

[R] one-class SVM in kernlab

2007-09-12 Thread Roberto Perdisci
Am I using ksvm() correctly? I also tried to use the formula ~. exluding the attribute 'type' from the dataset, and the results are exactly the same. Unfortunately reading ?ksvm didn't help me much. thank you, regards, Roberto __ R-hel