similar to: UseMethod infelicity

Hi all, I´m re-starting (as my name indicates) my little knowlegde of CRAN R. I made this function time before but I don´t know where is the error because nothing appears as an error but the histogram plot doesn´t appear. Should I install some special library to run sapply? pru<-function(){ randz<-matrix(rnorm(200000),100,2000) H<-matrix(0,100,2000) for (j in 2:2000){ for (i in

Hi all, I want to make the following: I want to run a linear regression on each column of a matrix "estima" on the correspondent column on the matrix "estima2". You see I want to regress estima[,1] on estima2[,1] this way to all columns.... At the same time I want to make a regression adding each time a new observation. You see, the first regression will regress only one

Using > methods("plot") [1] plot.Date* plot.HoltWinters* plot.POSIXct* [4] plot.POSIXlt* plot.TukeyHSD plot.acf* [7] plot.data.frame* plot.decomposed.ts* plot.default [10] plot.dendrogram* plot.density plot.ecdf [13] plot.factor* plot.formula* plot.hclust* [16] plot.histogram* plot.isoreg* plot.lm [19] plot.medpolish*

Hi, I'm using R to do a time series analysis. In the model, I use the lags of some variables. such the lags of the variables have different length, I just can't use them directly in the lm function. Intuitively, I feel that "subset" might be useful, but I do not know how to use it. Anyone can give me an example syntax? Thanks. Harry [[alternative HTML version deleted]]

Een ingesloten tekst met niet-gespecificeerde tekenset is van het bericht gescrubt ... Naam: niet beschikbaar Url: https://stat.ethz.ch/pipermail/r-help/attachments/20070820/920567a6/attachment.pl

Hello! I'd like to create an assets-based economic indicator using data from a national household survey. The economic indicator is to be the first principal component from a principal components analysis, which (given the source of the data) I believe should take in consideration the sampling weights of the observations. After running the PCA with svyprcomp(), from the survey package, I

Any better (more efficient, built-in) ideas for computing coef[1]+coef[2]*x+coef[3]*x^2+ ... than polynom <- function(coef,x) { n <- length(coef) sum(coef*apply(matrix(c(rep(x,n),seq(0,n-1)),ncol=2),1,function(z)z[1]^z[2])) } ? Ben -- 318 Carr Hall bolker at zoo.ufl.edu Zoology Department, University of Florida http://www.zoo.ufl.edu/bolker

Dear list, Please find below an example of odd behaviour of the by function. It occurs both under GNU/Linux R 2.6.2 and Windows R 2.7.0alpha. Respective sessionInfo()'s are given below. I hope I do not overlook anything. testFactor <- factor(sample(LETTERS[1:6], size = 42, replace = TRUE)) testMatrix <- matrix(rnorm(42 * 6), nrow = 42) testDf <- as.data.frame(testMatrix) by(data

format(c(1:2, NA)) gives the last value as "NA" rather than preserving it as NA, even if na.encode = FALSE (which does the 'expected' thing for character vectors, but not numeric vectors). This was already brought up in 2008 in https://bugs.r-project.org/show_bug.cgi?id=12318 where Gregor Gorjanc pointed out the issue. Documentation was added and the bug closed as

utils::citation() offers the option auto=(NULL|TRUE) to obtain the citations for a package from the package DESCRIPTION file or a CITATION file, as described below: |auto| a logical indicating whether the default citation auto-generated from the package 'DESCRIPTION' metadata should be used or not, or |NULL| (default), indicating that a 'CITATION' file is used if it exists, or

In theory, dlnorm(x, ...) == dnorm(log(x), ...)/x, x>0. Unfortunately, when sd=0, dlnorm and plnorm return NaN, while dnorm returns (if(x != mean)0 else Inf) and pnorm returns (if(x<mean)0 else 1). [A numerical optimization, maxLik{maxLik}, reported the NaNs for me.] help('dlnorm') says, "dlnorm is calculated from the definition (in ?Details?). [pqr]lnorm are based on the

Dear list, "+" (and "-") being defined for difftime class, I expected mean() to return something sensible. This is only half-true : > mean(c(1:5, 5:1),na.rm=TRUE) [1] 3 > mean(as.difftime(c(1:5, 5:1),unit="mins"),na.rm=TRUE) Time difference of 3 mins Fine so far. However : > mean(c(1:5, NA,5:1),na.rm=TRUE) [1] 3 > mean(as.difftime(c(1:5,

Hello, The main reason for resetting the tagbuf in `print.default()` and other entry points to the print routine is that it is currently not reset on exit. Creating a context to reset it on exit to its last value might work. This should be done in the entry points rather than in print-value-rec though, since callers of the latter might write to the tagbuf. Another solution to this problem is

Hello: sort(c('A', 'b', 'C')) seems to produce different answers in R interactive than in "R CMD check", at least under both Fedora 13 and Windows 7 with Windows 7 sessionInfo() copied below: In interactive, the result is c('A', 'b', 'C'); with R CMD check, it is c('A', 'C', 'b'). This produced

Hi, I'm having problems with the Zelig package - when using the below R displays the follwing message (I'm running R i386 2.15.3 for Windows and have updated all the Zelig packages): z.out<-zelig(Surv(psurv2, pcens2) ~ ren_sup3 + age, data=data_urgent, model="coxph") ** The model "coxph" is not available with the currently loaded packages, ** and is not an

Hello, I am writing as an administrator, not as an R user, so forgive me if I am not completely knowledgeable about R. I have a user who is creating an R package for windows from a Linux environment using the crossbuild environment by Jun Yan and A.J. Rossini. The packages she generated worked fine until she tried to install in R 1.9.1 for Windows. Now when she installs with

Hi Folks, According to the PDF file blogit.pdf in the Zelig documentation: "Use the bivariate logistic regression model ["blogit"] if you have two binary dependent variables (Y1,Y2), and and wish to model them jointly as a function of some explanatory variables. Each pair of dependent variables (Yi1,Yi2) has four potential outcomes, (Yi1=1,Yi2=1), (Yi1=1,Yi2=0),

Currently S3 objects nested in generic vectors cause the tag buffer to be reset.? This feels sub-optimal for those objects that don't have a print method: > list(a=list(b='hello')) $a $a$b???????????????? ### <<<< notice "$a$b" [1] "hello" > list(a=structure(list(b='hello'), class='world')) $a $b?????????????????? ###

Dear R-user: After the fitting the Tobit model using zelig, if I use the following command then I can get the regression coefficents: beta=coefficients(il6.out) > beta (Intercept) apache 4.7826 0.9655 How may I extract the "Naive SE" from the following output please? > summary(il6w.out) Call: zelig(formula = il6.data$il6 ~ il6.data$apache, model =

hello all I am using the R package Zelig for some tobit regression with robust standard errors. I have got R version 2.9.2 (2009-08-24) and Zelig Version: 3.4-5 when i do demo(robust) It ends like this way data(coalition) > # Fit the model with robust standard error > user.prompt() Press <return> to continue: > z.out3 <- zelig(Surv(duration, ciep12) ~ polar + numst2 +