ption to this is Thomas
Lumley's biglm package, which processes the data in chunks. We need
more tools like these. Ultimately, you'll need to find some method of
analysis that is pretty smart memory-wise, and this may not be easy.
Best of luck,
Jay
-
Original
n the same computer). Both seem quite easy to use,
essentially only needing one command to initiate the "cluster" and
then one command to do something like apply() in parallel. It takes a
little planning of your application, but the "painfully obvious"
parallel problem should be pa
ize +
UpLow + Plantation + Range, data = sites1pabuild)
OVEN.pa <- predict(OVENrpart.pa, newdata = sites1patest)
"L3", "Plantation", and "Range" are all character variables, the rest are
numeric.
Thanks!
Jay
[[alternative HTML version deleted]]
___
ps the best advice is to download the source code
RcmdrPlugin.x.tar.gz from CRAN (e.g. RcmdrPlugin.IPSUR), extract
the archive, and take a look at how other people have done it.
I hope that this helps, and good luck! :-)
Jay
***
G. Jay Ker
...[snipped]
>
> The /inst directory: in here you will need to put a file "menus.txt"
>
This should have been "the /inst/etc" directory, as you noted in your
original post.
Jay
__
R-help@r-project.org mailing list
t; Thank you
> Liviu
There was a closely related discussion last April:
https://stat.ethz.ch/pipermail/r-devel/2009-April/053094.html
and IIRC this was fixed for R version 2.10.
Hope this helps,
Jay
--
*******
G. Jay Kerns, Ph.D.
Associate Professor
De
solution:
library(prob)
urnsamples(0:2, size = 3, replace = TRUE, ordered = FALSE)
Regards,
Jay
***
G. Jay Kerns, Ph.D.
Associate Professor
Department of Mathematics & Statistics
Youngstown State University
Youngstown, OH 44555-0002 USA
Office: 1035 C
ave a 'scale'
parameter at all.
Just be careful, and you should be fine. And IMHO, given that the PDF
of the density is shown it is reasonably clear as-is.
Best,
Jay
***
G. Jay Kerns, Ph.D.
Associate Professor
Department of Mathematic
, you may want to check out Chapter 5 (and also some of
4).
>and second, what could I
> do when I have some independent variables that are not only numerical but
> categorical too, i.e. mixed (categorical and numerical), can I still use a
> logistic regression?
Easy pea
7L, 4L, 2L), m7 = c(14L, 5L, 2L, 11L, 20L, 3L, 13L, 17L,
14L, 20L, 10L, 15L, 6L, 10L, 16L, 3L, 8L, 7L, 4L, 2L), m8 = c(12L,
15L, 4L, 11L, 18L, 3L, 12L, 16L, 14L, 20L, 9L, 12L, 16L, 10L,
17L, 3L, 18L, 9L, 7L, 3L)), .Names = c("m1", "m2", "m3", "m4",
"m5&
the joint pmf
xtabs(probs ~ Var1 + Var2, data = Probspace)
# sorted by probability
Probspace[order(Probspace$probs), ]
See the prob package for other utilities in this vein. You can make
the labels better with c("7 – 7.50", "7.50 – 8.50", "8.50 – 10.00"),
etc, in the firs
; --
> Stephen
Thanks for the brief response, Stephen; he has seen (most of) this before.
Unbelievable.
Jay
***
G. Jay Kerns, Ph.D.
Associate Professor
Department of Mathematics & Statistics
Youngstown State University
Youngstown, OH 44555-
"prob") # see section 8
?rolldie
?addrv
?marginal
Then turn off the computer, pick up a pencil, and solve the problems
by hand. Wait a minute, maybe that should be step 1?
Have fun,
Jay
***
G. Jay Kerns, Ph.D.
Associate Professor
D
>
> Debbie
>
Check the CRAN Task View on Probability distributions:
http://cran.r-project.org/web/views/Distributions.html
and in particular, check out package VGAM.
HTH,
Jay
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/
being able to figure it out in 2*10^68 years
doesn't count, but within a couple months is acceptable.)
3) does the answer change if there is a
remove(.Random.seed)
command right before the save.image() command?
Any thoughts would be appreciat
Dear Duncan,
On Thu, May 14, 2009 at 4:39 PM, Duncan Murdoch wrote:
> On 5/14/2009 3:36 PM, G. Jay Kerns wrote:
>> Question:
>> 1) can you tell me what my original set.seed() value was? (I wouldn't
>> be able to figure it out, but maybe someone can)
>
> The onl
nning of the .Rnw (but after the set.seed)... and did not
distribute this information to the students... that would make it
much harder, yes?
Any ideas that are even better than this?
Conceivably, some of my students will be searching these archives in
the future; please feel free to res
a in 2008 for the R-wiki:
http://tolstoy.newcastle.edu.au/R/e5/devel/08/10/0481.html
There were no responses and I ran out of time to continue working on
it myself. If you are interested in proceeding along these lines then
I have some more ideas and would be willing
>
> Jay: why not post your R-books how to on the wiki itself???
Because I thought that it would be better to write the instructions in
R-wiki language that anybody could modify rather than post a PDF by
me.
Here is what I had in mind:
http://wiki.r-project.org/rwiki/doku.
every possible variant of
questions (and answers) about this, and you haven't said _exactly_
what you are looking for. In the absence of an already posted
solution, please specify exactly what you want and I'll wager an R
Ninja could dispatch it in moments.
Regards,
Jay
***
Jason,
On Fri, May 29, 2009 at 5:58 PM, Jason Rupert wrote:
>
> Jay,
>
>
> Thanks much for the reply. I think you are right about the prob.
> Unfortunately, I was not able to find the old emails I had discussing the use
> of the more powerful setdiff that essentially i
Jason,
(moved back to R-help)
On Sat, May 30, 2009 at 3:30 PM, Jason Rupert wrote:
>
> Jay,
>
>
> I really appreciate all your help help.
>
> I posted to Nabble an R file and input CSV files more accurately
> demonstrating what I am seeing and the output I
ions ->
Continuous -> Exponential -> Sample...
If you were wanting to normalize the row sums (or investigate the
sampling distribution of some other statistic, for that matter) then
check out the Sampling Distributions... menu in RcmdrPlugin.IPSUR.
Good luck,
Jay
> Basically I need to p
; checking for X... no
> configure: error: X11 not found but required, configure aborted.
> ERROR: configuration failed for package 'rgl'
> ** Removing '/usr/local/lib/R/site-library/rgl'
>
I had the same problem a few days ago, and followed the
recommendations of t
I will have to boot into Windows later this evening to investigate,
but in the meantime you may want to post your sessionInfo(). I have
run RcmdrPlugin.IPSUR on R-2.10.1, Windows XP, as recently as
yesterday during class without any problems.
Gotta go for now.
Jay
--
**
done that
because A) it confuses students and B) John went to nonzero trouble to
make the "remove" capability functional; I think it's cool.
I will think about it some more. Maybe there is a way to reorganize
it to not confuse students and simultaneously avoid conflicts with
other
like this using R is like shooting a squirrel with a ballistic
missile. Not that there's anything wrong with that; I don't like
squirrels.
Cheers,
Jay
P.S. Just kidding.
On Tue, Oct 5, 2010 at 12:06 PM, Greg Snow wrote:
> I would do it like this:
>
> library(TeachingDemos)
>
t loading it?
I ask because the vignette is built by making a special choice for
set.seed, and the workspace that ships with the package might be
interacting in an unexpected way.
Please let me know if IPSUR is the culprit.
Regards,
Jay
__
R-help@
, from the command line:
library(Rcmdr)
with(rock, scatter3d(area, peri, shape))
I hope that this helps,
Jay
__
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/post
)
You could also use copulas, but those depend on contributed packages
(and you can read more about them on the CRAN Task View for
probability distributions).
Hope this helps,
Jay
__
G. Jay Kerns, Ph.D.
Youngstown State University
http://people.ysu.edu/~gkerns/
On Sat, Nov 6, 2010 at 8:22 PM, michael wrote:
> Jay,
>
> Yes I'm looking for unif(0,1) and your method works just fine. I
> suppose your method should work for dimensions greater than 2, am I right?
>
> Michael
>
Yes, but it gets that much more tricky to spe
you store all of
those in a data frame A, say, then you can get all possible samples
from df1 with something like
apply(A, 2, function(x) df1[x]))
The good news is that the result will also be a matrix. Each row will
be a possible sample of size 6 from df1.
Hope this helps,
Jay
P.S. Note th
aming
scheme for menus is also updated to be consistent with the other
plugins in particular, and naming conventions in general.
Cheers,
Jay
P.S. If you are thinking to print parts of IPSUR yourself then I
recommend the publisher-quality PDF linked from the Downloads secti
lly, it depends on what you mean by 'get through'.)
>
> Also, do you have other suggestions for texts, manuals, web sites, etc. that
> would introduce statistics and R simultaneously?
Have you seen this?
http://rwiki.sciviews.org/doku.php?id=books:intrstat
Good luck,
Jay
anks, Joh
>
[snip]
Have a look at tk_choose.files in the tcltk package.
library(tcltk)
?tk_choose.files
Good luck,
Jay
*******
G. Jay Kerns, Ph.D.
Associate Professor
Department of Mathematics & Statistics
Youngstown State University
Youngstown,
Dear Thibault,
On Tue, Apr 12, 2011 at 5:06 PM, Thibault Grava
wrote:
> Hello Jay,
>
> I know that an old post but I really need to calculate KMO for my data and
> this is the only thread I found on the subject.
>
> I'm a very newby to R so sorry for the odd questions.
&g
ectly, the same thing happened to me when my Plugin
removed a menu which was used (that is, added to) by some other Plugin
(in particular, RcmdrPlugin.IPSUR deleted a menu to which
RcmdrPlugin.HH added a submenu.)
The solution in my case was "don't do that".
Can't really
hi,
i'm trying to understand r data structures. i see that vectors, matrix,
factors and arrays have a "dimension."
there seems to be no mention of dimensionality anywhere for lists or
dataframes. can i consider lists and frames to be of fixed dimension 2?
s.htm
-
On Aug 16, 2012, at 11:49 AM, Schumacher, Jay S wrote:
>
>
> hi,
> i'm trying to understand r data structures. i see that vectors,
> matrix, factors and arrays have a "dimension."
> there seems to be no
quire as to the
length of a list, so might, if forced to, admit that lists have one
dimension, length. But I do not think it is helpful to think of lists as
having dimension. Certainly, lists do not have two dimensions.
--
Don MacQueen
Lawrence Livermore National Laboratory
7000 East Ave., L-627
quot;, such as TODO lists, agendas,
capturing-archiving, calendar integration, time tracking,... These
could possibly address some of the other issues you mentioned.
Good luck, and Happy New Year.Jay
*********
G. Jay Kerns, Ph.D.
Mathematics and Statistics
Youngstown Stat
rlgamma in the
actuar package are correct.
There may be other people who can identify numerical issues that slow
down the convergence rate even more.
Best,
Jay
P.S. the MGF of gamma is all you need for the moments of loggamma...
E(X^k) = (1 - k / 2.6)^(-25)
On Sat, Jun 14, 2008 at 4:54 AM, mich
you need the
columns collapsed, then you can use the
apply(B, 1, paste, sep = "", collapse = "")
command that Johannes suggested. Details are in the prob package vignette,
vignette("prob")
I hope that this helps,
Jay
* fix: As it happens, your particular q
list(D), FUN = sum)
The output will be a data frame with the unique rows, and a column at
the end labeled "x" with the frequency of each unique row.
Once you get this you can convert to a list, manipulate, etc. I am
sure that there exist faster/better methods.
Best,
Jay
On Dec 21, 200
comments are welcomed. I would appreciate hearing
about your experiences with it in the classroom and elsewhere.
The audience for this package would include teachers and students of
elementary probability, or simply anyone wanting to dabble with
p
n assumption of independence
between X1 and X2. Otherwise, you would seem to need to know the joint
distribution of X1 and X2, which isn't mentioned in your problem
description. Is there additional information about the problem that
would suggest adding p2 rather than multiplying?
I hope this help
le functions.
After installing Rcdmr and RcmdrPlugin.IPSUR (version >= 0.1-5), do
library(RcmdrPlugin.IPSUR)
and after the Commander restarts take a look at the "Distributions"
menu. You can plot any of the above. In addition, you will have
uld simply be P( X < x | Y=y ) *
f(y), where f(y) is the marginal pdf of Y (a dnorm).
Note that the above is assuming that y is a fixed constant; if not,
then you may want to check out the Ryacas package.
I hope that this helps,
Jay
***
the inverse CDF approach, one can
generate bivariate samples that have any given marginal distribution.
See Nelson for details.
Best,
Jay
>
> thanks
> A.S. Qureshi
>
> __
> R-help@r-project.org mailing list
> https://stat.ethz.ch/m
nzi.psych.upenn.edu/R/R-devel/archive/26683.html
and as a consequence of that discussion:
library(prob)
setdiff(A,B)
Best,
Jay
*******
G. Jay Kerns, Ph.D.
Associate Professor
Department of Mathematics & Statistics
Youngstown State University
Y
Dear Kimmo,
It doesn't appear that anyone has yet mentioned pareto.chart() in the
qcc package; it may serve your purposes. Please note that it does not
require the data frame to be ordered beforehand.
x <- DATA[[2]]
names(x) <- DATA[[1]]
library(qcc)
pareto.chart(x)
Best wishes,
Ja
but am I missing something? and in
particular, is there an elegant way to check in the case that the mode
of the vector is not already known?
Thanks in advance for any insight you may have.
Best,
Jay
*******
G. Jay Kerns, Ph.D.
Associate Pro
Dear Steven,
> length(x)
> Does this cover all your use
> cases?
Yes, and thanks again to everybody else who later replied. I had
falsely imagined something so much more complicated...!? Next time,
I will wait 8*runif(1) before posting. :-)
Best,
Jay
> HTH
>
&g
Dear Everybody,
Thanks, length() is the answer.
Best,
Jay
On Mon, Nov 24, 2008 at 1:15 PM, Ted Harding
<[EMAIL PROTECTED]> wrote:
> On 24-Nov-08 17:41:25, G. Jay Kerns wrote:
>> Dear R-help,
>> I first thought that the empty set (for a vector) would be NULL.
>>
>
; I need to generate a matrix that looks like this:
>> A B C
>> A B D
>> A B E
>> A C B
>> A C D
>> A C E
>> A D B
>> A D C
>> A D E
>
> Hi,
>
> Does this do what you want?
>
> expand.grid(letter
Dear Kingsford,
You are quite right, my mistake:
urnsamples(ID, size = 3, replace = FALSE, ordered = TRUE)
Thanks.
Jay
On Wed, Dec 3, 2008 at 12:04 AM, Kingsford Jones
<[EMAIL PROTECTED]> wrote:
> However, I believe Brandon was trying to get the permutations of size
> 3,
s approach to
> a variable number of arguments to expand.grid.
>
> Hope this helps,
>
> baptiste
>
Here is another way:
library(prob)
urnsamples(1:3, size = 2, ordered = FALSE, replace = TRUE)
You can convert to a matrix with as.matrix(), if desired.
Regards,
Jay
--
***
ve in mind, it may be useful to think carefully about
_why_ it is desired to have all possible combinations. In some
circumstances, it is good enough to randomly generate combinations and
draw inferences from a sampling distribution of some sort associated
wi
"prob" argument accordingly.
> And what # of
> combinations should I use as the cutoff for doing a sampling instead
> of testing every possible combination?
There isn't a definitive answer to that question. Assuming that it is
desired to find an 'average
G*Power 3 is free software for Mac and PC, see
http://www.psycho.uni-duesseldorf.de/abteilungen/aap/gpower3/
Jay
On 9/16/07, MATTHEW BRIDGMAN <[EMAIL PROTECTED]> wrote:
> Is there a way to calculate power for repeated
> measures ANOVA (2 groups x 7 observations)? I have
> searche
Hi Rainer,
The distr package can calculate the distribution for you:
library(distr)
X <- Binom(size = 7, prob = 0.3)
Y <- Binom(size = 11, prob = 0.5)
Z <- X + Y
d(Z)( 0:18 ) # the pmf
r(Z)( n = 5 ) # random variates
Please note, however, that size and prob must be of length 1.
Best
I am working with a RCB 2x2x3 ANCOVA, and I have noticed a difference in the
calculation of sum of squares in a Type III calculation.
Anova output is a follows:
> Anova(aov(MSOIL~Forest+Burn*Thin*Moisture+ROCK,data=env3l),type=3)
Anova Table (Type III tests)
Response: MSOIL
el,"cB vs. A&C")
linearHypothesis(model,"cA vs. C")
#Sum of squares of the individual contrasts do not add to the main
effect of c
John J. Wiley, Jr.
PhD Candidate
State University of New York
College of Environmental Science and Forestry
Department of Environmenta
variances are similar. Is there
another metric I am missing for a continuous covariate?
Cheers,
John
From: John Fox [j...@mcmaster.ca]
Sent: Tuesday, October 09, 2012 2:59 PM
To: John Jay Wiley Jr.
Cc: r-help@r-project.org
Subject: RE: [R] car::linearHypo
: Tuesday, October 09, 2012 4:00 PM
To: John Jay Wiley Jr.
Cc: r-help@r-project.org
Subject: RE: [R] car::linearHypothesis Sum of Sqaures Error?
Dear John,
> -Original Message-
> From: John Jay Wiley Jr. [mailto:jwile...@syr.edu]
> Sent: Tuesday, October 09, 2012 3:37 PM
> To: J
like to know if there is any solution to fix it without giving the read
access for these two files /etc/passwd and /etc/group in the apparmor profile
as I did with ubuntu 18.04 and R 3.6.3. Thanks!
Best Regards!
Jay Gu
[[alternative HTML version deleted]]
__
Hi Ivan,
I'm running the R within docker container. Do you have any idea about it?
Thanks!
Best Regards!
Jay Gu
-Original Message-
From: Ivan Krylov
Sent: Wednesday, August 9, 2023 3:15 AM
To: Gu, Jay via R-help
Cc: Gu, Jay
Subject: Re: [R] R library highcharter function high
101 - 167 of 167 matches
Mail list logo
