similar to: Trigamma function

Hi, Is there any R library that is capable of handling polygamma function (Hurwitz zeta function also works)? I am aware of digamma(0 and trigamma(), but could not find more advanced versions. I'd appreciate any help. Hakan Demirtas

Hi, 2 questions: Question 1: example of what I currently do: for(i in 1:6){sink("temp.txt",append=TRUE) dput(i+0) sink()} x=scan(file="temp.txt") print(prod(x)) file.remove("C:/R-2.5.0/temp.txt") But how to convert the output of the loop to a vector that I can manipulate (by prod or sum etc), without having to write and append to a file? Question 2: >

Hello, I was wondering if the functions deriv3(), deriv() etc. could be extended to handle psigamma() and its special cases (digamma(), trigamma() etc.). From the error message it seems that 'psigamma' needs to be added to the derivatives table. This might be easy since psigamma() has a deriv argument. Additionally, this error message is also obtained when requesting for the Hessian of

Hello, I was asked to try the following code on R, gamma.mles function (xx,shape0,rate0) { n<- length(xx) xbar<- mean(xx) logxbar<- mean(log(xx)) theta<-c(shape0,rate0) repeat { theta0<- theta shape<- theta0[1] rate<- theta0[2] S<- n*matrix(c(log(rate)-digamma(shape)+logxbar,shape/rate-xbar),ncol=1) I<- n*matrix(c(trigamma(shape),-1/rate,-1/rate,shape/rate^2),ncol=2)

Hi, I'm currently in the process of writing an R-installation SOP for my company. As part of that process I'm using the recommendations from the 'R Installation and Administration' document, section 3.2, "Testing an installation". This is done on an XP machine, using the latest binary of 2.11.0. The binary is downloaded and then installed from the installer. I then

Hey, R-listers I was recommended to try using smooth.spline function for estimating 2-Dimensinal curve given a data set. So will you please tell me where to get this R function? Or which package provides this function? Thanks for your point. Fred

Hey, I want to draw several plots sequently, but have to make them dispaly in one figure. So how to achieve this? Thanks. Fred

Hi, R-listers I tried to plot several graphs in a sigle x-y coordinate settings, like the following: |(y) s | ****** s | ***** s | sssssssssssssssssss |_______________________________(x) where "*" and "s" denote two diffrent plots. However, when I used plot(data1); % data1 is the data points of "*"

Hey, All In principal component analysis (PCA), we want to know how many percentage the first principal component explain the total variances among the data. Assume the data matrix X is zero-meaned, and I used the following procedures: C = covriance(X) %% calculate the covariance matrix; [EVector,EValues]=eig(C) %% L = diag(EValues) %%L is a column vector with eigenvalues as the elements percent

Dear R-listers, I have a dxr matrix Z, where d > r. And the product Z*Z' is a singular square matrix. The problem is how to get the left inverse U of this singular matrix Z*Z', such that U*(Z*Z') = I? Is there any to figure it out using matrix decomposition method? Thanks a lot for your help. Fred

Hey, all Now I am using EM algorithm to do some optimization. Within that E-step, I have to calculate some multivariate integration given some parameter values Theta and function form of a probability density function f. So I want to know if there are some package in R to do such numerical integration given such function form and parameter values? Thanks for your kind support on this. Have a

Full_Name: Chanseok Park Version: 1.8.1 OS: linux-gnu Submission from: (NULL) (130.127.112.183) digamma with any negative value does not give a right answer. It gives -1.797693e+308 for any negative arguments. For example, digamma(-1.1) gives -1.797693e+308. The right answer should be 10.15416 This bug can be easily fixed by using the following digamma identity. digamma(x) = digamma(1-x) -

Assume F(x) is the cdf of stardard normal c.d.f, and want to get the integration of F(x)^2 over (-infinite, +infinite). So whats the value of this integration? And is there some function to achieve this? Thanks. -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info",

Hi, I am running R 2.13.0 on a Windows 7 machine. I ran the script: testInstalledBasic('devel') and received the following warning message: running tests of consistency of as/is.* creating ?isas-tests.R? running code in ?isas-tests.R? comparing ?isas-tests.Rout? to ?isas-tests.Rout.save? ...running tests of random deviate generation -- fails occasionally running code in

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

Summary: proc mixed vs. lme The objective of this summary is to help people to get more familiar with the specification of random effects with proc mixed or lme. Very useful are the examples of Ramon Littell's book: "SAS System for Mixed Models (1996)" (http://ftp.sas.com/samples/A55235) The same data set's are kindly made available by Douglas Bates in the

Dear All, In the N-dimensional space, give a data point A and a curve f, how to write the explicit expression for calculating the minimal distance between A and f? Or have to use some nonlinear optimization method to calcualte it? Thanks for your point. Fred [[alternative HTML version deleted]]

Hi, I have just seen Stefan Steinhaus' web page : http://www.uni-franfurt.de/~stst/ncrunch.html I think it would be nice to include "R" as well. I have taken Forrest Young's email on stat-lisp list and changed the stuff for "R" :) Here it is: (someone please check this so we can also send it to Stefan Steinhaus.

I've rolled up R-2.6.0.tar.gz a short while ago. This is a development release which contains a number of new features. In particular, the handling of data with a large number of identical strings should be more memory-efficient. Also, a number of mostly minor bugs have been fixed. See the full list of changes below. You can get it (later today) from

I've rolled up R-2.6.0.tar.gz a short while ago. This is a development release which contains a number of new features. In particular, the handling of data with a large number of identical strings should be more memory-efficient. Also, a number of mostly minor bugs have been fixed. See the full list of changes below. You can get it (later today) from