Re: [R] SNG uz malo regresije

2016-01-05 Thread Petar Milin
> On Jan 5, 2016, at 12:01 PM, Bojana Dinic wrote: > > najpre, srecna Nova godina i sve najbolje tbi i tvojoj porodici. Nadam se da > ces se za praznike lepo odmoriti. > Mene su vec zadesile neke obaveze, pa sam htela da te pitam da li znas neki > paket, makro,

[R] Hierarchical Cluster Analysis with large dataset

2013-11-03 Thread Petar Milin
Hello! Can anyone give me advice on running Hierarchical Cluster Analysis on large datasets? For example, 8x1. Calculating distances on such a dataframe seems impossible even on very powerful computer. Also, any other advice that would lead to reduction of dimensionality, i.e.,

[R] How to save very large matrix?

2013-10-29 Thread Petar Milin
Hello! I have a very large matrix of results: 5x10. I saved it as RDS, but I would also need to save it as txt or csv. Is there a way to do it? Now, with write.table I am receiving an error: Error in .External2(C_writetable, x, file, nrow(x), p, rnames, sep, eol, : long vectors not

Re: [R] How to save very large matrix?

2013-10-29 Thread Petar Milin
Hello, On Oct 29, 2013, at 10:16 PM, Prof Brian Ripley rip...@stats.ox.ac.uk wrote: On 29/10/2013 20:42, Rui Barradas wrote: Hello, You can use the argument to write.csv or write.table append = TRUE to write the matrix in chunks. Something like the following. That was going to be my

[R] Getting Approximate Estimates after Deleting Factors out from fastbw()

2013-09-12 Thread Petar Milin
Hello! I am using relatively simple linear model. By applying fastbw() on ols() results from rms package I would like to get subtable Approximate Estimates after Deleting Factors. However, it seems this is not possible. Am I right? I can only get coefficients for variables kept in the model

[R] Number of pieces in segmented package

2013-01-29 Thread Petar Milin
Hello ALL, I just started to use segmented package for some piecewise linear fitting problem. However, author noted: The package is not concerned with estimation of the number of the breakpoints. Although the BIC has been suggested, in general nonstatistical issues related to the understanding

Re: [R] Parametric effects in GAM

2012-10-03 Thread Petar Milin
that particular effect. Then, the same for the second interaction numericB:factorA. Unfortunately, I am not sure that my thoughts would pass sanity check. Can you comment, please? Many thanks again! Best, Petar On 02/10/12 20:36, Petar Milin wrote: Hello! Can anyone give a tip how to plot parametric

[R] Parametric effects in GAM

2012-10-02 Thread Petar Milin
Hello! Can anyone give a tip how to plot parametric effects in an Generalized Additive Model, from mgcv package? Thanks, PM __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide

[R] Printing localized pdf file from Sweave snippets

2012-09-14 Thread Petar Milin
Hello! I am struggling for quite some time with proper printing of local characters in pdf plot, via Sweave snippets in Rnw file. When I am working directly within R, all is fine and I can get local character properly, like: pdf('figs/fig-relativeEntropy0.pdf', h=6, w=6, encoding='CP1250')

Re: [R] Printing localized pdf file from Sweave snippets

2012-09-14 Thread Petar Milin
\u010dina (open-sea)' Duncan Murdoch On Fri, 2012-09-14 at 08:17 -0400, Duncan Murdoch wrote: On 14/09/2012 7:07 AM, Petar Milin wrote: Hello! I am struggling for quite some time with proper printing of local characters in pdf plot, via Sweave snippets in Rnw file. When I am working

[R] Which BUGS should one use?

2012-08-30 Thread Petar Milin
Hello ALL! Some times ago I started to learn and play with Bayesian stuffs. Many advice use of WinBUGS for Bayesian inference Using Gibbs Sampler. However, WinBUGS is discontinued, and now, development is under OpenBUGS. I wasn't lazy, so I installed both and tried out. In more than 90% of cases

[R] Residuals from binomial glm

2012-05-13 Thread Petar Milin
Dear all! Sometimes, when using plain lm(), if I have highly correlated numeric predictors (covariates), I firstly ran lm() to get residuals of one from the other. This way, they became uncorrelated (orthogonal), and I can check if the residualized one contributes to prediction, over and above the

[R] Installing texlive dependencies

2012-04-18 Thread Petar Milin
Hello ALL! I am running Fedora 16 x86_64. Due to some dubious problems with the TeXLive (2007, which is a default), I removed it. That removed R as well, and some other dependent packages. Then, I installed TeXLive 2011 from CTAN. However, when I wanted to install R, from Fedora's repositories, it

Re: [R] Installing texlive dependencies

2012-04-18 Thread Petar Milin
, but placing a symbolic link at /usr/share/texmf, pointing at /usr/local/texlive/2011 might do the trick. It's a bit of a hack, but that often works for me. Best, Randy On Apr 18, 2012, at 7:31 AM, Petar Milin wrote: Hello ALL! I am running Fedora 16 x86_64. Due to some dubious

[R] Are they fully identical: WinBUGS and OpenBUGS; R2WinBUGS and R2OpenBUGS

2011-10-20 Thread Petar Milin
Hello ALL! I am running Linux, Fedora 15 64-bits, and R on it. I need to use WinBUGS and R2WinBUGS, but as far as I read, WinBUGS is closed project, to be continued with/as OpenBUGS. Thus, I have found R2OpenBUGS on OpenBUGS Contributed Code (http://openbugs.info/w/UserContributedCode), not on

[R] Strange pdf() in Fedora 15

2011-09-15 Thread Petar Milin
Hello useRs! Recently, I migrated from Debian i386 to Fedora 15 64-bit, after a long and happy experiences with Debian distributions (first Ubuntu, later Debian itself). I installed R with yum install R, and then necessary packages. Everything seems to work fine, for now, except that exporting

[R] Checking and building package

2011-06-03 Thread Petar Milin
be in src/. What should be done if one wants to prepare for CRAN? Thanks! Best, Petar -- # Petar Milin Department of Psychology, University of Novi Sad, Serbia Laboratory for Experimental Psychology, University

Re: [R] Checking and building package

2011-06-03 Thread Petar Milin
! Petar On 03/06/11 22:32, Uwe Ligges wrote: On 03.06.2011 21:46, Petar Milin wrote: Hello! I am truing to compile an R-package having c-code. I put foo.c in src/ folder and useDynLib(foo) Where foo is the name of your package, I hope. Does R CMD INSTALL yourpackage generate a packagename.so

[R] Compiling C-code in Windows

2011-05-31 Thread Petar Milin
Hello ALL! I am an Linux user (Debian testing i386), with very dusty Win-experience. Nevertheless, my colleagues and I are making some package in R, and I built C-routines to speed up things. I followed instruction how to compile C for R (very useful link:

[R] Number of dimension in Multidimensional Scaling

2010-12-09 Thread Petar Milin
Hello! Very often one can hear that MDS usually ends with two-dimensional solution. Of course, there are methods, like Scree-test (proposed by Kruskal and Wish, 1981), to determine optimal number of dimensions. However, I am trying to find references to this two-dimensional gold-standard. Can

[R] How to obtain prediction intervals for random effects?

2010-11-12 Thread Petar Milin
Hello! I wonder how can one get upper and lower limits of a prediction interval -- exact values? They are shown on caterpillar plot using ranef() with argument postVar=TRUE, but I would like to know them. A while ago, some discussions were opened on Confidence Intervals for Random Effect

[R] What solve() does?

2010-09-01 Thread Petar Milin
Hello! Can anyone explain me what solve() function does: Gaussian elimination or iterative, numeric solve? In addition, I would need both the Gaussian elimination and iterative solution for the course. Are the two built in R? Thanks! PM __

[R] Coinertia randtest

2010-08-23 Thread Petar Milin
Hello! I dunno why, but I cannot make randtes.coinertia() from ade4 package working. I have two nice distance matrices (Euclidean): dist1 1 2 3 4 5 6 7 2 2.5776799 3 1.7892825 1.0637487 4 1.0557991 2.4270728 2.0626604 5 1.6745483

[R] Coinertia randtest

2010-08-18 Thread Petar Milin
Hello! I dunno why, but I cannot make randtes.coinertia() from ade4 package working. I have two nice distance matrices (Euclidean): dist1 1 2 3 4 5 6 7 2 2.5776799 3 1.7892825 1.0637487 4 1.0557991 2.4270728 2.0626604 5 1.6745483

[R] Finding a scalar value...

2010-08-16 Thread Petar Milin
Hello! I need to find a simple scalar value: Scal = ((1/R) - T) / (P - T), where R, T, and P are vectors in a data.frame. Please, can anyone tell me how to solve that in R? Best, PM __ R-help@r-project.org mailing list

Re: [R] Finding a scalar value...

2010-08-16 Thread Petar Milin
11:48, Petar Milin wrote: Hello! I need to find a simple scalar value: Scal = ((1/R) - T) / (P - T), where R, T, and P are vectors in a data.frame. Please, can anyone tell me how to solve that in R? Best, PM __ R-help@r-project.org mailing list https

Re: [R] QQ plotting of various distributions...

2009-09-27 Thread Petar Milin
(test_data)) abline(0,1) It is helpful to assess the sampling variability, by creating repeated sets of test_data, and plotting all of these along with your observations to create a confidence envelope. The SuppDists provides Inverse Gauss. On Thu, Sep 17, 2009 at 11:46 AM, Petar Milin pmi

Re: [R] QQ plotting of various distributions...

2009-09-27 Thread Petar Milin
Thanks all! I did not want to cause any trouble and, God forbid, offense. I thought, I asked a simple question to improve my understanding and R-skills. It seems that there ain't single gospel truth about QQs. :-) Thanks, again! Best, PM Juliet Hannah wrote: I think it's helpful to show the

[R] QQ plotting of various distributions...

2009-09-17 Thread Petar Milin
Hello! I am trying with this question again: I would like to test few distributional assumptions for some behavioral response data. There are few theories about true distribution of those data, like: normal, lognormal, gamma, ex-Gaussian (exponential-Gaussian), Wald (inverse Gaussian) etc. The

Re: [R] Very basic question regarding plot.Design...

2009-09-09 Thread Petar Milin
(ols2), above), error message returns: Error in predictDesign(fit, adj, type = x, non.slopes = non.slopes) could not find function Varcov Best, PM Frank E Harrell Jr wrote: Petar Milin wrote: I would like to have a line on this plot, instead of two points: x1 = rnorm(100, 10, 2.5) x2

Re: [R] Very basic question regarding plot.Design...

2009-09-09 Thread Petar Milin
Sorry, I hope this will be the last (for now, at least). Following your advices, I did: n - 100 anxiety - c(rnorm(n, 10, 2.5), rnorm(n, 26, 3.2), rnorm(100, 25, 3.1), rnorm(100, 10, 2.6)) m.status - c(rep('married', n*2), rep('not married', n*2)) gender - c(rep('male', n), rep('female',

[R] Very basic question regarding plot.Design...

2009-09-08 Thread Petar Milin
Hello ALL! I have a problem to plot factor (lets say gender) as a line, or at least both line and point, from ols model: ols1 - ols(Y ~ gender, data=dat, x=T, y=T) plot(ols1, gender=NA, xlab=gender, ylab=Y, ylim=c(5,30), conf.int=FALSE) If I convert gender into discrete numeric predictor, and

[R] Applying qqmath using gamma distribution...

2009-09-04 Thread Petar Milin
Hello ALL! Can anyone tell me how to specify qqmath() function with distribution qgamma? In help it is bit vague how to pass shape and scale parameters for gamma. Thank you in advance. Best, PM __ R-help@r-project.org mailing list

Re: [R] Applying qqmath using gamma distribution...

2009-09-04 Thread Petar Milin
)) Is that correct? (gammafit() from mhsmm package calculates parameters of gamma function) Thank you, again. Best, PM Duncan Murdoch wrote: On 9/4/2009 8:36 AM, Petar Milin wrote: Hello ALL! Can anyone tell me how to specify qqmath() function with distribution qgamma? In help it is bit vague how

[R] Updated R on Debian testing machine...

2009-03-18 Thread Petar Milin
Hello! I recently switched to Debian testing OS and explanation at the: http://cran.at.r-project.org/bin/linux/debian/ is a little bit fuzzy. Can anyone give step-by-step how-to have updated R on Debian testing and/or unstable. Thanks in advance! PM __

[R] Warning messages about replacing previous import...

2008-02-21 Thread Petar Milin
Hello! I have posted this message before, but no answer. Now, after a few upgrades, I still have the same warning messages. I am a bit worried, since I do not know does this affect or could affect R. I use R under Linux (Ubuntu 6.06). After start, and after loading few libraries, I got

[R] Two warning messages after loading package 'coda'...

2008-01-30 Thread Petar Milin
Hello! Recently I updated my R and all packages, under Linux (Ubuntu 6.06). After that, I started to receive warning messages: 1: In namespaceImportFrom(self, asNamespace(ns)) : replacing previous import: cov2cor 2: In namespaceImportFrom(self, asNamespace(ns)) : replacing previous import: