similar to: model selection vs. H0 based testing

Hello, I want to do simulations with an external simulation model. For this purpose I wrote a R-script which makes the pre- and post-processing of the data (generate scenarios, analyze results). The R-script writes the input file for the simulation with e.g. write(x, file="test.dat", ncolumns=10) The problem is, however, that this command (and even constructions with cat) will add a

To use the modify the solution from Tony and I so that you can pass the name of the function, rather than the function itself, like this: x <- 7 fx <- function(y) print(x*y) f <- function(fun, x) { fun <- get(fun) environment(fun) <- environment() do.call("fun",list(3)) } f("fx",2) --- Date: Thu, 11 Mar 2004 08:22:45 +0100 From:

Hi all, we found, that R is not only *the* superior data analysis, graphics and statistics package, but it is also a general tool for implementing, running and teaching ecological models. Last week we held our second course "Modelling in Limnology" using R. Some material (code and docs) was produced for this purpose and may be interesting to someone else, so I put the tutorial

== Summary == * Problem: BibTeX entries extracted from R packages via citation() require too much manual editing to be of general use. * Proposal: Date: fields should be made mandatory in package DESCRIPTION files, perhaps beginning with warnings from R CMD check * Proposal: Package authors should be encouraged to use a (new) Contributors: field in the DESCRIPTION file rather than packing all

I have a problem, that the plot of mathematical formula sometimes crashes R on Linux (SuSE 6.2 and 7.0 on a Dual PIII/500 machine with 1GB RAM). The problem occurs for example in demo("graphics"), i.e. in the last part were the mathematical symbols are plotted. There are several circumstances to reproduce this segmentation fault, but the following is one of the most simple examples. The

I have carried out a spatial interpolation (Kriging) using surf.gls() The contour plot looks very plausible and the image plot is o.k too. The problem is, hovever, that I need the vectorized isoplethes, as in contour(). Unfortunately, I found no solution to get the internal data from contour(). Is there a possibility to extract contour lines as x-y-polygon from a grid without hacking plot3d.c?

Dear useRs, There is a new version 0.1-2 of the package surv2sample available on CRAN. Users of the previous versions should update because a bug in the function cif2.ks has been fixed. General information about the package: surv2sample provides various two-sample tests for right-censored survival data. Three main areas and corresponding methods are: * comparison of two survival

Dear useRs, There is a new version 0.1-2 of the package surv2sample available on CRAN. Users of the previous versions should update because a bug in the function cif2.ks has been fixed. General information about the package: surv2sample provides various two-sample tests for right-censored survival data. Three main areas and corresponding methods are: * comparison of two survival

Hello, I was wondering if anyone could please help me with some simple questions regarding ANCOVA and the assumption of homogeneity of slopes. The standard design of ANCOVA assumes the homogeneity of regression coefficients of the different groups. This assumption can be tested using the factor ?? covariate interaction, which should subsequently be removed. However if this assumption is not

Dear all, I've been delighted to just notice that Cohen's formulas for Effect Size 'w' and the associated power have been implemented in the 'pwr' package (thanks to St?phane Champely and others).. There is one aspect, though, that perplexes me. I'm doing some last minute post hoc analyses, meaning that my sample size (N=3404) has been long fixed, and I'm

Hello, I try to synthesize harmonic functions with a subset of frequencies determined by fast fourier transform. I wrote two different functions, a looped version and a matrix multiplication version. As an example the looped version takes 5 sec and the matrix-version takes 3-4 sec (R 1.4, Athlon 1.2 GHz, 256 MB, Win ME), but the latter needs huge amount of memory due to the matrices. So, none of

I've just found the lavaan package, and I really appreciate it, as it seems to succeed with models that were failing in sem::sem. I need some clarification, however, in the output, and I was hoping the list could help me. I'll go with the standard example from the help documentation, as my problem is much larger but no more complicated than that. My question is, why is there one latent

Hi there, I want to generate a random point pattern using the Neyman-Scott cluster process using spatstat package in R. After running the following procedures, why i can not see any figures? > nclust <- function(x0, y0, radius, n) {return(runifdisc(n, radius, x0, y0))} > nclust function(x0, y0, radius, n) {return(runifdisc(n, radius, x0, y0))} > X <- rNeymanScott(10, 0.2, nclust,

?tryCatch -- Sent from my phone. Please excuse my brevity. On October 20, 2017 7:37:12 AM PDT, Evangelina Viotto <evangelinaviotto at gmail.com> wrote: >Hello I?m need fitt growth curve with data length-age. I want to >evaluate >which is the function that best predicts my data, to do so I compare >the >Akaikes of different models. I'm now need to evaluate if changing the

Dear R developers, I've got an information from Prof. Ripley regarding a bug found with AdressSanitizer in one of our packages. It is now fixed, thank you for this information. Now, I would like to run AddressSanitizer myself before submitting the patched package to CRAN. Is there a recommendation of a suitable Linux distribution with gcc 4.8, ideally an ISO image or (even better) a

Dear R developers, I recently observed a NOTE on several MaxOS X package checks: sh: line 1: gs: command not found !!! Error: Closing Ghostscript (exit status: 127)! /usr/bin/texi2dvi: thumbpdf exited with bad status, quitting. See for details: http://www.r-project.org/nosvn/R.check/r-release-macosx-ix86/simecol-00check.html or

Hello I?m need fitt growth curve with data length-age. I want to evaluate which is the function that best predicts my data, to do so I compare the Akaikes of different models. I'm now need to evaluate if changing the initial values changes the parameters and which do not allow to estimate the model. To do this I use the function nls(); and I randomize the initial values (real positive number).

Dear R developers, we have currently a 'mysterious' test problem with one package that successfully passed the tests on all platforms, with the only exception of Solaris x86 where obviously one of our help examples breaks the CRAN test. As we don't own such a machine I want to ask about a possibility to run a few tests on such a system: r-patched-solaris-x86 An even more recent

Dear R developers, I observed that the html help page index entry "Read overview or browse directory" for package vignettes is missing in recent R-devel. This happened on two independent computers (WinXP Prof. SP3, German) with R-devel compiled from sources svn rev. 48125 resp. 48128 It's the same for my own and also for more prominent packages as well (e.g. grid). The

I am using lm simulation in R and try to find the AICc and Akaike weight of the model. I searched out that using package "AICcmodavg" AICc is easily to get. I wonder how can I get the "Akaike weight", any function I can use to generate it? Thanks in advance. Sunny [[alternative HTML version deleted]]