1. Martin Maechler's comments should be taken as replacements
for anything I wrote where appropriate. Any apparent conflict is a
result of his superior knowledge.
2. 'eigen' returns the eigenvalue decomposition assuming the
matrix is symmetric, ignoring anything in m[upper.tri(
On 2/23/2010 5:07 AM, Cedrick Johnson wrote:
library(rant)
I'm sorry: I'm not familiar with the "rant" package. I couldn't find
it on CRAN.
I recall just over a year ago (Fall08) when I was new to the R
language. I realize there's a steep curve in any language, especially
with R and the m
I have a concern with the sentiment expressed: When one is new
to a language, it is often not easy to know where to start in the
available documentation. It's hard to RTFM when one does not know which
FM2R.
I currently subscribe to email help for three other open source
collabo
Hi, Ivan:
p.s. If you'd like to generate reports using LaTeX, I suggest you
also try "(LaTeX <- ???LaTeX)".
###
Are you trying to read standard text or csv files? If yes, then try
"?read.table" at a command prompt.
Or are you trying to connect directly from R
Hi, Ivan:
Are you trying to read standard text or csv files? If yes, then try
"?read.table" at a command prompt.
Or are you trying to connect directly from R to a database system?
If yes, which database system? This can be done from R, but as far as I
know, finding the functions an
raves
Prof Brian Ripley wrote:
On Tue, 9 Feb 2010, spencerg wrote:
Can one write a zip file from R?
I want to create a file with a name like "dat.zip", being a zip
file containing "dat.csv". I can create "dat.csv", then call
"system('zip -r9 d
Can one write a zip file from R?
I want to create a file with a name like "dat.zip", being a zip
file containing "dat.csv". I can create "dat.csv", then call
"system('zip -r9 dat.zip dat.csv')". Is there a better way?
I can use "gzfile" to write a gz file, but I don't kno
Have you considered the "fda" package and the companion book,
"Functional Data Analysis with R and Matlab" (Springer, 2009) by Ramsay,
Hooker and Graves? This will NOT help you directly with bathymetry =
f(long, lat), but will help with b=f(x)+e AND with translations between
Matlab and R.
standard square root computation requires an iteration. Spencer
Peter Dalgaard wrote:
Mads Jeppe Tarp-Johansen wrote:
To R-helpers,
R offers the polyroot function for solving mentioned equations iteratively.
However, Dr Math and Mathworld (and other places) show in detail how to
solve men
To me, R is the language of choice for a rapidly increasing number
of people involved in new statistical algorithm development. If they
are happy with the tools they currently use, learning R may be a lot of
pain for little gain.
However, if they want to stay current with the late
Hi, Chuck:
Thanks very much, but why do I get "package 'expm' is not
available" from
install.packages("expm",repos="http://R-Forge.R-project.org";)?
Best Wishes,
Spencer Graves
Charles C. Berry wrote:
On Sun, 1 Nov 2009, spencerg
A question, a comment, and an alternative answer to matrix^(-1/2):
QUESTION:
What's the status of the "expm" package, mentioned in the email you
cited from Martin Maechler, dated Apr 5 19:52:09 CEST 2008? I tried both
install.packages('expm') and
install.packages("expm",repos="http://R-For
'isi' has only 14 values, 'age' has 31. The help page for
'smooth.basis' says, "argvals: a vector of argument values correspond to
the observations in array 'y'."
I think 'isi' and 'age' should have the same length.
I apologfize if the error message was, uh, cryptic.
Hi, Ben:
Which chapter in which of Ramsay's books?
For his most recent book ("Functional Data Analysis with R and
Matlab", with Giles Hooker and your's truly), there is one script file
for each chapter with names like "fdarm-ch01.R", ... "fdarm-ch11.R".
These script files "rep
Charles C. Berry wrote:
On Thu, 22 Oct 2009, Ben Bolker wrote:
Allan.Y wrote:
Hi everyone,
I am wondering if there exists a stepwise regression function for the
Bayesian regression model. I tried googling, but I couldn't find
anything. I know "step" function exists for regular stepwise
In addition to the taskview (below), you might consider using the "sos"
package, something like the following:
sa <- ???'spatial autocorrelation' # 58 matches
sc <- ???'spatial correlation' # 181 matches
s. <- sa|sc # union of the to 'findFn' objects
summary(s.) # 224 links in 57 packages
ins
There is a huge literature on time series analysis with many
functions and contributed packages in R, including many different ways
"to estimate a recursive model".
If you are working with "asset returns", I suggest you start with
Diethelm Würtz, Yohan Chalabi, William Chen, Andrew
library(sos)
tm <- findFn('text mining')
tm
This produced 15 matches, which you could also find using
"RSiteSearch('text mining', 'function')". The difference is that
findFn{sos} displays the results in a table sorted to place the package
with the most matches first. In this case, ther
library(sos)
a <- ???accelerometry # no matches
d3 <- ???'3D' # 411 matches; retrieved 20 pages 400 matches
d3 <- ???'3D'(99) # retrieved all 411 matches
summary(d3) # first 13 of 144 packages
findFn2xls(d3) # write d3.xls containing sheets PackageSum2, findFn, call
installPackages(d3) # install
RSiteSearch('exponential of a matrix')
produced 982 matches.
RSiteSearch('{exponential of a matrix}')
produced 13.
David Winsemius wrote:
On Sep 30, 2009, at 9:18 PM, Kon Knafelman wrote:
Hi Guys,
Im trying to find the exponential of a matrix.
Can someone please help me
Sylvester's formula (http://en.wikipedia.org/wiki/Sylvester%27s_formula)
applies to a square matrix A = S L solve(S), where L = a diagonal matrix and S =
matrix of eigenvectors. Let "f" be an analytic function [for which f(A) is well
defined]. Then f(A) = S f(L) solve(S).
We can
The best reference I know for this is something I wrote with Jim
Ramsay and Giles Hooker: Functional Data Analysis with R and Matlab
(Springer, 2009). Others may have better material.
After "install.packages('fda')", I suggest you try
"system.file('scripts', package='fda')", as s
RSiteSearch('java', 'fun') produced 144 matches for me just now.
The "sos" package can help organize study of those 144 matches:
library(sos)
j <- ???java # RSiteSearch('java', 'fun') -> a data.frame of class '
findFn'
summary(j) # 40 packages
installPackages(j) # install packages with
Searching help pages of contributed packages just got
easier with the release of the new "sos" package. This is a
replacement for and substantial enhancement of the existing
"RSiteSearch" package. To learn more about it, try
vignette("sos").
We hope you find this as useful as
Hi, Giovanni:
1. I've used primarily "sspir" and "dlm". There is apparently a
new release of "sspir", which I have not used. When I last tried
"sspir", it did NOT have a forecast function, while "dlm" did. The two
packages have functions with the same name but incompatible code. It
This may be overkill for your application, but you might be
interested in the "fda" package, for which a new book appeared a couple
of months ago: "Functional Data Analysis with R and Matlab" (Springer
Use R! series, by Ramsay, Hooker and Graves; I'm the third author).
The package inclu
ate a
functional time series. I may also be able to get traffic counts, so I
could model deflection as a function of traffic.
Thanks again to Middleton, Bowlingson, and Austenfeld.
Best Wishes,
Spencer
Barry Rowlingson wrote:
On Sun, Sep 6, 2009 at 1:07 AM, spencerg wrote:
What software exists for digitizing video to quantify the motion
of specific features in the image? I might be willing to use something
that's NOT in R, though I'd prefer something in R (or at least with an R
intereface).
Thanks,
Spencer Graves
__
Bryan Hanson wrote:
Looks like the discussion is no longer about R Style, but S3 vs S4?
To that end, I asked more or less the same question a few weeks ago, arising
from the much the same motivations. The discussion was helpful, here's the
link:
http://www.nabble.com/Need-Advice%3A-Consider
S3 is very easy to change; S4 is very difficult. This provides
advantages and disadvantages for both. Some people swear by one and
curse the other -- from both sides.
S4 is newer, and I have had problems in the past finding out what
S4 methods are available and finding acceptabl
You might also try
library(RSiteSearch)
de <- RSiteSearch.function("differential equation")
# This is different from
des <- RSiteSearch.function("differential equations")
# Combine as
de. <- de | des
# display:
HTML(de.)
Hope this helps.
Spencer Graves
Bert Gunter wrote:
For n
Dear Mike:
I don't know.
1. What specific error message do you get?
2. Your example is too long for me to parse, especially with
the color being stripped before I saw it.
3. Have you tried using "debug(cro.etest.grab)", then
walking through t
Hi, Terry:
Thanks for the comments. I too vastly prefer S3 to S4. Your
comparison is based on much greater experience than mine.
Could you please check the link you sent? I couldn't get it to
work.
Thanks again.
Spencer
Terry Therneau wrote:
For 90 percent of
Have you looked at RExcel and the RExcelInstaller package?
There is now a companion book: Heiberger and Neuwirth (2009) R
Through Excel: A Spreadsheet Interface for Statistics, Data Analysis,
and Graphics (Springer)? Both Amazon and the Springer web site say it's
not yet availabl
Have you looked at Heiberger and Neuwirth (2009) R Through Excel:
A Spreadsheet Interface for Statistics, Data Analysis, and Graphics
(Springer)?
Both Amazon and the Springer web site say it's not yet available.
However, the Springer booth at the Joint Statistical Meetings was
se
There may also be a difference in reliability, which would not so
easily be measured by an individual user. I've selected the closest
geographically until it seemed to be down, then tried the second
closest, etc. This could be automated centrally, but then you'd have to
deal with the hum
Have you considered the "timeDate" package?
Spencer
Denis Chabot wrote:
Thank you very much Duncan.
I'll follow your suggestion.
Why do I want to do what the designer did not think anyone would want
to do? I have data acquisition equipment taking measurements every 15
min or so for d
You also need to know which days are trading vs. non-trading
days. To find functions to provide this additional information, I first
tried the following:
library(RSiteSearch)
td <- RSiteSearch.function('trading days')
HTML(td)
This produced nothing. Therefore, I divided the task
I can't identify the problem. The package author and maintainer,
Paul Gilbert, might be able to help.
Have you tried "debug(SS)"? This will allow you to walk through
the function line by line looking at things, etc. This often produces
enlightenment.
Alternatively, have
Hello:
1. I tried the example. The first thing I noticed was that
"library(fSeries)" loaded with the following comment:
# The new version of 'fSeries' has been renamed to 'timeSeries'
2. With this, I tried "library(timeSeries)". Then your first
example worked without error
1. What does "i" in your formula represent? Have you worked the
examples in the "nls" help page, and do you understand how it works?
"nls" tries to do vector computations.
2. Unfortunately, "nls" often quits with errors like "singular
convergence". A standard way around that p
There are many more things available than just 'nls'. An overview
is available at "www.r-project.org" -> CRAN (select your favorite
mirror) -> "Task Views" (on the left, third line under CRAN) ->
"Optimization Optimization and Mathematical Programming".
For a possibly more di
I have not used 'sna'.
Have you tried using "debug" to walk through the code line by
line, examining and even changing things at will?
For example, how big is "rho", passed as starting values to
"optim"? If that matches the size of your adjacency matrix, it could
expose a
Have you plotted the data? There are two standard sources of
non-convergence problems like this: First, there may not be enough
information in your data to estimate all four parameters. Second, if
that's not the case, then your starting values are not appropriate.
I routinely u
Your script ran for me without error under the following configuration:
> sessionInfo()
R version 2.9.1 (2009-06-26)
i386-pc-mingw32
locale:
LC_COLLATE=English_United States.1252;LC_CTYPE=English_United
States.1252;LC_MONETARY=English_United
States.1252;LC_NUMERIC=C;LC_TIME=English_United Sta
Hello:
For any package on CRAN, you can download the "*.tar.gz" source
file. When this is unzipped, the result will typically include a
subdirectory "src", containing the source code for routines written in
Fortran or C or C++.
You didn't say where one could find the function
Dear Andriy:
1. The help file for "fit.NH" says the first argument should be
a "vector of data". Consider the following modification of the example
on that help page:
> rs <- matrix(rseries, 10, 202, dimnames=list(letters[1:10], 1:202))
> rs. <- fit.NH(rs)
Error: cannot allocate
Your example is not self-contained, which makes it more difficult
to answer your question.
However, I believe you should be able to answer the question by
studying the examples on the help pages for "gnls" and possibly
"corCompSymm", using the following additional commands:
# T
Hello, Yuhan:
If I wanted to get something sensible today, I'd do ordinary least
squares using lm(y~x), the garchFit on the residuals. This will give
you a reasonable answer except that the confidence intervals from "lm"
will not be accurate. I'd want to do normal probability plots of
I have not seen a reply to this post, so I will offer a couple of
comments:
To look for other garch capabilities in R, I tried the following:
library(RSiteSearch)
gch <- RSiteSearch.function('garch')
summary(gch)
HTML(gch)
This found 111 help pages in 18 packages containi
I have not seen a reply to this question, so I will offer a
comment; someone who knows more than I may correct or add to my comments.
There are many different kinds of splines. Perhaps the most
common are B-splines, which sum to 1 inside their range of definition
and are 0 outsid
I just tried the following:
library(RSiteSearch)
sound <- RSiteSearch.function('sound')
HTML(sound)
This identified 117 help pages in 40 packages containing the word
"sound". This included 30 help pages in a package called "sound", which
is "A Sound Interface for R". It also includ
Hi, Mario:
The following suggest to me that these data may not be available
in any contributed R package:
library(RSiteSearch)
fbi <- RSiteSearch.function('FBI Homicide')
hits(fbi) # 0
hom <- RSiteSearch.function('Homicide')
hits(hom) # 4
HTML(hom)
The 4 items exposed via the l
I agree that the documentation is a primary source, made now more
accessible with the availability of the RSiteSearch package that allows
more structures searches of the help pages in contributed packages than
previously available. On the other hand, it can sometimes be difficult
to get t
The documentation for the "nlme" package was improved a few years
ago by an informal process of this nature. When I first started
following r-help, I answered many questions suggesting in part the the
person read some portion of Pinheiro and Bates (2000). At that time,
none of the help p
Hi, Justin:
In addition to the example below, Thomas Petzoldt also has two
vignettes included with his "simecol" package. These might help you
more easily use his package.
Have you looked for what you want using the "RSiteSearch.function"
in the "RSiteSearch" package? Conside
The "1L", etc., forces the number to be integer:
> is.integer(1)
[1] FALSE
> is.integer(1L)
[1] TRUE
Hope this helps.
Spencer
Paul Smith wrote:
On Thu, Jun 11, 2009 at 9:56 PM, Brecknock, Peter wrote:
Apologies if this is not the correct list for this question.
The Rglpk pac
Dear In-Sun:
I have not seen a reply, so I will offer some suggestions, hoping
I can help without understanding all the details.
1. Have you run that code with "options(warn=2)"? It
produced over 50 warnings for me, and "options(warn=2)" will convert the
warnings to erro
Sundar Dorai-Raj taught me to do the following:
Rdir <- "c:\appropriatepath\R"
Rfiles <- dir(Rdir, pattern='\\.R$', full.names=TRUE)
invisible(lapply(aTR, source))
The "invisible" suppresses the garbage while still displaying
error messages.
H
It looks to me like you need to add "x" to the arguments in your
call to "constrOptim", something like the following:
constrOptim(beta_i, myFunction, NULL, ui, ci, mu = 1e-04, control = list(),
method = "Nelder-Mead",outer.iterations = 100, outer.eps = 1e-05, x=x)
The first "x"
"loc.gov" lists two books on this topic:
* Mosteller and Wallace (1964) Inference and disputed
authorship, The Federalist (Addison-Wesley)
* Mosteller and Wallace (1984) Applied Bayesian and
classical inference : the case of the Federalist papers, 2nd ed. of
Inf
You example is NOT self contained, which means that any potential
respondent must guess what you mean by "a function with a variable of
almost 200,000". The following clarifies this:
> start0 <- rep(1, 20)
> msLE2 <- function(x)sum(x^2)
> nlminb(start=start0, msLE2, control = list(x
Dear Erin:
I do not have a direct answer to your question, but the following
quickly identified for me just now 244 help pages in contributed
packages that mentioned either "path analysis" or "SEM", sorted to put
first the package with the most such hits:
pa <- RSiteSearch.function(
I'm not sure what you are asking, especially since I do not have
access to "regaccdis". However, will something like the following do
what you want?
caeLvls <- c(1, 5, 10)
for(i in 1:3)
coef[i,2:4] <-
coef(plm(ff,data=regaccdis,na.action=na.omit,model="pooling",subset=(regaccdis$caedois==cae
s,
Spencer
Gabor Grothendieck wrote:
Your google search or
daily data site:https://stat.ethz.ch/pipermail/r-sig-finance/
both seem to give reasonable results.
On Mon, May 25, 2009 at 9:45 AM, spencerg wrote:
Hi, Gabor, et al.:
What's the best way to search the R-SIG-Finance
Frank E Harrell Jr wrote:
spencerg wrote:
Dear Frank, et al.:
Frank E Harrell Jr wrote:
Yes; I do see a normal distribution about once every 10 years.
To what do you attribute the nonnormality you see in most cases?
(1) Unmodeled components of variance that can generate
Hi, Gabor, et al.:
What's the best way to search the R-SIG-Finance archives? Are
R-SIG-* included in any of the main R search engines (at
"http://www.r-project.org/search.html";)? Gmane has an option for
"gmane.comp.lang.r.finance", but when I tried it just now, I got, "No
such group
Dear Frank, et al.:
Frank E Harrell Jr wrote:
Yes; I do see a normal distribution about once every 10 years.
To what do you attribute the nonnormality you see in most cases?
(1) Unmodeled components of variance that can generate errors
in interpretation if ignored, even
Hello:
Have you looked at Pfaff (2008) Analysis of Integrated and
Cointegrated Time Series with R, 2nd ed. (Springer)?
I have not read this book, but the title and table of contents
sounds like it contains many alternative answers to your question, with
the "best" among those b
earCen"), row.names
= c(4L,
18L, 26L, 40L, 67L, 75L, 94L, 104L, 118L, 126L, 140L, 167L, 175L,
194L, 204L, 218L, 226L, 240L, 267L, 275L, 294L, 304L, 318L, 326L,
340L, 367L, 375L, 394L, 404L, 418L, 426L, 440L, 467L, 475L, 494L,
504L, 518L, 526L, 540L, 567L, 575L, 594L, 604L, 618L, 626L, 640L,
The first exaample on the "lmer" help page uses a formula
"Reaction ~ Days + (Days|Subject)". Here, "Subject" is the name of a
column in the data.frame "sleepstudy", with levels "308", "309", ... .
Does this answer your question? If no, please provide commented,
minimal, self-con
Have you tried the "maxLik" package? If that is not adequate, you
can check "CRAN Task View: Optimization and Mathematical Programming"
(http://cran.fhcrc.org/web/views/Optimization.html). Or try the new
"RSiteSearch.function" in the "RSiteSearch" package. If none of these
adequate, ple
If you have the "RSiteSearch" package installed, you can do the
following:
library(RSiteSearch)
nrow(nll <- RSiteSearch.function("nonlinear regression with latent"))
HTML(nll)
This just produced 8 hits for me. If this doesn't solve your
problem, you might try other search terms.
The following searches for help pages in contributed packages
including terms "stationarity" or "unit root":
library(RSiteSearch)
st <- RSiteSearch.function('stationarity')
ur <- RSiteSearch.function('unit root')
ur. <- st|ur
nrow(st) # 68
nrow(ur) # 122
nrow(ur.)# 180
HTML(st&ur) # Show
Dear Dieter:
Thanks for the correction. I failed to test the code as written
before I posted it.
Spencer Graves
Dieter Menne wrote:
spencerg prodsyse.com> writes:
de <- RSiteSearch.function("differential equation")
des <- RSiteSearch.function(&quo
There is also the relatively new "RSiteSearch" package. It's
"RSiteSearch.function" searches only help pages of contributed packages
but returns the result in a data.frame (of class "RSiteSearch") sorted
to put the most interesting package first with help pages sorted within
packages. If
My favorite tool for finding things like this is
"RSiteSearch.function" in the "RSiteSearch" package. For the objects
you mention, I get the following:
library(RSiteSearch)
hits(a.s <- RSiteSearch.function("auto.stats")) # 0
hits(sx <- RSiteSearch.function("saving.x")) # 0
hits(rn <- R
There are 17 different help pages in 5 different packages citing
"Agresti and Coull". This is quickly displayed using the "RSiteSearch"
package as follows:
library(RSiteSearch)
HTML(RSiteSearch.function("Agresti and Coull"))
I have not checked all these 17, but they doubtless hel
I do not understand the term "mexval statistics".
I think you want to look for "anova.glm", fitting several models
leaving each term out one at a time in succession and then using
"anova.glm" to compare your general model with each submodel in
succession. If that does NOT give you
"Search" on the left under "www.r-project.org" provide 5 different
search engines devoted to R. The R-help archives for 4/20/2009 (+/-1
depending on your local time zone and mine) contains an interesting
discussion of this issue, Subj: "Re: [R] Two or more dimensional root
(Zero) finding
In addition to seeing the code by typing the name of the function
(and copying it from there into a file), you can also enter
"debug(fitdistr)", for example. Then the next time you use "fitdistr",
either directly or indirectly, it puts you in the environment of that
function, and you can
1. Have you worked through the examples in the maxLik help page? Your
example is sufficiently complicated that I hesitate to try it myself,
especially since I see characters in your email that are not simple
ASCII. If you can get the examples in the maxLik help page to work,
identify the diffe
You might get some information about what is available in R using
the new "RSiteSearch" package. For example:
library(RSiteSearch)
ipm <- RSiteSearch.function('interior point method')
hits(ipm) # 39
SQP <- RSiteSearch.function('SQP')
hits(SQP) # 2
sqp <- RSiteSearch.function('sequential
I suggest you try to translate your constraints into an
unconstrained constrained problem using logarithms, then do "nonlinear
mixed effects modeling" as described in chapters 6-8 of Pinheiro and
Bates (2000). To do this, I would first start with the simpler linear
estimation problem to g
Dear Doug, et al.:
What would you recommend for analyzing a longitudinal abundance
survey of 22 species, when the species were not selected at random? A
prominent scientist tried to tell me that mixed-effects modeling is
inappropriate in that case because the species were selected
pur
Dear Avraham:
For problems with many parameters to estimate, I highly recommend
Pinheiro and Bates (2000) Mixed-Effects Models in S and S-Plus
(Springer). This book includes numerous examples showing how to use the
"nlme" package. The value of this book is greatly enhanced by the
estimated in the same way but where the effect of past
betas in the transition taken into account. But how do I estimate that
matrix, is that done with a MLE,SUR or some other statistical teqnique.
Im also assuming in this example that a[t] are time invariant, which gives
W[a] = 0
Appr
Have you looked at the vignette in the deSolve package?
(deS <- vignette('compiledCode')) # opens a "pdf" file
Stangle(deS$file) # writes an R script file to "getwd()"
In spite of the name, this vignette includes an example entirely
in R. By comparing it with
Have you considered the following:
solve(qr(A), B)
I have not tried this with a small toy example, and the "qr"
documentation in the Matrix package seems to suggest it. This solves
the optimization problem you mentioned, as noted in
"http://en.wikipedia.org/wiki/Linea
Another way of making functions visible to examples is to list the
package name as "Depends" in the package "DESCRIPTION" file.
Alternatively, if you precede use in "\examples" with
"library(packageName)", you should also list it as "Suggests" in
"DESCRIPTION".
Hope this helps.
Sp
To see what's available in other packages, try the following:
library(RSiteSearch)
AFT <- RSiteSearch.function('AFT model')
summary(AFT) # 24 help files found in 8 different packages
HTML(AFT) # opens a table with 24 rows in a web browser.
There may be nothing here that will help you,
I don't understand. The help page for "maxLik" says that if the
"grad" argument is 'NULL', "numeric gradient will be used."
If this does not answer your question, please read the posting
guide "www.R-project.org/posting-guide.html" and provide commented,
minimal, self-contained, r
Hi, Stephen:
Have you discussed this with any of your professors? With a
little luck, you might find the right prof to work with who could help
you select an ecology journal and write an article for that journal
giving an overview of your package. I suggest you think in terms of a
2-
Have you worked through "vignette('dlm')"? Vignettes are nice
because they provide an Adobe Acrobat Portable Document Format (pdf)
file with a companion R script file, which you can get as follows:
(dlm. <- vignette('dlm'))
Stangle(dlm.$file)
The first of these two lines opens
-1.304531849 0.062833294 -0.143141379
[6] -0.005995477
$d
[1] 1e-04
$method
[1] "Richardson"
$method.args
$method.args$eps
[1] 1e-04
$method.args$d
[1] 1e-04
$method.args$zero.tol
[1] 1.781029e-05
$method.args$r
[1] 4
$method.args$v
[1] 2
attr(,"class")
[1] "Darray"
To look up "observer agreement", you might consider the
"RSiteSearch" package. The "RSiteSearch.function" looks only for
matches in help pages of contributed packages.
library(RSiteSearch)
oa <- RSiteSearch.function("observer agreement")
attr(oa, "hits") # 4 functions matching this term
Have you considered genD{numDeriv}?
If this does not answer your question, I suggest you try the
"RSiteSearch" package. The following will open a list of options in a
web browser, sorted by package most often found with your search term:
library(RSiteSearch)
pd <- RSiteSearch.fun
Dear Berwin: Thanks for the elegant correction. Spencer
Berwin A Turlach wrote:
G'day all,
On Sat, 09 May 2009 08:01:40 -0700
spencerg wrote:
The harmonic mean is exp(mean(logs)). Therefore, log(harmonic
mean) = mean(logs).
Does this make sense?
I think yo
not a mathematician, so
any recommendations are welcome! Thanks! -Jamie
spencerg wrote:
Are all your numbers positive? If yes, have you considered using
logarithms?
I would guess it is quite rare for people to compute likelihoods.
Instead I think most people use log(likelihoods).
Are all your numbers positive? If yes, have you considered using
logarithms?
I would guess it is quite rare for people to compute likelihoods.
Instead I think most people use log(likelihoods). Most of the
probability functions in R have an option of returning the logarithms.
1 - 100 of 111 matches
Mail list logo
