[R] confidence interval around Aalen-Johansen if no events

2023-09-06 Thread Andreas Leha
(using the bpcp package). Many thanks in advance! Best, Andreas PS: Some code that illustrates my problem: - library("dplyr") library("survival") ## read data sdat <- structure(list(time = c(2, 189, 182, 2, 184, 179, 179, 159, 18, 177, 177, 182,

[R] unused argument(s) (Header = 1) help!

2023-08-09 Thread Andreas Noviyanto
.out, file="lans", + set="01", ext="xlsx", header=1) Error in ABE(path.in = path.in, path.out = path.out, file = "lans", : unused argument (header = 1) > ABE(path.in=path.in, path.out=path.out, file="lans", + set="02", ext

Re: [R] Use of ellipsis

2022-05-29 Thread Andreas Matre
Thank you very much Ivan and Bert! I used the eval(substitute()) workaround suggested by Ivan and it worked perfectly. Andreas Matre __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help

[R] Use of ellipsis

2022-05-26 Thread Andreas Matre
think that's relevant, since as far as I can tell, there should be no non-standard evaluation here. [1] https://stat.ethz.ch/pipermail/r-help/2010-February/228881.html Andreas Matre __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more

Re: [R] PSOCK cluster and renice

2017-12-03 Thread Andreas Leha
Hi Henrik, Thanks for the detailed in fast reply! My guess would be that the confusion comes from the different use of nice and renice. The workraund you provided work fine! Thanks a lot. Best, Andreas Henrik Bengtsson writes: > Looks like a bug to me due to wrong assumptions ab

[R] PSOCK cluster and renice

2017-12-03 Thread Andreas Leha
: > cl <- parallel::makeCluster(2, renice = 19) nice: ‘+19’: No such file or directory Kind regards, Andreas __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guid

[R] greport

2017-10-17 Thread Andreas Betz
package. Thank you Andreas This email message is for the sole use of the intended r...{{dropped:11}} __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http

Re: [R] encoding/locale problem with ssh -X

2017-06-23 Thread Andreas Leha
Hi Paul, Thanks for following this up! It used to be R version 3.3.1 I updated to R version 3.4.0 Now everything seems to work! Many thanks! Best, Andreas On 23/06/17 03:02, Paul Murrell wrote: > Hi > > What version of R do you have (on the remote machine) ? > > I can

[R] encoding/locale problem with ssh -X

2017-06-21 Thread Andreas Leha
ForwardX11 yes | ForwardX11Timeout 596h | IdentityFile ~/.ssh/id_rsa | IdentitiesOnly yes | ForwardAgent yes | ServerAliveInterval 300 `---- Thanks in advance for your help! Best, Andreas __ R-help@r-project.org mailing

[R] [R-pkgs] New Package: pinbasic - Fast and Stable Estimation of the Probability of Informed Trading (PIN)

2016-10-26 Thread Andreas Recktenwald
be a good starting point. Development version of the package is available at: https://github.com/anre005/pinbasic -- Diplom-Kaufmann Andreas Recktenwald Statistik & Ökonometrie Rechts- und Wirtschaftswissenschaftliche Fakultät Universität des Saarlandes Campus C3 1, Raum 2.06 66123 Saarbrück

[R] outlier labels incorrectly assigned with ggplot2 box plot

2016-10-20 Thread Andreas Nord
atedly and erroneously (and seemingly randomly). Labelling outliers by their numerical values (i.e. by changing 'subject' to 'variable1' in the 'safe_ifelse function) does not cause problems. I assume I am missing something obvious here -

[R] ifelse for creating discriminating variable based on two conditions

2016-10-14 Thread Andreas Nord
Dear list, Apologies for a likely naive question. I am trying to create a discriminating dummy variable using 'ifelse' based on conditions in two variables. Specifically, I want to assign levels in a new factor as '0' or '1' based on a user-defined cut off. I.e. something similar to: >da

[R] [R-pkgs] New package 'forestinventory: Design-Based Global and Small-Area Estimations for Multiphase Forest Inventories'

2016-07-11 Thread Hill Andreas
vities: We will publish a vignette illustrating the use of the various estimators as soon as possible. Best, Andreas Hill Alexander Massey --- ETH Z�rich Andreas Hill Forstl. Ingenieurwesen, Heinimann CHN K 75.1 Universit�tstrasse 16 8092 Z�rich Telefon:

Re: [R] why must a named colClasses in read.table be in correct order

2015-07-08 Thread Andreas Leha
Hi Henrik, Thank you very much for looking into this. And thanks for the patch! Yes, let's hope this is a typo that gets fixed. Regards, Andreas Henrik Bengtsson writes: > Thanks for insisting; I was wrong and I'm happy to see that there is > indeed code intended for n

Re: [R] why must a named colClasses in read.table be in correct order

2015-07-08 Thread Andreas Leha
cter"))) --8<---cut here---end--->8--- What am I missing? Best, Andreas On 09/07/2015 02:21, Henrik Bengtsson wrote: > read.table() does not make use of names(colClasses) - only its values. > Because of this, ordering is critical, as you noted. It shou

[R] why must a named colClasses in read.table be in correct order

2015-07-08 Thread Andreas Leha
colClasses = cclasses) ## <--- error read.table(textConnection(kkk), sep="\t", header = TRUE, colClasses = cclasses[order(names(cclasses))]) --8<---cut here---end--->8--- Thanks, Andreas

Re: [R] aov + pairwise.t.test: only selected pairs

2014-06-22 Thread Andreas Maunz
with pooled variance are not supported in multcomp: http://www.inside-r.org/packages/cran/multcomp/docs/contrMat Am I wrong? If not, how could I do repeated t-tests for specific contrasts? Thank you Andreas On Wed, Jun 11, 2014 at 10:28 AM, ONKELINX, Thierry < thierry.onkel...@inbo.be> wrote

[R] multcomp: Contrasts for Tukey

2014-06-18 Thread Andreas Maunz
I want to use the glht function. In the linfct argument, I provide my linear functions to be tested. I provide contrasts there, such as (for example) c("Group1 - Group2 = 0", "Group1 - Group3 = 0") to compare group no 1 against the other groups. I want to do Tukey's test in this fashion. Is it po

[R] aov + pairwise.t.test: only selected pairs

2014-06-11 Thread Andreas Maunz
does not allow me to select specific pairs, it always analyzes all pairs. How can I analyze specific pairs only? Thanks Andreas [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help

Re: [R] filename of current device

2014-05-19 Thread Andreas Leha
On 18/05/2014 18:08, Prof Brian Ripley wrote: > On 14/05/2014 23:04, Andreas Leha wrote: >> Hi all, >> >> how do I find out about the filename, the currently open (let's say pdf) >> device is writing to? >> >> If I find 'dev.cur()' returnin

Re: [R] filename of current device

2014-05-14 Thread Andreas Leha
Hi David, David Winsemius writes: > On May 14, 2014, at 3:44 PM, Andreas Leha wrote: > >> Hi David, >> >> Thank you for following up on that. >> >> David Winsemius writes: >> >>> On May 14, 2014, at 3:04 PM, Andreas Leha wrote: >&

Re: [R] filename of current device

2014-05-14 Thread Andreas Leha
Hi David, Thank you for following up on that. David Winsemius writes: > On May 14, 2014, at 3:04 PM, Andreas Leha wrote: > >> Hi all, >> >> how do I find out about the filename, the currently open (let's say pdf) >> device is writing to? >> >>

[R] filename of current device

2014-05-14 Thread Andreas Leha
for my poor organization is highly appreciated... Regards, Andreas __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

Re: [R] transition from depends to imports

2014-05-09 Thread Andreas Leha
Duncan Murdoch writes: > On 09/05/2014, 5:39 PM, Andreas Leha wrote: >> Hi Duncan, >> >> Thank you for your follow-up and fast response! >> >> Duncan Murdoch writes: >> >>> On 09/05/2014, 5:08 PM, Andreas Leha wrote: >>>> Hi all, >

Re: [R] transition from depends to imports

2014-05-09 Thread Andreas Leha
Hi Duncan, Thank you for your follow-up and fast response! Duncan Murdoch writes: > On 09/05/2014, 5:08 PM, Andreas Leha wrote: >> Hi all, >> >> I am having trouble to do a particular transition from depends to >> imports in one of my packages. >> >> Th

[R] transition from depends to imports

2014-05-09 Thread Andreas Leha
at it is usable? Many thanks in advance. Regards, Andreas __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

[R] fitting a family of curves

2014-04-23 Thread andreas betz
referral to packages. Thank you Andreas [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and

Re: [R] Problems with a R-packages

2014-01-24 Thread Andreas Rybicki
ackages’ > " If it does not take too much of your time, maybe you could have a look at it? AND: "lossDev" got no mark under "installed version", but that is maybe obvious after the a/m R-message. I think that I have to look for prof. help here in Germany. Anyway, thank

[R] Problems with a R-packages

2014-01-22 Thread Andreas Rybicki
tion-window, there appears no remark in the column "installed version", it is still blank. Next I tried to install "rjags" and "logspline". Same result regarding download. Even if this question is simple, any help highly appreciated. I am available via Skype too, to s

Re: [R] quickly extract response from formula

2013-11-01 Thread Andreas Leha
resp <- if (length(frm)==3) frm[[2]] else NULL > (or call stop(), or warning(), ...) Will do. Thanks. - Andreas > > Bill Dunlap > Spotfire, TIBCO Software > wdunlap tibco.com > > >> -Original Message- >> From: r-help-boun...@r-project.org [mailto:

Re: [R] quickly extract response from formula

2013-11-01 Thread Andreas Leha
Hi David, thanks for your quick answer! David Winsemius writes: > On Oct 31, 2013, at 1:27 PM, Andreas Leha wrote: > >> Hi all, >> >> what is the recommended way to quickly (and without much burden on the >> memory) extract the response from a formula? > >

[R] quickly extract response from formula

2013-10-31 Thread Andreas Leha
uot;) data("ALL") y <- pData(ALL)$sex x <- t(exprs(ALL)) mf <- cbind(as.data.frame(x), y=y) extractResponseFast <- function(formula, data) { y <- eval(as.symbol(as.character(formula)[2]), environment(formula)) y } extractResponseFast(y~., data=

Re: [R] No speed effect by using RcppArmadillo compared to R in matrix operations

2013-10-26 Thread Andreas Recktenwald
Hi, another option if you're using Linux AND an Intel processor would be linking R against Intel MKL (Math Kernel Library). Under Linux you can get a (free) non-commercial licence for it. Here I'm using an Intel(R) Core(TM) i5-3210M CPU @ 2.50GHz laptop processor with R 3.0.2 build with i

[R] Heatmap with mouseover

2013-10-05 Thread Andreas Maunz
Dear all, what is the recommended way to create heatmaps (or plots in general) with mouseover-effect in R, so that I can embed it into a web app? Options I have found so far are a) rCharts or b) sendplot. What are your experiences and recommendations? Thanks Andreas [[alternative HTML

[R] R function for censored linear regression

2013-09-24 Thread Andreas Wittmann
my idea is to use survreg from the survival package library(survival) s1 <- survreg(Surv(y, stat) ~ x1 + x2, data, dist='gaussian') summary(s1) my question is, is this the right approach for my aim? Is it right, that here each censored observations could have its own l

Re: [R] rgl snapshot on headless server

2013-09-17 Thread Andreas Maunz
server that would do it (tried XF4VNC and TigerVNC). Thanks for your help anyway! Andreas On Sun, Sep 15, 2013 at 1:18 PM, Duncan Murdoch wrote: > On 13-09-15 6:02 AM, Andreas Maunz wrote: > >> The other write* options seem to limited to me. I am looking for a way >> to obtain

Re: [R] rgl snapshot on headless server

2013-09-15 Thread Andreas Maunz
extures/sunsleep.png", package="rgl"), back="filled" ) foo<-readline('Enter to cont') I can skip through the steps and check at each with xwud, but it just manages to draw the white background. The two last steps fail completely, I only receive a black screen.

Re: [R] rgl snapshot on headless server

2013-09-11 Thread Andreas Maunz
file format? Thanks Andreas On Tue, Sep 10, 2013 at 6:38 PM, Duncan Murdoch wrote: > On 10/09/2013 10:58 AM, Andreas Maunz wrote: > >> Hi all, >> >> I have a shiny app, in which I want to use rgl's snapshot function. I am >> running Xvfb on my server so that

[R] rgl snapshot on headless server

2013-09-10 Thread Andreas Maunz
GL and send the results to the browser. But the rgl.snapshot or rgl.postscript functionality do not work, i.e. they produce black or empty images. I assume this is due to Xvfb. Any chance I can create snapshots? Thanks Andreas [[alternative HTML version deleted]] __

Re: [R] Problems installing FAME package on Windows 7.

2013-08-23 Thread Andreas Dibiasi
un the function separatly and specify lib.loc as .libsPaths() the function seems to work. My .libPaths() are the following: [1] "C:/APPS/R/R-3.0.1/library" "C:/R/R-3.0.1/library" could this be an issue? which other environmental variable could be responsible? Best, Andreas Dibiasi O

[R] Problems installing FAME package on Windows 7.

2013-08-20 Thread Andreas Dibiasi
I have problems to install the FAME package. I have installed the fame client 9.3.32.2 and the environmental variable is set. However, for some reason, the warnings "error: argument "lib.loc" > is missing, with no default" pops up under the i386 installation and "LoadLibrary failure: %1 is not a v

Re: [R] suppress startup messages from default packages

2013-07-15 Thread Andreas Leha
Hi Martin, Martin Morgan writes: > On 07/15/2013 06:25 AM, Duncan Murdoch wrote: >> On 15/07/2013 8:49 AM, Andreas Leha wrote: >>> Hi Helios, >>> >>> "Helios de Rosario" writes: >>> >>> >> Hi all, >>> >>

Re: [R] suppress startup messages from default packages

2013-07-15 Thread Andreas Leha
Hi David, David Winsemius writes: > On Jul 15, 2013, at 4:11 AM, Andreas Leha wrote: > >> Hi all, >> >> several packages print messages during loading. How do I avoid to see >> them when the packages are in the defaultPackages? > > I'm pretty sure

Re: [R] suppress startup messages from default packages

2013-07-15 Thread Andreas Leha
all messages. There are hacks to do that even for the very persistent messages [fn:1]. My question is how to suppress these messages, when it is not me who calls library() or require(), but when the package is loaded during R's startup through the defaultPackages option. Regards, Andreas F

[R] suppress startup messages from default packages

2013-07-15 Thread Andreas Leha
lehash")) | rm(old) ` I get as last line when starting R: , | filehash: Simple key-value database (2.2-1 2012-03-12) ` Another package with (even more) prints during startup is tikzDevice. How can I avoid to get these messages? Regards, Andreas

[R] problems with setClass or/and setMethod

2013-06-10 Thread andreas betz
= "Trajectories", y = "Trajectories", ... = "Trajectories") omitted in the method definition cannot be in the signature Did anything change in the transition to R-3.0? is there any other, more recent introduction to S4 classes recommended? Thank you for your help. A

Re: [R] R web application development

2013-03-19 Thread Andreas Recktenwald
posting guide http://www.R-project.org/posting-guide.html > and provide commented, minimal, self-contained, reproducible code. Hi, have you tried the 'shiny' and 'shiny-server' packages by Rstudio? http://www.rstudio.com/shiny/ -- Dipl.-Kfm. Andreas Recktenwa

[R] Function ar() in package stats: AIC procedure and time series length

2013-03-11 Thread Andreas Klein
, -0.00688,  0.00022,  0.00472, -0.00028,  0.00362,  0.00122,  0.00852, -0.00118) y <- c(-0.0046,  0.0018,  0.0088,  0.0051, -0.0020,  0.0018,  0.0083, -0.0012,  0.0096, -0.0014, -0.0035, -0.0107, -0.0052, -0.0027, -0.0095,  0.0008,  0.0070, -0.0018,  0.0031, -0.0019,  0.0097, -0.0022) #Calculations

[R] Ordered Probit/Logit with random coefficients

2013-01-21 Thread Andreas Karpf
Hello, I searched everywhere but I didn't find what I want, that is why I as the question here. Threads discussing this issue on this mailing list are already quite old. Does anybody know of a function in R which allows to estimate ordered probit/logit model with random coefficients. The only mix

[R] Parameter estimates for each observation (ordered choice)

2013-01-09 Thread Andreas Karpf
I have several demographic variables with which I want to explain the ordered choice of individuals within a survey in an ordered choice (probit or logit, this is not important) framework. Standard ordered choice estimations of course just give me aggregate/average parameter estimates. For my task

Re: [R] sqldf Date problem

2012-11-04 Thread Andreas Recktenwald
-17 6 1/18/2011 2011-01-18 7 1/19/2011 2011-01-19 8 1/20/2011 2011-01-20 9 1/21/2011 2011-01-21 10 1/22/2011 2011-01-22 11 1/23/2011 2011-01-23 12 1/24/2011 2011-01-24 13 1/25/2011 2011-01-25 On Sat, Nov 3, 2012 at 4:22 PM, Andreas Recktenwald wrote: Dear R-help readers, i've created

[R] sqldf Date problem

2012-11-03 Thread Andreas Recktenwald
Dear R-help readers, i've created a database for quotes data (for 4 years; 2007 -- 2010) with the sqldf package. This database contains a column "Date" in the format mm/dd/. The table in the database is called "main.data" and the database itself "Honda". I tried to get the Data just f

[R] sqldf Date problem

2012-11-03 Thread Andreas Recktenwald
Dear R-help readers, i've created a database for quotes data (for 4 years; 2007 -- 2010) with the sqldf package. This database contains a column "Date" in the format mm/dd/. The table in the database is called "main.data" and the database itself "Honda". I tried to get the Data just

[R] segfault in gplots::heatmap.2

2012-08-30 Thread Andreas Leha
egfault offending <- dget("offending.txt") heatmap.2(x=offending) #+end_src Interestingly, I get a segfault when loading any of c("XLConnect", "venneuler", "xlsx"), which all depend on rJava. But loading rJava on its own did not produce a segfault.

[R] Generate distribution from asymmetrical confidence interval and mean

2012-03-28 Thread Andreas Emanuelsson
more and less right skewed in each case, and I need the distributions to simulate joint uncertainty of a model using these inputs (Monte Carlo). Best regards Andreas [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

Re: [R] julian() and numerical noise

2012-03-24 Thread Andreas Eckner
Very helpful, thanks! Andreas On 3/24/2012 10:57 AM, R. Michael Weylandt wrote: It might be easier to just work with a character representation rather than fighting the floating point demons: fnc<- function(date){ as.Date(paste(date, "01", sep = "."), format = &

Re: [R] julian() and numerical noise

2012-03-24 Thread Andreas Eckner
converting a dates in .MM format to a CHRON object. Here is an example for future reference (the output might be highly system-dependent): > library(chron) > date <- 2012.02 > y <- floor(date) > m <- (date %% 1) * 100 > julian(2, 1, 2012) - julian(m, 1

Re: [R] julian() and numerical noise

2012-03-23 Thread Andreas Eckner
Sorry, forgot to mention: R version 2.14.2 The function julian() is part of the "chron" package of the base distribution. On 3/23/2012 2:47 PM, Uwe Ligges wrote: On 23.03.2012 15:45, Andreas Eckner wrote: Hi, does anybody know if the following behavior of julian() is intentional?

[R] julian() and numerical noise

2012-03-23 Thread Andreas Eckner
ubject to numerical noise in the 'day' argument, but not in the 'month' and 'year' argument? Another example: > julian(2 - 1e-15, 1, 2012) - julian(1 - 1e-15, 1, 2012) [1] 30 Thanks, Andreas __ R-help@r-project.org

[R] Multiple data frames in single csv output

2012-03-22 Thread Andreas Emanuelsson
e any other neat way to join data.frames but ignoring all the indexing? (maybe turning them into text and printing them as a csv file)? Best regards Andreas Emanuelsson [[alternative HTML version deleted]] __ R-help@r-project.org mailing list

Re: [R] 3-parametric Weibull regression

2012-01-31 Thread Pfrengle Andreas (GS-SI/ENX7)
eibull exponent ($scale) isn't modelled in dependence of the factors, but constant for the whole dataset? a <- 1 / fit1$scale b <- function(Gruppe, d) exp(fit1$coeff[1] + fit1$coeff[2] * Gruppe + fit1$coeff[3] * d + fit1$coeff[4] * Gruppe * d) x <- seq(0,5,.02) plot(x, dweibull(x - o

[R] 3-parametric Weibull regression

2012-01-26 Thread Pfrengle Andreas (GS-SI/ENX7)
t a failure by survregDtest(survreg.distributions$weibull3p) that x0 is missing (since the transformation function is expected to take only one argument). Any ideas? Or is there maybe a package somewhere that supports 3-parametric weibull regressio

Re: [R] R on Android

2012-01-02 Thread Andreas Borg
Ken Hutchison schrieb: If stranded on a desert island and getting R fully functional on Droid would save you, I'd advise looking for food; else be prepared to name your phone Wilson. Ken Hutchison I propose to add that to fortunes :-D Andreas Borg -- Andreas Borg Abte

Re: [R] What's the baseline model when using coxph with factor variables?

2011-12-02 Thread Andreas Schlicker
half Of David Winsemius Sent: Thursday, December 01, 2011 9:36 AM To: a.schlic...@nki.nl Cc: r-help@r-project.org Subject: Re: [R] What's the baseline model when using coxph with factor variables? On Dec 1, 2011, at 12:00 PM, Andreas Schlicker wrote: Hi all, I'm trying to fit a Cox regressi

Re: [R] help in dbWriteTable

2011-12-02 Thread Andreas Borg
, Andreas arunkumar schrieb: hi I need some help in dbWriteTable. I'm not able to insert the rows in the table if the column order are not same in the database and in the dataframe which i'm inserting. Also facing issue if the table is already created externally and inserting it th

[R] What's the baseline model when using coxph with factor variables?

2011-12-01 Thread Andreas Schlicker
Hi all, I'm trying to fit a Cox regression model with two factor variables but have some problems with the interpretation of the results. Considering the following model, where var1 and var2 can assume value 0 and 1: coxph(Surv(time, cens) ~ factor(var1) * factor(var2), data=temp) What is t

[R] Question on density values obtained from kde2d() from package MASS

2011-11-24 Thread Andreas Klein
Hello, I am a little bit confused regarding the density values obtained from the function kde2d() from the package MASS because the are not in the intervall [0,1] as I would expect them to be. Here is an example: x <- c(0.0036,0.0088,0.0042,0.0022,-0.0013,0.0007,0.0028,-0.0028,0.0019,0.0026,-0

[R] Block length for Bivariate Stationary Bootstrap for Inference for Correlation

2011-11-18 Thread Andreas Klein
Hello, I am searching for a method, algorithm or some solution to the following problem: I need an estimator for the block length of two series to bootstrap the correlation coefficient with the stationary bootstrap procedure. The problem here: the two series have different block lengths! Here

[R] Function not found, maybe respective package has to be put in environment?

2011-11-13 Thread Andreas Klein
Hello everybody, I have a problem and would like to start with an example: library(snow) library(tseries) fn <- function(x) adf.test(x) clusterApply(cl=cl, x=x , fun=fn) R cannot find the function adf.test() because it is inside the function fn(). This problem does not occur when, for example

[R] Repeated Measures Regression

2011-10-07 Thread Andreas Pedroni
this be suitable? Is it a problem, if the Variable Risk is autocorrelated within Subjects? Thank you very much for any help and sorry if this is a too trivial question for the list. Kind regards Andreas Pedroni __ R-help@r-project.org mailing list

[R] Robust covariance matrix with NeweyWest()

2011-09-28 Thread Andreas Klein
Dear R-users, I would like to compute a robust covariance matrix of two series of realizations of random variables: ###Begin Example### data <- cbind(rnorm(100), rnorm(100)) model <- lm(data ~ 1) vcov(model) library(sandwich) NeweyWest(model) #produces an error ###End Example### NeweyWest()

[R] Gui editor / viewer for large data

2011-07-14 Thread Andreas Borg
element to view and potentially edit data - which only loads as much of the data as it displays at one point in time, - thereby able to handle datasets that are too large to fit into main memory Thanks for any suggestions, Andreas -- Andreas Borg Medizinische Informatik UNIVERSITÄTSMEDIZIN de

Re: [R] How to remove last 3 letters

2011-07-13 Thread Andreas Borg
estions are highly appreciated __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. -- Andrea

[R] Testing the proportional odds assumption of an ordinal generalized estimating equations (GEE) regression model

2011-06-28 Thread Andreas Karlsson
Dear list members, I am estimating an ordinal generalized estimating equations (GEE) regression model on repeated measurements data. Is there any way (preferably in R) to test the proportional odds assumption of this model? Thanks in advance. Andreas

[R] problems with plots in loop (corrected Email)

2011-06-14 Thread Andreas Betz
smooth.spline(Zeit, guess) der <- 100*(predict(fit, Zeit, deriv = 1))$y The problem appears to be hidden in the mfrow statement. I spent quite a bit of time on this problem, thus any help or new ideas would be very much appreciated. Thank you Andreas Betz Scientist andreasbetz@ Dea

[R] problems with plots in loop

2011-06-14 Thread Andreas Betz
ooth.spline(Zeit, guess) der <- 100*(predict(fit, Zeit, deriv = 1))$y The problem appears to be hidden in the mfrow statement. I spent quite a bit of time on this problem, thus any help or new ideas would be very much appreciated. Thank you Andreas Betz Scientist andreasb..

Re: [R] Compiling C-code in Windows

2011-05-31 Thread Andreas Borg
dll C:\trial.c' under Linux, and expecting working DLL? Of course not, because Windows and Linux libraries differ by more than the filename extension. Once everything is set up, compiling in windows is pretty easy thanks to the toolset. Good luck, Andreas -- Andreas Borg Medizinische

Re: [R] Options for print()

2011-05-05 Thread Andreas Borg
Also take a look at sprintf, which offers everything the C-equivalent has. sprintf returns a string which can be sent to the console via cat. Andreas baptiste auguie schrieb: Hi, Try this, cat(format("The TITLE", width=80, justify="centre")) HTH, baptiste On 5 May 2

Re: [R] user input

2011-04-20 Thread Andreas Borg
Hi Ivan, there are also data.entry() and edit(), which are more tailored to data frames and might look awkward for entering a single value. You could also take a look at CRAN, maybe there is a GUI package that does what you want. Andreas Ivan Calandra schrieb: Dear users, I have looked

[R] Sweave stops when opening X11 device fails

2011-04-20 Thread Andreas Borg
I think it would be nice if the script just kept on running without plotting on the screen if opening X11 fails, just like it does when no X11 is available at all. Best regards, Andreas -- Andreas Borg Medizinische Informatik UNIVERSITÄTSMEDIZIN der Johannes Gutenberg-Universität Institut f

Re: [R] Where did my packages go ?

2011-04-19 Thread Andreas Borg
installation. Andreas eric schrieb: Running linux 10.04 ubuntu. Looks like Ubuntu automatically updated to version 2.13. I was running version 2.12 until today. But now I'm getting error messages when I use the require or library command and one of the packages I've downloaded in the past. F

Re: [R] mapply to lapply

2011-04-18 Thread Andreas Borg
:nrow(cells), function(rowInd, i, j), [your function goes here], i=i, j=j) Viele Grüße, Andreas Alaios schrieb: Dear Andreas, I would like to thank you for your reply. I have tried two alternatives but none of the two worked out: F2[i+1,j+1]<-sum(lapply(1:nrow(cells), function(rowInd) Fwit

Re: [R] mapply to lapply

2011-04-18 Thread Andreas Borg
d,2],b=cells[rowInd,4],c=cells[rowInd,1],d=cells[rowInd,3],e=cells[rowInd,5])) Best regards, Andreas Alaios schrieb: Dear all, I would like to ask your help concerning converting a mapply function to lapply. The reason is that I would like to use mclapply which requires lapply syntax.

Re: [R] How do I make this faster?

2011-04-11 Thread Andreas Borg
Hi Hasan, I'd be happy to help you, but I am not able to run your code. You use commandArgs to retrieve arguments of the R program, but which ones do you actually provide? Best regards, Andreas Hasan Diwan schrieb: I was on vacation the last week and wrote some code to run a 50

Re: [R] random sampling with levels and with replacement

2011-04-08 Thread Andreas Borg
which(df$V2=="good") isBad <- which(df$V2=="bad") sampleGood <- df[sample(isGood, replace=TRUE),] sampleBad <- df[sample(isBad, replace=TRUE),] summary(rbind(sampleGood, sampleBad)) Please include a more specific example with test data (for "final" in this case)

Re: [R] Assign Names of columns in data.frame dinamically

2011-04-01 Thread Andreas Borg
Hi Marcos, I'd be glad to help, but your example is not usable for anyone but yourself - I don't have the files that your code reads and I don't have any idea what the resulting data frames look like. Please provide examples of the data and the result. Andreas Marcos

Re: [R] a for loop to lapply

2011-03-30 Thread Andreas Borg
ease make sure your examples are self-contained and executable (i.e. they don't throw an error and they do not depend on undefined variables, such as dimx etc. in this case) Best regards, Andreas -- Andreas Borg Medizinische Informatik UNIVERSITÄTSMEDIZIN der Johannes Gutenberg-Universität I

[R] (no subject)

2011-03-14 Thread Andreas Emanuelsson
ass[j] : NA/NaN argument Thank you very much for your time if you can help me! /Andreas [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://ww

[R] (no subject)

2011-03-11 Thread Andreas Emanuelsson
e column? Any quick suggestions? Andreas [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide comm

[R] panel data

2011-02-21 Thread Andreas Sveding Otterström
Hello I have a question about how to regress panel data in R. What are the appropriate commands and preparations that need to be made when regressing panel data? /Andreas Otterstrom [[alternative HTML version deleted

Re: [R] RSQLite to input dataframe

2011-01-04 Thread Andreas Borg
with dbWriteTable. There might be a way to read files directly from SQLite, but I don't know about that. Best regards, Andreas -- Andreas Borg Medizinische Informatik UNIVERSITÄTSMEDIZIN der Johannes Gutenberg-Universität Institut für Medizinische Biometrie, Epidemiologie und Inf

[R] survival package - calculating probability to survive a given time

2010-12-10 Thread Andreas Wittmann
lity to get the probability for surviving 20 units of time ind <- s$time >= 20 sum(ind) / length(ind) Thanks and best regards Andreas __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide htt

Re: [R] How to find out if a data frame has automatic row names?

2010-12-07 Thread Andreas Borg
That does it. Thanks! Check ?.row_names_info Peter Ehlers Andreas -- Andreas Borg Medizinische Informatik UNIVERSITÄTSMEDIZIN der Johannes Gutenberg-Universität Institut für Medizinische Biometrie, Epidemiologie und Informatik Obere Zahlbacher Straße 69, 55131 Mainz www.imbei.uni

[R] How to find out if a data frame has automatic row names?

2010-12-07 Thread Andreas Borg
mes(df2) [1] "1" "2" > # when converting to a matrix, it makes a difference > as.matrix(df) a b [1,] 1 3 [2,] 2 4 > as.matrix(df2) a b 1 1 3 2 2 4 Thanks for any suggestion, Andreas -- Andreas Borg Medizinische Informatik UNIVERSITÄTSMEDIZIN der Johannes Gutenberg-

[R] Quadratic programming with semi-definite matrix

2010-12-04 Thread Andreas Jensen
ction is not positive definite! Thanks in advance, Andreas Jensen __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.

Re: [R] Performance tuning tips when working with wide datasets

2010-11-29 Thread Andreas Borg
inal data. I am not sure if this works with other database engines (there are a couple of database interfaces on CRAN); when I tried to use RPostgreSQL, it created the whole result set at once, leading to the same memory problem. Maybe that behavior can be changed by some config variable. Best re

[R] read.table / type.convert: Unexpected column class

2010-11-23 Thread Andreas Borg
x27; is treated as "other decimal or hexadecimal constant" while the leading whitespace before "2" is stripped. Actually, leaving just the trailing space after '3' without the comment gives the same result. Is there any reason for this behaviour? Best regards,

[R] Closing unreferenced result sets in dbi / RSQLite

2010-11-04 Thread Andreas Borg
hich is especially annoying when working with a temporary data base where everything is discarded on disconnection. Is there any way to create a new reference to the pending result set or to close result sets which are not bound to a variable? Thanks for any suggestion, Andreas -- Andreas

[R] gamma glm - using of weights gives error

2010-10-24 Thread Andreas Wittmann
ence in the results? What can i do to locate the root of the error and how can i avoid this error? Can i to the regression without the weights argument? Thanks and best regards Andreas __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/

[R] RWeka - Error in model.frame.default - evaluate_Weka_classifier

2010-10-19 Thread Andreas Jansson
fier function? Finally, some info on my system: R version 2.10.1 (2009-12-14) RWeka version [Most recent as of 2010-10-19, can't find the exact number anywhere, sorry!] Java version 1.6.0_22 (from the sun-java6-jdk package) OS: Ubuntu 10.04.1 Many thanks, Andreas Jansson __

  1   2   3   >