similar to: 'postscript' command within a function

> hi all, > I am using optim. I am getting the following error message: > > Error in optim(par = start.vals[, h], fn = post.func.pois, gr = post.grad. > pois, : > L-BFGS-B needs finite values of fn > > If I look at optim typing '> optim' it seems that the error comes from > inside .Internal(optim), so I wonder how can I see the code for .Internal(

hello, I am using LDA and QDA function of MASS library. I understand Fisher LDA is a method non-probabilistic in nature, so I wonder what happens when I try to predict my test set examples as in: > fit <- lda(labels~., data=train.table, prior=c(.5,.5)) > pred <- predict(fit, data=test.table, prior=c(.5,.5)) Specifically I ask this because in my problem there are 700 examples

Hello all, I am using libraries (mva,cluster) to produce dendrograms. With 1000 examples the dendrogram gets too crowded, and i am wondering whether there is an option (which i cannot find) to set the number of leaf nodes, like in matlab, and return the plot and the assignment map examples -> leaf nodes. Any suggestion is appreciated. Thanks Edo

Dear R users, I'm having trouble using aggregate() and would greatly appreciate your advice. I am using R 1.5.1 on Windows 2000. I want to call my function in the following way extract.data(x=dat[, "Age", "Year"]) where extract.data() uses aggregate() to count the number of cases for each combination of "Age" and "Year". I've defined

Hi, I am using the following R version: > R version 2.4.1 (2006-12-18) > Copyright (C) 2006 The R Foundation for Statistical Computing > ISBN 3-900051-07-0 I believe I found a bug in the 'hist' function, when 'probability=TRUE'. I looked in the archives and I came across problems with the 'hist' functions (e.g., bug PR# 944, posted in 2001), however, a

Dear users, is anybody of you aware of a R command to perform laplace transform or even its inversion? Thank you very much. Luca

My computer is a pentium 4 running at 2.4 GHz. My R is 1.7.1 I have written a program in R that calls C++. The program spends most of the time in C++ ( > 90% ). R basically deals with output and input. How slower can this be compared with the program I would get from rewriting everything in C++? Thank you.

Hi, I am running a loop to plot multiple plots. In s-plus, it shows multiple pages in the graphic window to allow checking on each plot. but in R, the next plot always overwrite the previous one, so i can only have the last plot produced, is there a way to have multiple pages in the graphic window just like S-plus does? Thanks

installed from the binary image (http://cran.r-project.org/bin/ macosx/) on an intel mac 10.4.10 throws three warnings: ----------------- R version 2.5.1 (2007-06-27) Copyright (C) 2007 The R Foundation for Statistical Computing ISBN 3-900051-07-0 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or

Hello there friends & foes, I've had if with M$ windows! Next week I'll sabotage my shop (how 'bout password protecting dualboots on WS) and force LINUX down my users' throats ;-) (oh, don't I like this tone... makes me feel Billy G.) Oops, showstopper! ACL enabled network filesystem, where is it? I can't go with alpha code nfs4, so is there a patch to the smb

hallo i'm running mepis (debian unstable) and if i install wine from the regular (non-wine specific) deb repositories, I get wine 20041019. if I install from wine.sourceforge, I get 200501011. with the former, the application seems to lock up much more often. with the latter, the ability to play the movie (e.g. to see which audio stream to include) is disabled. i see from this screenshot

Hi, Suppose I have a multidimensional array: tmp <- array(1:8, c(2,2,2)) is there a function out there that, given a one-dimensional array index, will return the separate indices for each array dimension? for instance, tmp[8] is equivalent to tmp[2,2,2]. I'd like to derive the vector (2,2,2) from the index 8. thanks, Brad Buchsbaum

I added the winehq debian repository to apt-get | synaptic (I am running mepis 3.3) and would like to install those rather than the debian official unstable builds but since the winehq repository is also marked as part of 'unstable' distro I can't find a way to distinguish the wine / winelib provided by winehq from the official ones (which would need the many winelib-* packages)

I have written the following function to calculate the weighted mean difference for univariate data (see http://www.xycoon.com/gini_mean_difference.htm for a related formula). Unsurprisingly, the function is slow (compared to sd or mad) for long vectors. I wonder if there's a way to make the function faster, short of creating an external C function. Thanks very much for your advice. gmd

I am participating in a project in the university, and appeared a doubt: if can I call one function from R as "solve" from C code. The intention of the research is to create a DLL. Is it possible and where can I find an example? Sorry for bad english! ___________________________________________________________________________________ Acesse nosso portal www.click21.com.br Porque

I've seen that there's an history for graphics. How can I save it manually? Can I save more than 4 plots? I'm using R1.6.2 under win98 [[alternative HTML version deleted]]

Hey, All Now I have a data set which is n-dimensional. And want to plot the Scatter Plot Matrix which is n by n. Does R have such function to achieve this? Thanks for your point. Fred

I am brand new to R and would like to know if there is a way to be able to create, save and retrieve a program in R like a script object in S-Plus. I am just not sure how to do this. I'm sure there are many ways to do this ranging from simple to complex, but the simplest way would be best for me. Thanks, in advance... __________________________________________________ Y! Web Hosting - Let

How do I get non-squared correlation coefficient in some more sensible way than sqrt(summary(lm(y~x))$r.squared)? Thanks Matej -- Matej Cepl, matej at ceplovi.cz, Finger: 89EF 4BC6 288A BF43 1BAB 25C3 E09F EF25 D964 84AC 138 Highland Ave. #10, Somerville, Ma 02143, (617) 623-1488 The difference between death and taxes is death doesn't get worse every time Congress meets -- Will

Dear List Members I'm using R to create a trellis plot using the library of Pinheiro & Bates (trellis graph for grouped data). Here is my question: How can I change the format of the axis (x or y) in terms of the number decimal points that should appear on the plot? In other words, I have that for a set of plots in the y-axis the values appear as 0.2,0.3, etc and for another set they