Re: [R] [Fwd: Re: How to apply functions over rows of multiple matrices]

2007-08-10 Thread Johannes Hüsing
> 1. matrices are stored columnwise so R is better at column-wise operations > than row-wise. I am seeing this by my code which contains more t() than what seems healthy. However, the summaries are patient-wise over repeated measurements. Out of convention, I am storing patients in rows and measur

[R] [Fwd: Re: How to apply functions over rows of multiple matrices]

2007-08-10 Thread Johannes Hüsing
[Apologies to Gabor, who I sent a personal copy of the reply erroneously instead of posting to List directly] [...] > Perhaps what you really intend is to > take the average over those elements in each row of the first matrix which correspond to 1's in the second in the corresponding > row of the

[R] How to apply functions over rows of multiple matrices

2007-08-09 Thread Johannes Hüsing
Dear ExpRts, I would like to perform a function with two arguments over the rows of two matrices. There are a couple of *applys (including mApply in Hmisc) but I haven't found out how to do it straightforward. Applying to row indices works, but looks like a poor hack to me: sens <- function(test,

Re: [R] write to clipboard under Linux

2006-11-01 Thread Johannes Hüsing
Prof Brian Ripley: > Linux does not have a clipboard but an X11 session has primary and > secondary selections. ?file says > [...] > so RTFM applied. Indeed it did here. Many thanks for the externsive answer despite the lucid help entry already being there! _

Re: [R] latex() in Hmisc and n.group

2006-10-29 Thread Johannes Hüsing
> This is a bug that has now been fixed. Until a new release of Hmisc > appears see the following to get a corrected version of latex( ): > http://biostat.mc.vanderbilt.edu/twiki/bin/view/Main/LatestRFunctions - > you will need getLatestSource('latex') It worked. Many thanks, and apologies for po

[R] write to clipboard under Linux

2006-10-28 Thread Johannes Hüsing
ein My interpretation is that I am not allowed to write into the clipboard from a program called from R. Is there a way to change this behaviour? -- Johannes HüsingThere is something fascinating about science. [EMAIL PROTECTED] One gets such wholesale returns of conjecture

[R] (kein Betreff)

2006-10-28 Thread Johannes Hüsing
Frank Harrell rote in a message dating from Oct 8th: > n.group is an argument to latex.default in the Hmisc package I must admit that I can't find it in the function head, which reads on my installation: function (object, title = first.word(deparse(substitute(object))), file = paste(title, "

Re: [R] acos(0.5) == pi/3 FALSE

2006-09-19 Thread Johannes Hüsing
Peter Dalgaard: > Ben Bolker <[EMAIL PROTECTED]> writes: >> 1. compose your response > I've always wondered why step 1. - often the time-consuming bit - is not > listed last. The advice applies to the situation when answering immediately would be your knee-jerk reaction. It is assumed tha

Re: [R] Dividing objects in classes using function sample()

2006-09-03 Thread Johannes Hüsing
> I want to divide n objects in k classes and need an output with all > (n+1)(n+2)/2 possibilities. That's the "set of compositions". You may use the partitions package and proceed from there (provided the brute-force method suggested by Gabor is not viable').

Re: [R] Dividing objects in classes using function sample()

2006-09-02 Thread Johannes Hüsing
> I want to divide n objects in k classes and need an output with all > (n+1)(n+2)/2 possibilities. That's the "set of compositions". You may use the partitions package and proceed from there (provided the brute-force method suggested by Gabor is not viable').

Re: [R] split a y-axis to show data on different scales

2006-08-15 Thread Johannes Hüsing
> The pro's and con's of using "scale breaks" were discussed by > Cleveland (1985) The Elements of Graphing Data (Wadsworth, pp. 85-91, > 149). I don't know what Cleveland said about this is the second edition Spencer Graves: > but I believe there are times when scale breaks are > appropr

Re: [R] S curve via R

2006-07-25 Thread Johannes Hüsing
> Hello sir: > How can I get "S curve" function via R? > For SPSS,the function is:y=exp(b0+b1/x) > I am not sure if this is the answer you want, but Scurve <- function(x, b0=0, b1=1) { exp(b0+b1/x) } should do what you request. Greetings Johannes _

Re: [R] change the name of file

2006-07-24 Thread Johannes Hüsing
> Of course file name Data_i.txt will be the same for changing i, > unfortunately. ?paste, e.g. paste("Data_", i, ".txt", sep="") __ R-help@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http

Re: [R] dotchart with log scale?

2006-07-23 Thread Johannes Hüsing
> Dear all, > I would like to draw a dot chart on a log scale. > What is the syntax for this? A barchart may use > log="x", but trying this with dotchart() leads > to an error message. ok, I extended dotchart() with a log option. The main reason I am sticking to dotchart() is adjusting of axes wit

[R] dotchart with log scale?

2006-07-20 Thread Johannes Hüsing
Dear all, I would like to draw a dot chart on a log scale. What is the syntax for this? A barchart may use log="x", but trying this with dotchart() leads to an error message. Greetings Johannes __ R-help@stat.math.ethz.ch mailing list https://stat.eth

Re: [R] Manipulating code?

2006-05-23 Thread Johannes Hüsing
> Johannes Hüsing <[EMAIL PROTECTED]> writes: [...] > I think a simple example of what you are trying to do might be needed. > I don't think so, as ... > But take a look at the help pages for assign() and get(). ... this seems to be what I was looking for. Many thanks

[R] Manipulating code?

2006-05-23 Thread Johannes Hüsing
Dear expeRts, I am currently struggling with the problem of finding cut points for a set of stimulus variables. I would like to obtain cut points iteratively for each variable by re-applying a dichotomised variable in the model and then recalculate it. I planned to have fixed names for the dichotom

Re: [R] fixed variance of lmer object

2006-03-03 Thread Johannes Hüsing
duction to R". Am I missing something here? If desired, I can write a patch for R-intro.texi. Greetings Johannes -- Johannes HüsingThere is something fascinating about science. [EMAIL PROTECTED] One gets such wholesale returns of conjecture

[R] fixed variance of lmer object

2006-03-01 Thread Johannes Hüsing
Dear expRts, could anybody nudge me toward some documentation about how to extract the variance-covariance matrix of the fixed parameters from an lmer object? Best Johannes -- Johannes HüsingThere is something fascinating about science. [EMAIL PROTECTED] One gets such wholesale

Re: [R] a series of 1's and -1's

2006-01-12 Thread Johannes Hüsing
> You could try to zip your data file and see whether there is a change in > file size that you feel is significant in which case the series is not > random (-: ... after converting the -1s and 1s to bits, of course. __ R-help@stat.math.ethz.ch mailing

Re: [R] reading contigency tables

2006-01-10 Thread Johannes Hüsing
Just a remark: > I need some help using read.ftable to read a contingency table. My columns > are organized as follows: > order--family--species--location--number of individuals As the first three variables are nested, I would expect this table to contain a lot of structural zeroes. I understand

[R] graphics: axis label

2006-01-10 Thread Johannes Hüsing
Hello, par(las=1) sets the orientation of the axis labels to horizontal. That is, the tick mark labels. How do I set the orientation of the axis label, which annotates the variable plotted along the axis, to horizontal? Sorry for asking such a basic question here, but I haven't found anything in t

RE: [R] normality test

2005-04-28 Thread Johannes Hüsing
Bert wrote: >>You probably have noticed that >> I'm quite new >> to statistics, but I'm working on that... >> >> > And you want to use Bayesian methods?! > I was always under the impression that it's mostly a matter of mindset if you go Bayesian or frequentist, not of your statistical skills. [..

Re: [R] error messages on R CMD check

2005-04-06 Thread Johannes Hüsing
Hallo Uwe, Uwe Ligges schrieb: [...] > In Johannes' case, the problem is different, because the error message > is not that clear. > Johannes, can you install and load the package? Is the DESCRIPTION file > correct? If so, you might want to send the package in a private message... I have found the

[R] error messages on R CMD check

2005-03-30 Thread Johannes Hüsing
Dear all, I am trying to wrap up a package. On entering R CMD check, I get the following error messages: [...] * checking S3 generic/method consistency ... WARNING Error in .try_quietly({ : Error in library(package, lib.loc = lib.loc, character.only = TRUE, verbose = FALSE) : package/names

Re: [R] R equivalent to funcall?

2005-03-17 Thread Johannes Hüsing
> On Thu, 17 Mar 2005 13:20:32 +0100 (CET), Johannes Hüsing > <[EMAIL PROTECTED]> wrote : [...] >>Is there a way that the dot-dot-dot argument of a >>function accepts a list as single arguments, such >>as funcall in several Lisp dialects? > > do.call() comes c

[R] R equivalent to funcall?

2005-03-17 Thread Johannes Hüsing
Dear all, I have a list of time series and want to plot them. Is there a way that the dot-dot-dot argument of a function accepts a list as single arguments, such as funcall in several Lisp dialects? Greetings Johannes __ R-help@stat.math.ethz.ch maili