similar to: strange loss of accuracies under Win NT

Hi R-users,   I have this code here: library(numDeriv)   fprime <- function(z) { alp  <- 2.0165;   rho  <- 0.868;   # simplified expressions   a      <- alp-0.5   c1     <- sqrt(pi)/(gamma(alp)*(1-rho)^alp)   c2     <- sqrt(rho)/(1-rho)   t1     <- exp(-z/(1-rho))   t2     <- (z/(2*c2))^a   bes1   <- besselI(z*c2,a)   t1bes1 <- t1*bes1   c1*t1bes1*t2 }   ## Newton

Hi r-users,   I hope somebody can help me to understand the error message.  Here is my code; ## Newton iteration newton_gam <- function(z) { n   <- length(z)   r   <- runif(n)   tol <- 1E-6   cdf <- vector(length=n, mode="numeric")   fprime <- vector(length=n, mode="numeric")   f   <- vector(length=n, mode="numeric")     for (i in 1:1000)   {

Does anyone know how to read SPSS .por (ie Portable) files? The foreign package only deals with SPSS .sav files and not with those with a "por" extension. Thanks, José Mr José Luis Iparraguirre D'Elia Senior Research Economist Economic Research Institute of Northern Ireland Floral Buildings 2 - 14 East Bridge Street Belfast BT1 3NQ Northern Ireland United Kingdom

I am having difficulty installing the package "lattice". First I tried downloading it from the CRAN site (in the normal way) : the message said it worked but when I typed library I got an error message ("there is no package"). Second I tried installing it from the local zip. I have reproduced the result below. utils:::menuInstallLocal() package 'lattice'

Dear all, Does anyone know if the PLM package (to run Panel Data Analysis) accepts quarterly data? The package vignette and documentation only use annual data -and the only time index available seems to work for years. José Mr José Luis Iparraguirre Senior Research Economist Economic Research Institute of Northern Ireland 2 -14 East Bridge Street Belfast BT1 3NQ Northern

R-helpers I have 120 small Excel sheets to read and I am using library(xlsReadWrite): one example below. I had hoped to read sheets by looping over a list of numbers in their name (eg Book1.xls, Book2.xls, etc). I thought I had seen examples which used eval-parse-paste in this way. However, I have not been able to get it to work.. 1. is this a feasible approach? 2. if not

Hi all - I am trying to tune an SVM model by optimizing the cross-validation accuracy. Maximizing this value doesn't necessarily seem to minimize the number of misclassifications. Can anyone tell me how the cross-validation accuracy is defined? In the output below, for example, cross-validation accuracy is 92.2%, while the number of correctly classified samples is (1476+170)/(1476+170+4) =

Ape is an R package for "analyses of phylogenetics and evolution". The first version (0.1) has been released on 27 August 2002 and is available on CRAN. >From the 'Description' file of version 0.1: Ape provides functions for reading, and plotting phylogenetic trees in parenthetic format (standard Newick format), analyses of comparative data in a

At 11:57 25/10/00 +0100, Brian Ripley wrote: >> Date: mer., 25 oct. 2000 12:38:55 +0200 >> From: Emmanuel Paradis <paradis@isem.univ-montp2.fr> > >> I think it would be nice to have par(ask=T) set by default in termplot(), >> like it is in plot.lm(). > >Well, it isn't really the default in plot.lm, the default for `ask' being > >interactive()

Hi, I am using lmer (from the package lme4) to predict a binary response variable (REL) from a bunch of fixed effects and two random effects (Speaker_ID and NPhead_lemma): fit <- lmer(REL ~ SPEAKER_GENDER + log(SPEECHRATE) + SQSPEECHRATE + ..... + (1|Speaker_ID) + (1|NPhead_lemma), family="binomial", data=data.lmer, method="Laplace", model=T, x=T) I

Hi, I have a TrixBox install with a Sangoma A200 and 4 FXO ports, there are 3 BT lines connected directly to these ports. One of the lines has BT FeatureLine Compact and this is the line I am having problems with, the other 2 lines are working perfectly, detecting CID, answering incoming calls and placing external calls via SIP devices. I am receiving a error log entry: chan_zap.c:

There seems to be a bug with var() when the argument is a vector with exactly three values of 1e30 (or close to this value). This does not happen with twice, four (or more) times this value, or another value. > var(rep(1e30, 3)) [1] 2.971056e+28 > var(rep(1.2e30, 3)) [1] 2.971056e+28 > var(rep(0.9e30, 3)) [1] 2.971056e+28 > var(rep(0.8e30, 3)) [1] 0 > var(rep(1e29, 3)) [1] 0 >

Does anyone know how to retrieve the residuals after running an ERS unit root test on a time series (urca package)? There are no 'residuals' or 'res' (or similar) slots in the ur.ers-class, but a plot instruction - ie plot(ur.ers(....)) renders amongst other things a residual plot. If an augmented Dickey-Fuller test is run, an object of ur.df-class is generated which contains a

I am using the function prop.test (base package), which returns a list with class "htest". All I want to do is to assign one of its values to a variable, but I do not want R to print the results and added warning message whenever I invoke the function. How can I prevent R from printing on using prop.test? Regards, Jos? Luis Mr Jos? Luis Iparraguirre Senior Research Economist Economic

I could not find a package to run panel data unit root tests in R (there's a STATA module, though - PANELUNIT- and routines for Sarno and Taylor's MADF test, Levin-Lin-Chu's test, Im-Pesaran-Shin's Test, etc, to do this). Hence, I am toying with the idea of having a go at writing up one for R, but would like to know whether anyone is already working on this. José Mr

Hi: I am a newer in using R for data mining, and find the "e1071" pakage an excellent tool in doing data mining work! what frustrated me recently is that when I using the function "svm" and using the "cross=10" parameters, I got all the "accuracies" of the model greater than 1. Isn't that the accuracy should be smaller than 1? so I wander how, the

Hi, I have a Data Set x and I want to check with a Kolmogorow-Smirnow-Test, if x comes from a Subbotin Distribution, whose density function is: function(x,location,scale,tail) # Exponential power (=Subbotin) { const<- 2*scale*tail^(1/tail) *gamma(1+1/tail) z<- (x-location)/scale exp(-1/tail*abs(z)^tail)/const } How can i do this? Thank you --

Hello everybody! I've got a network with a working Windows2000 Server (acting as a file server and being responsible for user authentification) and some Win2000 Professional Clients. I want to do the following: I want to integrate other Workstations to the existing Domain and I want to integrate a Linux computer. Therefor I want to use Samba to share disk space with the Win Clients. But there

Hi, I think it would be nice to have par(ask=T) set by default in termplot(), like it is in plot.lm(). Best, Emmanuel Paradis -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !)

Dear R-users, Thanks to Friedrich Leisch, my document "R pour les d?butants" is available on CRAN in the contributed documentation section: http://cran.r-project.org/doc/contrib/Rdebuts.pdf. It is written in French. "R pour les d?butants" ("R for beginners") gives a starting point for people with no experience with R (or S). I tried to explain the basics in the