differently. It looks like Pop uses apt so you
might try apt-get for R and see if that works or consider compiling
from source.
Collin Lynch.
On Mon, Oct 14, 2019 at 7:36 AM Dirk Eddelbuettel wrote:
>
>
> Adam,
>
> You may find this blog post and the video instruct
Magi, is there a reason that you need to run the script via R? If your
plan is to download the data via python than then process with R, you
might consider using the Rpy2 package to link them. This would allow
you to call the downloading code from python and then have python feed
the data to R.
kland.ac.nz
> > Sent: Mon, 27 Jul 2015 15:12:43 +1200
> > To: cfly...@ncsu.edu
> > Subject: Re: [R] VIF threshold implying multicollinearity
> >
> >
> > On 27/07/15 13:36, Collin Lynch wrote:
> >
> >> The following sources discuss the issues generall
The following sources discuss the issues generally and may be a goof
pointer to the literature on VIF. Particularly the Schroeder paper.
@article{Yi:Evaluation,
AUTHOR = {Youjae Yi},
TITLE = {On the Evaluation of Main Effects in Multiplicative
Regression Models.},
JOURNAL =
Hi Ignacio, If I am reading your code correctly then the top while loop is
essentially seeking to select a random set of names from the original set,
then using unique to reduce it down, you then iterate until you have built
your quota. Ultimately this results in a very inefficient attempt at
samp
right although the exact function of your fitness operator is not
clear to me. I recommend looking at "An Introduction to Genetic
Algorithms" by Melanie Mitchell. That is a good resource for more
general GA advice.
Sincerely,
Collin Lynch.
On Fri, May 22, 2015 at 3:09 AM, R
Hi Paul a quick search popped up these:
http://astrostatistics.psu.edu/datasets/R/html/index.html
http://finzi.psych.upenn.edu/
http://r.789695.n4.nabble.com/Online-R-documentation-td1009656.html
Are they what you are looking for?
Collin.
On Thu, Apr 16, 2015 at 5:02 PM, paul wrote:
> The
))
> hist(out2)
> mean( out2 <= 0.05 )
> binom.test( sum(out2 <= 0.05), B, p=0.05 )
>
> This simulates from a continuous exponential (skewed) and shifts to
> get the means (shifted location is a common assumption, though not
> required for the actual test).
>
> On T
Live: OO#.. Dead: OO#.. Playing
>> Research Engineer (Solar/BatteriesO.O#. #.O#. with
>> /Software/Embedded Controllers) .OO#. .OO#.
>> rocks...1k
>>
>> --
Greetings, I am working on a project where we are applying the
Kruskal-Wallace test to some factor data to evaluate their correlation with
existing grade data. I know that the grade data is nonnormal therefore we
cannot rely on ANOVA or a similar parametric test. What I would like to
find is a me
Greetings, I am working on a project where we are applying the
Kruskal-Wallace test to some factor data to evaluate their correlation with
existing grade data. I know that the grade data is nonnormal therefore we
cannot rely on ANOVA or a similar parametric test. What I would like to
find is a me
I recommend rpy2. http://rpy.sourceforge.net/rpy2.html
It provides direct access to a running R instance with full support for R
functions including package loading. It has some minor issues with
graphics drivers making it best for programmatic and not interactive use
but it is excellent for mun
a desire to retain that information.
> >
> > One way (not extracting cyrillic text):
> >
> > library(XLConnect)
> > DF <- readWorksheetFromFile( "exampX.xlsx", sheet="examp" )
> > pattern <- "^.*(\\d+) *\\* *(\\d+)[^\\d]*(\\d+) *\\*
Dr. Polanski, I would recommend something else. Given the messy nature of
your data I would suggest using a language like Python or Perl to extract
it to an appropriate format. Python has good regular expression support
and unicode support. If you can save your data as a csv file or even text
li
You might consider using something other than R to clean the file and even
to load it. I regularly use python to preprocess data for R and often feed
it to R directly via the rpy2 interface. If the dates are delimited by
some feature (e.g. ") you could potentially use the python csv library to
lo
If you do not need a pure R solution, you might also find it helpful to
blend languages. For scraping and munging tasks such as this I generally
turn to python to do extraction then feed data to R for analysis via rpy.
On Thu, Nov 20, 2014 at 8:57 PM, Spencer Graves <
spencer.gra...@structuremoni
Hi Katharina, what gam package are you using? With mgcv you can inspect
the results of the output variables to check whether the fixed field is
true which would indicate whether the df is fixed or floating. I'm not
sure if this is applicable to what you want.
My rough translation of your German
I'll echo this and expand. Hui it is possible, indeed likely, that you
are running a 32bit version of Java. In that case the error may be
attributed to a miscommunication between the two.
You can check you java version by running "java -version" in the command
prompt. If it is 64 bit it will te
Greetings, I'm working on some analyses where I need to calculate wilcox
tests for paired samples. In my current literature search I've found a
few papers on sample size determination for the wilcox test notably:
Sample Size Determination for Some Common Nonparametric Tests
Gottfried E. Noether
J
Hi Mike, I recently had this issue and didn't find any package that
implemented these tests directly for the gam object. I found it simplest
just to pull the residuals from it and run tests like shapiro.test
directly.
Best,
Collin.
On Thu, 5 Dec 2013, Mike.lang wrote:
> Dear all
Erin, at first glance I would say that this is an R error. When Rpy2
detects an error it will pass it through errors as library errors like
this. At first glance it appears that your r code is not loading the
requisite packages as it cannot find them. That I suspect is what is
causing your coord
> You can use the "offset" function as part of a formula in "lm" (and
> other model fitting functions) to set a specific slope or set of
> slopes. Using this up front will give you the correct residuals,
> standard errors, etc. This is better than trying to modify a fitted
> regression object.
G
> Even if this is possible, won't all the other estimates (i.e., standard error
> of betas) produced be junk since they aren't derived from the associated
> estimators?
That was actually my primary concern. It looks like the offset is the
solution.
Thanks.
Collin.
>
> Michael
>
Dila, take a look at this:
http://r.789695.n4.nabble.com/How-to-replace-all-lt-NA-gt-values-in-a-data-frame-with-another-not-0-value-td2125458.html
Does that help?
Best,
Collin.
On Tue, 12 Nov 2013, dila radi wrote:
> Hi all,
>
> I have a data set with missing value. I w
You might need to implement it as a nested pair of for loops using rbind.
In essence iterate over the rows in df1 and each time find the matching
row in df2. If none is found then add the df1 row by itself to the
result. If one is then remove it from df2 and rbind both of them. Once
done just me
> The default functional link for mgcv::gam is "log", so I doubt that
your theoretical understanding applies to GAM's in general. When Simon
Wood wrote his book on GAMs his first chapter was on linear models, his
second chapter was on generalized lienar models at which point he had
written over
falls under any analogue to the Gauss-Markov Theorem.
Thank you in advance for any help.
Sincrely,
Collin Lynch.
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R
hey seem to be in order.
>
> Thanks,
> Erin
>
>
>
> On Wed, Oct 30, 2013 at 10:51 PM, Collin Lynch wrote:
>
> > Erin can you share the internal error details?
> >
> > As a first guess are the files executable by all? CGI requires world rwx.
> >
>
Erin can you share the internal error details?
As a first guess are the files executable by all? CGI requires world rwx.
Best,
Collin.
On Wed, 30 Oct 2013, Erin Hodgess wrote:
> Hi again.
>
> I'm putting together a little project with R, python, and a website. So I
> have an H
; File "", line 1, in
> File "/usr/local/lib/python2.7/site-packages/rpy2/robjects/__init__.py",
> line 216, in __getitem__
> res = _globalenv.get(item)
> LookupError: 'buzz' not found
> >>>
> root@erinminfo [/home/erinminf/public_ht
Erin, one question, can you access the defined functions by key?
In lieu of:
> x = R.r.buzz(3)
Can you do:
x = R.r['buzz'](3)
Alternatively if you need only one or two custom functions have you
considered just defining them via python as in:
PStr = """
function(LM) {
S <- summary(LM);
pr
uot;)$loglik
[1] 11.15125
fitdistr(E + 1, "lognormal")$loglik
[1] -0.8575117
fitdistr(E, "exponential")$loglik
[1] -73.18107
Is this a sign of error, in my system or on my part? And if not am I
correct in interpreting lognormal as the best choice?
r the mgcv GAMs, can
anyone recommend a good test implementation? I am familiar with the
ncvTest method implemented in the car package but that applies only to
lms.
Thank you,
Collin Lynch.
__
R-help@r-project.org mailing list
Thank you Arun!
Collin.
On Mon, 7 May 2012, arun wrote:
> Hi Collin,
Look in the package 'pwr' for 'pwr.r.test'.
A.K.
- Original Message -
From: Collin Lynch
To: r-help@r-project.org
Cc:
Sent: Monday, May 7, 2012 1:44 AM
Subject: [R] Statistical pow
Great thanks Peter!
Collin.
On Mon, 7 May 2012, peter dalgaard wrote:
>
> On May 7, 2012, at 07:44 , Collin Lynch wrote:
>
> > My apologies for the statistical naivete of my question but...
> >
> > Is there an established method or calulating the statistical
My apologies for the statistical naivete of my question but...
Is there an established method or calulating the statistical power of a
correlation test? And if so is there a method in R for it?
Thank you,
Collin Lynch.
__
R-help@r
Am 23.09.2008 um 23:57 schrieb Peter Dalgaard:
> For this kind of problem I'd go directly for the binomial
> distribution. If the actual probability is 0, this is essentially
> deterministic and you can look at
>
> > binom.test(0,99,p=.03, alt="less")
>
> This means that you don't sample from the
Thank you Peter. That is incredibly helpfyul, and much much smaller!
Best,
Collin.
__
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
the specified significance levels.
My questions are those above, and, as a followup whether the caveat about
n being the number in each group means that I need to sample twice that
number in a single group. I don't believe so but I want to be sure.
Thanks in advance,
39 matches
Mail list logo