[R] clusterCrit package produces Nan

2014-01-29 Thread Paola Tellaroli
38.57547 6.877532 18.36552 < datissimi = som(test, grid=somgrid(xdim=2, ydim=4, topo='rectangular'), rlen=500, keep.data=T) > somissimo = datissimi$unit.classif > intCriteria(test, somissimo, "Silhouette")

[R] reference for variance in svykm

2013-05-28 Thread Paola Rebora
and I was wondering if it was ever fomally validated and published. Thanks, Paola [[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 ht

[R] error in comparisons using glht function

2012-11-05 Thread paola
hi everybody this is my last chance, I dont know what else to do I was trying to find my multiple comparisons using the code in R glht next is my function and the error I get vol.interallaovsol<-aov(scoresol~ffjob+ffage+ffstudies+ffjob:ffstudies+ffjob:ffage+ffstudies:ffage+ffjob:ffage:ffstudies

[R] Contrasts in manova

2012-11-03 Thread paola
ind anything, maybe I didnt use in a corect way, *any advice is welcome* I always get the error Error in solve.default(wcrossprod(model.matrix(model), w = wts)) : Lapack routine dgesv: system is exactly singular > Thanks in advance Paola -- View this message

[R] error in lm

2012-10-30 Thread paola
(fit, scope$add, scale = scale, trace = max(0, trace - : no addterm method implemented for "mlm" models Does anybody knows why I get this error? Thanks in advance Paola -- View this message in context: http://r.789695.n4.nabble.com/error-in-lm-tp4647840.html Sent from the R hel

Re: [R] Change color in forest.rma (metafor)

2011-08-29 Thread Paola Tellaroli
Thank you so much!!! Could you tell me also how to change the size of the chart? There is not enough space below the chart to add the arrows! 2011/8/28 Uwe Ligges-3 [via R] < ml-node+3774557-1567708350-262...@n4.nabble.com> > > > On 26.08.2011 15:50, Paola Tellaroli wrote: >

Re: [R] Change color in forest.rma (metafor)

2011-08-26 Thread Paola Tellaroli
I lied, that was not my last question: how can I add two arrows at the bottom with the words "in favor of A / B"? This is not specified in the pdf and with "text" I have the impression that I can't add text below the x-axis. 2011/8/26 Paola Tellaroli > Dear Prof.

Re: [R] Change color in forest.rma (metafor)

2011-08-26 Thread Paola Tellaroli
have one last question. I'm trying to copy the style of a forest plot that I've seen and I like (the one in the attached file, page 1034): can I do this in R? Best wishes, *Paola* 2011/8/25 Viechtbauer Wolfgang (STAT)-2 [via R] < ml-node+3768683-1225159815-262...@n4.nabble.com>

Re: [R] Change color in forest.rma (metafor)

2011-08-25 Thread Paola Tellaroli
Thank you for your attention and help! In this way I get the diamond coloured, but actually I would have the squares representing the values of the individual studies coloured. Is it somehow possible? *Paola* 2011/8/24 Viechtbauer Wolfgang (STAT) < wolfgang.viechtba...@maastrichtuniversity

[R] Change color in forest.rma (metafor)

2011-08-24 Thread Paola Tellaroli
"Events", "N"), cex=1.2 ) text(c(-8.7, -5.5, 8), 6.5, c("S", "A", "Log")) text(-15, 6, "Trials", pos=4) text(10, 6, "Hazard Ratio [95% CI]", pos=2) par(op) Even if I have specified "col="red", border=&qu

[R] help with nls fitting

2011-08-20 Thread Paola Lecca
uced when evaluating the model What could be the error? I tried with differnet intial values of v4_Kd, but I did not fix the problem. Here below the code. Thanks in advance, Paola. Data_pp3_mrna <- data.frame( p1 = protein_1,

[R] problemsn in using nls

2011-08-05 Thread Paola Lecca
Dear all, I tried to use nls, but I got the following error Error in numericDeriv(form[[3L]], names(ind), env) : Missing value or an infinity produced when evaluating the model Any suggestion? Thanks, Paola. The code I wrote is Data_pp2_mrna <- data.fr

[R] fit a 2-variables function to data

2011-08-05 Thread Paola Lecca
Dearl all, I have to fit a function y = f(x1, x2) to data experiemntal data describing the measured behavior of y. x1 and x2 are the independent variables. Could you suggest me wich R package can I use for this purpose? Thanks, Paola. -- *Paola Lecca, PhD* *The Microsoft Research

[R] use of modMCMC

2011-08-04 Thread Paola Lecca
of parameters. modMCMC performs a Markov Chain Monte Carlo simulation. I do not understand very well how modMCMC can be used in a context of parameter estimation. Could someone help me in understanding the use of this function and its utility for parameter fitting? Thank you very much in advance, Paol

[R] use of modMCMC

2011-08-04 Thread Paola Lecca
MC performs a Markov Chain Monte Carlo simulation. I do not understand very well how modMCMC can be used in a context of parameter estimation. Could someone help me in understanding the use of this function and its utility for parameter fiting? Thank you very much in advance, Paola. -- *Paola Lec

[R] Help with modFit of FME package 2

2011-08-01 Thread Paola Lecca
ses I gave to them. The problem is not due to the fact that I have only one equation (I tried also with more equations, but I still obtain this error). I would appreciate if someone could help me in understanding the reason of the error and in fixing it. Thanks for your attention, Paola Lecca.

[R] Help with modFit of FME package

2011-08-01 Thread Paola Lecca
ntion, Paola Lecca. Here the error: > summary(Fit) Parameters: Estimate Std. Error t value Pr(>|t|) pro1_strength1 NA NA NA Residual standard error: 2.124 on 10 degrees of freedom Error in cov2cor(x$cov.unscaled) : 'V' is not a square numer

[R] Help with modFit of FME package

2011-07-28 Thread Paola Lecca
your attention, Paola Lecca. Here the error: > summary(Fit) Parameters: Estimate Std. Error t value Pr(>|t|) pro1_strength1 NA NA NA Residual standard error: 2.124 on 10 degrees of freedom Error in cov2cor(x$cov.unscaled) : 'V' is not a squ

[R] Help with modFit of FME package

2011-07-28 Thread Paola Lecca
attention, Paola. Here the error: > summary(Fit) Parameters: Estimate Std. Error t value Pr(>|t|) pro1_strength1 NA NA NA Residual standard error: 2.124 on 10 degrees of freedom Error in cov2cor(x$cov.unscaled) : 'V' is not

[R] frequency table help

2010-01-25 Thread Paola Pesci
is is not what I need. Do you know how I can do? Thank you in advance. Paola [[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://