Displaying 20 results from an estimated 20000 matches similar to: "? extended rep()"
2009 Aug 12
3
Random sampling while keeping distribution of nearest neighbor distances constant.
Dear All,
I cannot find a solution to the following problem although I imagine
that it is a classic, hence my email.
I have a vector V of X values comprised between 1 and N.
I would like to get random samples of X values also comprised between
1 and N, but the important point is:
* I would like to keep the same distribution of distances between the X values *
For example let's say N=10 and
2020 Aug 04
2
Problem with intermediate certificate (tls cafile)
I have several samba servers on Debian 10 all using :
samba 2:4.9.5+dfsg-5+deb10u1 amd64
I use tls cafile, tls certfile and tls keyfile with certificates from
Sectigo (https://cert-manager.com)
And when checking my connexion from the samba server, or from outside,
I've got "unable to verify the first certificate" even if tls_cafile is
provided in smb.conf.
What is wrong
2009 Feb 05
3
"open-ended" plot limits?
Hi Folks,
Maybe I've missed it already being available somehow,
but if the following isn't available I'd like to suggest it.
If you're happy to let plot() choose its own limits,
then of course plot(x,y) will do it.
If you know what limits you want, then
plot(x,y,xlim=c(x0,x1),ylim(y0,y1)
will do it.
But sometimes one would like to
a) make sure that (e.g.) the y-axis has a
2008 Mar 12
1
[follow-up] "Longitudinal" with binary covariates and outcome
Hi again!
Following up my previous posting below (to which no response
as yet), I have located a report which situates this type
of question in a longitudinal modelling context.
http://www4.stat.ncsu.edu/~dzhang2/paper/glm.ps
Generalized Linear Models with Longitudinal Covariates
Daowen Zhang & Xihong Lin
(This work seems to originally date from around 1999).
They consider an outcome Y,
2010 May 12
6
A primitive OO in R -- where next?
Greetings All,
Out of curiosity, I've just done a very primitive experiment:
Obj <- list(Fun=sum, Dat=c(1,2,3,4))
Obj$Fun(Obj$Dat)
# [1] 10
That sort of thing (much more sophisticated) must be documented
mind-blowingly somewhere. Where?
Where I stand right now: The above (and its immediately obvious
generalisations, like Obj$Fun<-cos) is all I know about it so far.
Ted.
2009 May 17
2
Output of binary representation
I am interested in studying the binary representation of numerics
(doubles) in R, so am looking for possibilities of output of the
internal binary representations. sprintf() with format "a" or "A"
is halfway there:
sprintf("%A",pi)
# [1] "0X1.921FB54442D18P+1"
but it is in hex.
The following illustrate the sort of thing I want:
1.1001 0010 0001 1111
2010 Sep 08
6
'par mfrow' and not filling horizontally
Greetings, Folks.
I'd appreciate being shown the way out of this one!
I've been round the documentation in ever-drecreasing
circles, and along other paths, without stumbling on
the answer.
The background to the question can be exemplified by
the example (no graphics window open to start with):
set.seed(54321)
X0 <- rnorm(50) ; Y0 <- rnorm(50)
2010 Jun 18
4
Drawing sample from a circle
Hi, I would like to draw 10 uniformly distributed sample points from a circle with redius one and centered at (0,0). Is there any R function to do that?
Thanks,
[[alternative HTML version deleted]]
2008 Apr 27
2
Deb-4.0 Etch and sources.list for R
Hi Folks,
I'm running Debian-4.0 Etch, installed last September
from a DVD, and regularly updated as things arise.
I have R version 2.4.0 Patched (2006-11-25 r39997)
installed (initially at the time of first installation
of Debian, as provided by Debian), along with a variety
of packages.
I'd like to be able to connect to the CRAN repositories
for Debian R, for updates etc.
When I visit
2008 Oct 11
2
R vs SPSS contrasts
Hi Folks,
I'm comparing some output from R with output from SPSS.
The coefficients of the independent variables (which are
all factors, each at 2 levels) are identical.
However, R's Intercept (using default contr.treatment)
differs from SPSS's 'constant'. It seems that the contrasts
were set in SPSS using
/CONTRAST (varname)=Simple(1)
I can get R's Intercept to match
2010 Jul 14
3
Convergent series
What are some reliable R functions that can compute the value of a
convergent series?
David
--
David R. Bickel, PhD
Associate Professor
Ottawa Institute of Systems Biology
Biochem., Micro. and I. Department
Mathematics and Statistics Department
University of Ottawa
451 Smyth Road
Ottawa, Ontario K1H 8M5
http://www.statomics.com
Office Tel: (613) 562-5800 ext. 8670
Office Fax: (613) 562-5185
2010 Jun 15
2
Unspecified [upper] xlim/ylim?
Greetings!
I would like to be able to specify a fixed (say) lower limit
for plotting, while leaving the upper limit "floating, when
plotting. The context is that the maximum in the data to be
plotted is unpredictable, being the consequence of a simulation,
whereas I know that it cannot be less than (say) 0; and I want
to fix the lower limit at 0 in any plot, leaving the upper limit
to be
2007 Aug 15
3
Covariance of data with missing values.
I have a data matrix X (n x k, say) each row of which constitutes an
observation of
a k-dimensional random variable which I am willing, if not happy, to
assume to be
Gaussian, with mean ``mu'' and covariance matrix ``Sigma''. Distinct
rows of X may
be assumed to correspond to independent realizations of this random
variable.
Most rows of X (all but 240 out of 6000+ rows)
2007 Nov 30
2
Organising tick-marks in plot()
Hi Folks,
I'm advising someone who's a beginner with R,
and therefore wants the simplest answer possible.
The issue is to produce a plot using
plot(x,y,...)
where, on the X-axis, the tick-marks should be
on the lines of:
-- Range of X-axis: 0:1000
-- tick-marks labelled "0","200",...,"800","1000"
-- unlabelled tick-marks every 50 from 0 to
2008 Oct 12
1
png(): Linux vs Windows
Hi Folks,
Quick question. I have the following line in an R code file
which runs fine on Linux:
if(PNG) png(GraphName,width=12,height=15,units="cm",res=200)
I learn that, when the same code was run on a Windows machine,
there was the following error:
Error in png(GraphName,width=12,height=15,units="cm",res=200):
unused argument(s) (units = "cm")
Sorry to
2009 Mar 31
2
"digits" in round()
Hi Folks,
Compare
print(1234567890,digits=4)
# [1] 1.235e+09
print(1234567890,digits=5)
# [1] 1234567890
Granted that
digits: a non-null value for 'digits' specifies the minimum
number of significant digits to be printed in values.
how does R decide to switch from the "1.235e+09" (rounded to
4 digits, i.e. the minumum, in "e" notation) to
2007 Dec 09
2
Adding info from summary(lm(...)) to plot
Hi Folks,
Say I have 2 continuous variables X,Y.
I can of course plot (X,Y) with
plot(X,Y,pch="+",col="blue")
say, and add the regression line from lm(Y~X)
by extracting the coefficients 'a' of Intercept
and 'b' of X from Y.lm <- lm(Y~X).
Now, however, I want to have not only a general
explanatory title such as
main="Plot of Y against X"
2007 Oct 19
6
r achives
sorry but how do i accsess r archives
__________________________________________________
[[alternative HTML version deleted]]
2009 Sep 13
2
ountour "resolution"
Greetings all!
I'm wanting to plot contours of a function, and I don't
want to evaluate it at a dense grid of points (i.e. I don't
want a huge array of values).
Say I have a vector of x-values such as x <- 0.1*c(1:10),
and the same for y <- 0.1*(0:10).
I then evaluate a 10*10 matrix z of values of z = f(x,y).
so I can then do
CL <- contourLines(x,y,z)
and get a list of
2008 May 20
1
contr.treatments query
Hi Folks,
I'm a bit puzzled by the following (example):
N<-factor(sample(c(1,2,3),1000,replace=TRUE))
unique(N)
# [1] 3 2 1
# Levels: 1 2 3
So far so good. Now:
contrasts(N)<-contr.treatment(3, base=1, contrasts=FALSE)
contrasts(N)
# 1 2
# 1 1 0
# 2 0 1
# 3 0 0
whereas:
contr.treatment(3, base=1, contrasts=FALSE)
# 1 2 3
# 1 1 0 0
# 2 0 1 0
# 3 0 0 1
contr.treatment(3, base=1,