similar to: loglin( tab, margin, start = bad.start ) kills R (PR#7123)

Kai, Apologies for the double message, it didn't go to the list last time. On Thu, Oct 22, 2015 at 9:59 AM, Kai Nitschke < kai.nitschke at uniklinik-freiburg.de> wrote: > > ?.Call? calls the C/C++ function ?C_LogLin?. But when I am running it line > by line I get > the following error on line 23/24: > Error: object 'C_LogLin' not found > > Hence, my

Hi everyone, I have a question regarding a C function of the "stats" package in R. I tried to understand the ?loglin? basic function of the ?stats? package implemented in R. The implemented function itself runs without any problem (perhaps see sample). When I ran it line by line it stopped at the lines 23-24 of the loglin-function; (the following line): z <- .Call(C_LogLin,

Hi. I am trying to modify kmeans function. It seems that is failing something obvious with the workspace. I am a newbie and here is my code: myk = function (x, centers, iter.max = 10, nstart = 1, algorithm = c("Hartigan-Wong", + "Lloyd", "Forgy", "MacQueen")) + { + do_one <- function(nmeth) { + Z <- switch(nmeth, { + Z

In running the following example of a loglinear model for the Titanic data, I was surprised to see NaN reported for the Pearson chisq > loglin(Titanic, margin=list(1:3, 4)) 2 iterations: deviation 2.273737e-13 $lrt [1] 671.9622 $pearson [1] NaN $df [1] 15 $margin $margin[[1]] [1] "Class" "Sex" "Age" $margin[[2]] [1] "Survived" Tracing it back,

Summary: loglin will fail if there is a global variable called c. I have fixed this for 0.64.2 and 0.65. On Sat, 22 May 1999, Cor en Aylin wrote: > Dear Brian, > > Using the mass library for R (the most recent port dated 16 may 1999) I Could you give version numbers in such reports please: at the moment there is an official version (5.3pl037) of MASS for R, plus a pre-release of MASS

Dear R gurus, I am looking for a way to fit a predictive model for a contingency table which has counts. I found that glm( family=poisson) is very good for figuring out which of several alternative models I should select. But once I select a model it is hard to present and interpret it, especially when it has interactions, because everything is done "relative to reference cell". This

Greetings! I want to have the coefficients that R uses in shapiro.test() for the Shapiro-Wilk test for a prticular sample size, i.e. the a[i] in W = Sum(a[i]*x[i])/(Sum(x[i] - mean(x))^2) (where the x[i] are sorted). Two questions: Q1: Is there a readymade R function from which I can extract these? Q2: I was wondering if I might be able to modify the code for the function shapiro.test() so

# Why does expressing one function require(ctest) t.test # return only function (x, ...) UseMethod("t.test") <environment: namespace:ctest> # but expressing another function shapiro.test # returns more complete code? function (x) { DNAME <- deparse(substitute(x)) x <- sort(x[complete.cases(x)]) n <- length(x) if (n < 3 || n > 5000)

After a year my R programming style is still very "C like". I am still writing a lot of "for loops" and finding it difficult to recognize where, in place of loops, I could just do the same with one line of code, using "sapply", "lapply", or the like. On-line examples for such high level function do not help me. Even if, sooner or later, I am getting my R

Hello, I am trying to call a FORTRAN subroutine from R. The Fortran code is @: http://lib.stat.cmu.edu/apstat/206 It performs a bivariate isotonic regression on a rectangular grid (m X n) matrix. I used the g77 compiler and successfully created a dll file and it also loads successfully from R. But somehow the programs fails to run properly. (I do get the correct result when I compile the

Hello, I have a C++ code for 2D random walks with traps and I want to convert it in a R code with its syntaxs, can anyone help??????? It's easy for me to adapt the body but I want help with the beginig (variable declaration) and th end exporting the output to a file ( like write.table() or sink() ) Thank you... #include <iostream> #include <math.h>#include

Hello, here is a working reproducible example which crashes R using kmeans or gives empty clusters using the nstart option with R 15.2. library(cluster) kmeans(ruspini,4) kmeans(ruspini,4,nstart=2) kmeans(ruspini,4,nstart=4) kmeans(ruspini,4,nstart=10) ?kmeans either we got empty always clusters and or, after some further commands an segfault. regards, Detlef Groth ------------ [R] Empty

> library(Kendall) > Kendall(1:3,1:3) WARNING: Error exit, tauk2. IFAULT = 12 <<<<<< tau = 1, 2-sided pvalue =1 I believe Kendall tau is well-defined for this case and the reported value is correct; isn't it a bug to give a warning? (And if, e.g., the pvalue is not well-defined in this case, wouldn't it be better to return NA or NaN or something?) Also,

aggregate.ts does not behave in the same way as the equivalent method aggregate.rts in S-PLUS. In particular it - changes the start of the time series - tends to have a length which is 1 shorter For example: R> x <- ts(1:10) R> aggregate(x, nfreq=0.5, FUN=min) Time Series: Start = 2 End = 8 Frequency = 0.5 [1] 2 4 6 8 S> x <- rts(1:10) S> aggregate(x, nf=0.5, fun = min) [1]

From: Maros Zatko <mzatko@redhat.com> --- fish/rl.c | 158 ++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++ fish/rl.h | 32 +++++++++++++ 2 files changed, 190 insertions(+) create mode 100644 fish/rl.c create mode 100644 fish/rl.h diff --git a/fish/rl.c b/fish/rl.c new file mode 100644 index 0000000..bb8fd62 --- /dev/null +++ b/fish/rl.c @@ -0,0 +1,158 @@ +/* guestfish -

Dear McLeod, I am a PhD student working with the Kendall package in R. I used this package to make the so-called TauKr matrix correlation analysis (Hemelrijk, 1990). However, lately I've been getting this warning: Error exit, tauk2. IFAULT = 12. From the forums I've been looking into, I noticed that when this warning pops up, it means that the p-value of the correlation may be wrong

Hello everybody, I observed some strange behaviour of R (0.64.1/Linux 2.0.25) when trying to use FORTRAN-libraries. I first downloaded "ALGORITHM AS 304" from StatLib. All I changed at the code, was to substitute "REAL"-routines by "DOUBLE PRECISION". Compiling went o.k.: > g77 -fpic -O2 -shared -o rtest.o fisher.f ("g77 --version" gives 2.7.2.1)

Hello, I am an R newbie, trying to use lsoda to solve standard Lotka-Volterra competition equations. My question is: how do I pass a parameter that varies with time, like say, phix <- 0.7 + runif(tmax) in the example below. # defining function lotvol <- function(t,n,p){ x <- n[1]; y <- n[2] rx <- p["rx"]; ry <- p["ry"] Kx <-

Hi. It seems that nearly every time I try to use factanal I get the following response: >faa2db1<-factanal(mretdb1,factors=2,method="mle",control=list(nstart=25)) Error in factanal(mretdb1, factors = 2, method = "mle", control = list(nstart = 25)) : unable to optimize from these starting value(s) > In the case cited above, mretdb1 is synthetic data created

I''ve created a simple script (below) which calls the effect.appear script in order to make a group of items appear at the page load. I would like to have the images randomly appear at different times; e.g. the 3rd image might start appearing 2 seconds after load, the 6th image immediately after load, the first image 1 second after load, etc... How can I achieve this affect? ---