search for: stirlings

Hi anyone coded up Stirling numbers in R? [I need unsigned Stirling numbers of the first kind] cheers Robin -- Robin Hankin Uncertainty Analyst National Oceanography Centre, Southampton European Way, Southampton SO14 3ZH, UK tel 023-8059-7743

Hi there - I am using R1.5.1 on WinNT and the latest MASS (Venables and Ripley) library. Running the following code: >minimod<-glm(miniSF~gtbt*f.batch+log(mxjd),data=gtbt,family="binomial") >summary(minimod,cor=F) Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) 0.91561 0.32655 2.804 0.005049 ** gtbtgt 0.47171

I have some big combinations like: 4444444444444444444444444444 choose 784645433 Can R compute these? Is there any package that does stirlings approximation in R? -- Thanks, Jim. [[alternative HTML version deleted]]

Hello all. I have two small questions in one post, for the sake of brevity. 1. I have some objects that I want to delete. I have the line: rm (c (channelheader, paste ("channel", 1:3, sep=""))) I have tried a few variations, including list=, but cannot figure it out. In SAS, I can use a ':' as a wildcard. Is there any equivalent in R? 2. Is there any possible was to

Oliver Stirling wrote: > Oliver Stirling wrote: > >> My icecast server is behind a http proxy server, does anyone know of a >> program I can use to get icecast to send data to that will forward it >> to the yp dir? (I'm on linux btw) >> >> Oli >> >> --- >8 ---- >> List archives: http://www.xiph.org/archives/ >> icecast project

Is there any way to read these files (standard .mat files, created by matlab version 7 onwards are compressed)? I know that R.matlab doesn't read them (it even says in the file MatlabServer.m "Matlab v7 saves compressed files, which is not recognized by R.matlab's readMat()" (lines 47-8)). I know I should be able to make R call Matlab and transfer data (not that I managed to

Hi Daniel, > Do you have experience how JEdit compares to other > R-capable editors like e.g. WinEdt? Does it offer something like WinEdt's > integration with the Windows R-gui? I have never tried out WinEdit but after reading your message I took a quick glance at their homepage. It would not be useful for JEdit to implement the provided interaction with R-GUI since JEdit has to

My icecast server is behind a http proxy server, does anyone know of a program I can use to get icecast to send data to that will forward it to the yp dir? (I'm on linux btw) Oli --- >8 ---- List archives: http://www.xiph.org/archives/ icecast project homepage: http://www.icecast.org/ To unsubscribe from this list, send a message to 'icecast-request@xiph.org' containing only

Dear R users, i try to integrate lgamma from 0 to Inf. But here i get the message "roundoff error is detected in the extrapolation table", if i use 1.0e120 instead of Inf the computation works, but this is against the suggestion of integrates help information to use Inf explicitly. Using stirlings approximation doesnt bring the solution too. ## Stirlings approximation lgammaApprox <- function(x) { 0.5*log(2*pi)+(x-(1/2))*log(x)-x } integrate(lgamma, lower = 0, upper = 1.0e120) integrate(lgammaApprox, lower = 0, upper = 1.0e120) > integrate(lgamma, lower = 0, upper = 1.0e120) 1.3740...

Hi One of my packages needs a look-up table of pre-calculated numbers in the data directory. I would like to have the matrix as large as possible. What is the largest size matrix that would be an acceptable datafile in an R package? [ The table is a square, upper triangular matrix consisting of logs of Stirling numbers calculated by Maple. As discussed on the List a few days ago (thanks

On 17 Jun 2003 at 14:18, Oliver Stirling wrote: > On a slightly different note I was wondering if it was possible to add > to the yp directory on a different port than http? My station is a > student radio station in the UK and we are behind the uni firewall and > as far as I know I can't get icecast to use the proxy. > > However with shoutcast we can report to port 666

Hi all, I have created a new edit mode file for the editor JEdit (http://www.jedit.org) based on the file from Zed Shaw. The difference between his and my file is that mine can be used as a spell checker. It does not highlight all functions but provides syntax highlighting for the base package (keyword types 1 and 2) and 26 additional packages (boot, class, cluster, ctest, eda, foreign, grid,

Hi all, This is probably a very simple question but I cannot figure out how to do it. I run the fourthcorner method with my data and would like to adjust the p values for multiple comparisons using Holm correction. When I run the fourthcorner I obtain the results in yellow. What do I need to do to be able to aply the Holm correction to those p values? > library(ade4) > four1 <-

I tried to intall R in a iMac with Mac OS E1- 9.0.4 (Mac OS ROM VE-4.6.1). When I tried to run R for the first time I got a message about "Carbonlib-- CF RangeMake" not beeing found. Please, could you help me to tackle the problem. Thanks very much. Barahona -- Enrique Barahona Fernandez ESTACION EXPERIMENTAL DEL ZAIDIN, CSIC Profesor Albareda, 1 18008, GRANADA (Spain) Tel. + 34

Stefan Neufeind wrote: > On 17 Jun 2003 at 14:18, Oliver Stirling wrote: > > >>On a slightly different note I was wondering if it was possible to add >>to the yp directory on a different port than http? My station is a >>student radio station in the UK and we are behind the uni firewall and >>as far as I know I can't get icecast to use the proxy. >>

Hi guys, I want to get Rails to run optimally on Fedora Core 4 + Apache 2.0.54. Can anyone point me in the right direction please? Best Regards, Ashish. Senior Consultant SAP For Business Ltd., 37 Springfield Road, Stirling, Scotland, United Kingdom FK7 7QN. +44-774-810-7492. -------------- next part -------------- An HTML attachment was scrubbed... URL:

Dear R community, I was using R 2.11.1 without internet connection on Windows XP and whenever I type ?mean for example, R would freeze and crash out... Is this something that can be fixed? I would like to use the internel help file if possible... Thanks. ########################################### Assistant Professor Steve Su School of Mathematics and Statistics Faculty of Engineering,

Hi Folks! I have a question regarding an issue on which I´ve read a few threads but can´t resolve --> Labelling every other tick mark on the x-axis using factors. Here´s an example: db<-data.frame(date=as.factor(c(1:50)),var=rnorm(50)) barplot(db$var,names.arg=levels(db$date),las=2,cex.axis=0.8,cex=0.8) This labels all the tickmarks in my barplot, which is the default behavior. I would

We've been going through some growing pains fumbling through porting some custom C modules that used to live under apache 1.3 into the 2.0 spec, even though none of us really considers ourselves a 'C developer'. We've basically got them all working now under the 2.0 API, and the question was posed about whether or not we should consider attempting to upgrade everything to run

I have found a strange scoping problem in Sweave. The following Rnw file doesn't produce the same output in Sweave as it does if I produce an R file using Stangle and execute that: \documentclass[12pt]{article} \begin{document} <<R>>= election <- data.frame(A=1:3, B=9:7, C=rep(0,3)) partytotal <- rep(0, ncol(election)) for (i in 1:ncol(election)) { partytotal[i] <-