Hello - R-experts,
This may sound simple to many ---
We can specify the x limit by saying xlim = c(0, 20) for example.
Then the graph will show the range of x -axis between 0 - 20.
But the coordinate gap will be automatic like 0 then 5 then 15 and 20.
Is there any way by which we can set it
Hello - R-experts,
This may sound simple to many ---
We can specify the x limit by saying xlim = c(0, 20) for example.
Then the graph will show the range of x -axis between 0 - 20.
But the coordinate gap will be automatic like 0 then 5 then 15 and 20.
Is there any way by which we can set it in
I have made a few more improvements:
expand.grid.id - function(id, ...) {
vars - list(...)
nv - length(vars)
lims - sapply(vars,length)
stopifnot(length(lims) 0, id = prod(lims), length(names(vars)) == nv)
res - structure(vector(list,nv), .Names = names(vars))
if (nv 1) for(i in nv:2) {
Hi, thanks for your e-mail, but it was not the case. I found the lines
with the readline/readline.h information. This was because I did not
have the realine package installed. After installing it I was able to
install the R-2.0.0 using the .rpm binary version, however, I'm still
not able to
Hello,
I am trying to translate a matlab program into R, and I have some problems
to find out how to translate the PPVAL function.
I know that ppval(pp,xx) returns the value at the points xx of the
piecewise polynomial contained in pp( constructed by spline), but I don't
know how to do the same
try the following:
plot(0:20, axes = FALSE)
axis(2)
axis(1, at = seq(0, 20, 2), labels = seq(0, 20, 2))
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catholic University of Leuven
Address: Kapucijnenvoer 35, Leuven,
Is It (Tobit Regression) need some residual (error) assumption ?
like OLS method :
1. Normal
2. No Autocorrelation
3. Homogeneity of Variance
4. No multicollinearity
How if the assumption not fulfilled ?
[[alternative HTML version deleted]]
Yes, the R2HTML route is probably the quickest. Its just one line
of code (plus the call to load in R2HTML). Try this where iris
is a data set built into R:
library(R2HTML)
HTML( iris, file(clipboard,w), append=FALSE )
Now paste the clipboard into Excel and from there into Word.
(If you
perhaps you are looking at something like
plot(0, xlim=c(0,20), xaxt=n)
axis(1, at=pretty(0:20, 10))
Stefano
-Messaggio originale-
Da: [EMAIL PROTECTED]
[mailto:[EMAIL PROTECTED] conto di Subhabrata
Inviato: venerdì 10 febbraio 2006 7.05
A: r-help
Cc: [EMAIL PROTECTED]
Oggetto: [R] Fw:
hello,
is it possible to do something like splinefun(x,y), but with a 3rd
dim?
for example, if i have a 2-dim system like:
x-1:100
y-rexp(100,1)
func-splinefun(x,sort(y))
func(n) # returns interpolated value of y (after sorting) given x=n
# i can check this by:
plot(x,sort(y))
Hello - R-experts,
This may sound simple to many ---
We can specify the x limit by saying xlim = c(0, 20) for example.
Then the graph will show the range of x -axis between 0 - 20.
But the coordinate gap will be automatic like 0 then 5 then 15 and 20.
Is there any way by which we can set it in
Typing latent class on the R project search page gives you at least four
packeges that do latent class models mmlcr, flexmix, gllm, depmix etc
Hth, ingmar
From: Lisa Wang [EMAIL PROTECTED]
Organization: UHN\RIS
Date: Thu, 09 Feb 2006 16:31:04 -0500
To: R-Help r-help@stat.math.ethz.ch
On 2/9/06, Berton Gunter [EMAIL PROTECTED] wrote:
Folks:
R 2.2.0 on Windows.
I find the following somewhat puzzling:
a-1; x-0:1; y-x
## following works fine:
plot(x,y ,main= bquote(n[1] == .(a) ))
## following produces an error:
plot(y~x ,main= bquote(n[1] == .(a) ))
Error in
Bill,
?fft
Do you have a specific question?
On 2/9/06, Bill Hunsicker [EMAIL PROTECTED] wrote:
R-help:
I need to do a fft on a data set. I was wondering if any guidance may
be available.
Regards,
Bill
Bill Hunsicker
RF Micro Devices
7625 Thorndike Road
Greensboro, NC 27409
On Fri, 10 Feb 2006, Subhabrata wrote:
Hello - R-experts,
This may sound simple to many ---
We can specify the x limit by saying xlim = c(0, 20) for example.
Then the graph will show the range of x -axis between 0 - 20.
But the coordinate gap will be automatic like 0 then 5 then 15 and
Hi, Erin:
I tried the example in the help file for 'decompose':
m - decompose(co2)
m$figure
plot(m)
Consistent with your post, I got the following:
class(m)
[1] decomposed.ts
Then I tried the following:
methods(class=decomposed.ts)
[1]
Reshape version 0.6
===
Reshape is an R package for flexibly restructuring and aggregating
data. It's inspired by Excel's pivot tables, and it (hopefully) makes
it very easy to get your data into the shape that you want. You can find out
more at http://had.co.nz/reshape
What's
Probably ?spline and the builtin spline() function is
already sufficient for you.
If you need more, use the (standard R) 'splines' package:
library(splines)
library(help = splines)
pecqueur == pecqueur [EMAIL PROTECTED]
on Thu, 9 Feb 2006 10:53:45 +0100 (CET) writes:
pecqueur
Hi all,
I've done a simple function for import files but can't do the same think
for export files.
Could you help me solve it?
Thanks in advance for your cooperation.
wd-C:\\Teste
import-function(wd0){
wd1-paste(wd,wd0,Input\\., sep = \\)
Chris Behr wrote:
I am currently running analyses in winbugs with two different formats of
data: matrix and list. The data in the list have different dimensions than
the matrix. Do I need to create a single format for entry into R which R
passes to WinBUGS?
It is not quite clear to me what
I didn't follow this thread entirely, but I did make a LaTeX recommendation and
I know that wasn't what you were asking for. But, if I may, let me respond to
the ideas you present below in an attempt to be somewhat persuasive.
IMHO, this are horrible inefficiencies of SPSS and other packages,
Hi
plot(1:100,1:100, log=xy)
abline(v=seq(0,100,10), lty=3)
abline(h=seq(0,100,10), lty=3)
or
plot(1:100,1:100, log=xy, axes=F)
axis(1, at=(seq(0,100,20)))
axis(2, at=(seq(0,100,10)))
abline(v=seq(0,100,20), lty=3)
abline(h=seq(0,100,10), lty=3)
isn't it nice?
See some other magigraphical
Gabor Grothendieck [EMAIL PROTECTED] writes:
Yes, the R2HTML route is probably the quickest. Its just one line
of code (plus the call to load in R2HTML). Try this where iris
is a data set built into R:
library(R2HTML)
HTML( iris, file(clipboard,w), append=FALSE )
Now paste the
I have written a function called mcsamp() that is a wrapper that runs
mcmcsamp() and automatically monitors convergence and structures the
inferences into vectors and arrays as appropriate.
But I have run into a very little problem, which is that mcmcsamp()
shortens the variable names. For
Hi,
I was wondering if there were a option to have command line completion
in R like in the Bash shell?
Best regards,
David
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
Hello!
I would like to get MLE for parameter lambda of Poisson distribution. I
can use fitdistr() for this. After looking a bit into the code of this
function I can see that value for lambda and its standard error is
estimated via
estimate - mean(x)
sds - sqrt(estimate/n)
Is this MLE? With
Hello,
Is there any function / package with which you can create your own CDF
files, i.e. you create your own probe sets and based on that a new CDF
file is created?
This doesn't seem possible with makecdfenv (which needs a CDF as input)
The package makecdfenv does this:
Description: This
I ran into the same issue and again it was not having proper PATH (with include
folders set) and LD_LIBRARY_PATH set.
-Original Message-
From: André Beló [mailto:[EMAIL PROTECTED]
Sent: Thursday, February 09, 2006 11:24 PM
To: Bhola, Gautam
Cc: r-help@stat.math.ethz.ch
Subject: Re: [R]
Andrew Gelman [EMAIL PROTECTED] writes:
I have written a function called mcsamp() that is a wrapper that runs
mcmcsamp() and automatically monitors convergence and structures the
inferences into vectors and arrays as appropriate.
But I have run into a very little problem, which is that
On 10 Feb 2006 12:39:43 +0100, Peter Dalgaard [EMAIL PROTECTED] wrote:
Gabor Grothendieck [EMAIL PROTECTED] writes:
Yes, the R2HTML route is probably the quickest. Its just one line
of code (plus the call to load in R2HTML). Try this where iris
is a data set built into R:
Gregor Gorjanc [EMAIL PROTECTED] writes:
Hello!
I would like to get MLE for parameter lambda of Poisson distribution. I
can use fitdistr() for this. After looking a bit into the code of this
function I can see that value for lambda and its standard error is
estimated via
estimate -
Hi All,
I have problem to plot histogram.
What I did is:
export a .csv file from PSQL database's table.
so, inside this .csv file it looks like:
31.845
24.598
29.1223
24.715
23.1847
24.2321
25.2995
23.4261
30.7873
..
Then, I use command:
score- read.csv('file.csv', header = FALSE,sep = ,)
R-help,
I'm using the following code:
plot(1985:2005, data, type = h, lwd = 10)
which should plot histogram like figure.
The bar tops are square and I was wondering if
it is intended to be so or not (I can of course
use barplot instead).
Thanks
__
jia ding wrote:
Then, I use command:
score- read.csv('file.csv', header = FALSE,sep = ,)
hist(score, main = score)
it gives error msg:
Error in hist.default(score, main = score) :
'x' must be numeric
Can any of you know about it explain me why?
Have a look at 'score' in R
Hi,
This following function works correctly. However, I would like to improve
it by adding to the table of result the titles of columns and rows. I
prefer to define vector column by : D = (Dj = debit[j] ,
j=1:length(debit)) and vector row by: J = (Ji = time [ i ], i = 1:
length(time))
Thank you all.
par(usr)
is the perfect solution.
It *is* in the help files, but was quite hard to find.
Thomas
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
Hello all,
I am trying to cluster similar rows of a data matrix together using
hierarchical clustering and display the results in the form of a heat map.
In a Heat map, is it possible to get the order in which the rows will be
clustered , starting with the first rows joined into a group and then
Dear R helper,
I mange to transform uniform sequences to mixture
normal distributions using the following cods:
K-5
prime-c(29) , where 29 is prim number
UN-seq(1:K)%*%t(sqrt(prime))
U1-UN-as.integer(UN)
e-norMix(mu=c(-0.825,0.275), sig2 = c(0.773,0.773),
w = c(0.25,0.75), name =
use JGR,see http://stats.math.uni-augsburg.de/JGR/
2006/2/10, David Ruau [EMAIL PROTECTED]:
Hi,
I was wondering if there were a option to have command line completion
in R like in the Bash shell?
Best regards,
David
__
R-help@stat.math.ethz.ch
Colby,
Function surf() in package labdsv uses the gam() function from mgcv
to do this in conjunction with akima. You might want to look at that
routine for an idea. Currently it fits the gam as z - gam(s(x) + s(y),
but it's possible in the mgcv version of gam to fit z - gam(s(x,y)) as
Hello:
I have problems generating big JPG plots in functions. the following code runs
fine in the script:
runs fine in script, generate 40k test1.jpg in the given directory
plotPlatesAlong - c(1:5)
plotDirPath - paste(dataPath, OutputPlots\\, sep=)
filename -
Most likely this is FAQ 7.22.
althogh how we are supposed to know what plotHeatTrellis does is beyond
me.
On Fri, 10 Feb 2006, [EMAIL PROTECTED] wrote:
Hello:
I have problems generating big JPG plots in functions. the following code
runs
fine in the script:
runs fine in script,
Peter Dalgaard wrote:
Gregor Gorjanc [EMAIL PROTECTED] writes:
Hello!
I would like to get MLE for parameter lambda of Poisson distribution. I
can use fitdistr() for this. After looking a bit into the code of this
function I can see that value for lambda and its standard error is
estimated
Hi,
I was trying to stack a topographic map readed in from esri shapefile with a
contour map and a vector map. However, the plotMap(maptools) and
contourplot(lattice) do not seem to work well on top of each other. Here is
part of my code.
xrange-range(225000:35)
Just to add a couple of thoughts to the previous suggestions.
If you really want output from things like summary.lm to have tabs instead of
spaces then you can type:
getAnywhere(print.summary.lm)
at the R prompt and it will show you the (not quite) source of code that R uses
to do the
I don't have a citation, but I think as long as the process is
stationary and not completely deterministic, the concept of a marginal
distribution is well defined and data from such a process will
eventially converge to that distribution. Of course, as the level of
dependence
DR == David Ruau [EMAIL PROTECTED] writes:
DR Hi, I was wondering if there were a option to have command
DR line completion in R like in the Bash shell?
Of course, via Emacs and ESS :-)
For ESS see: http://ess.r-project.org/
For Emacs see: http://www.gnu.org/software/emacs/emacs.html
Hi, Doug and Chris:
I just got the same error message with the lmList example in lme4:
(fm1 - lmList(breaks ~ wool | tension, warpbreaks))
Call:
Error in getResponseFormula(form) : Form must be a two sided formula
sessionInfo()
R version 2.2.1, 2005-12-20, i386-pc-mingw32
Hello R users,
I am trying to run a model with a binary response variable (nesting
success: 0 failure, 1 success) and 8 fixed terms. Nesting success was
examined in 72 cases in 34 territories (TER) during a 6 study years.
Territories are nested within 14 patches (PATCH). I want to run a model
At 11:42 10.02.2006 -0300, CENDOYA, Maria Gabriela wrote:
X-UIDL: PE5!!DnS!!np3!Sl2!
Hi Tom:
May be I didn't understand your question but, what I do to cut
and paste results from say summary.lm, in a word processors without losing
the nice shape of the R Console, is to choose the
Dear R community,
I'm trying to figure out if there is any way to split the printed output
of some commands, for example summary.lme, so that I can intersperse
comments in Sweave. I don't mind running the command numerous times and
masking various portions of the output, or saving the output as
At 09:25 10.02.2006 -0700, Gregory Snow wrote:
Just to add a couple of thoughts to the previous suggestions.
If you really want output from things like summary.lm to have tabs instead
of spaces then you can type:
getAnywhere(print.summary.lm)
at the R prompt and it will show you the (not
Hi,
I'm doing robust regression with the following command
rlm(dip~ind1+ind2-1,method=M,psi=psi,maxit=1000,acc=1e-15)
now when I ask for a summary
summary(rlm(dip~ind1+ind2-1,method=M,psi=psi,maxit=1000,acc=1e-15))
I get
Coefficients:
Value Std. Error t value
ind1 -0.0377
1. ?options (Look for scipen)
2. std error = Value/t value
3. Why do you care with these values?! Seems silly. It's probably off by
more than you're worrying about due to numerical error, deviation from
assumptions, arbitrariness of tuning constants, etc.
-- Bert Gunter
Genentech Non-Clinical
On Sat, 11 Feb 2006, Andrew Robinson wrote:
Dear R community,
I'm trying to figure out if there is any way to split the printed output
of some commands, for example summary.lme, so that I can intersperse
comments in Sweave. I don't mind running the command numerous times and
masking
At 06:09 10.02.2006 -0500, Doran, Harold wrote:
I didn't follow this thread entirely, but I did make a LaTeX
recommendation and I know that wasn't what you were asking for. But, if I
may, let me respond to the ideas you present below in an attempt to be
somewhat persuasive.
No, you are
Hi,
On Fri, 10 Feb 2006 10:56:32 -0800
Berton Gunter [EMAIL PROTECTED] wrote:
1. ?options (Look for scipen)
I'll check it
2. std error = Value/t value
3. Why do you care with these values?! Seems silly. It's probably off by
more than you're worrying about due to numerical error,
Thank you all for very useful and interesting responses. After reading the
comments and after some experiments, I added the following to a text I will
be handing out to the students (I would not mind comments):
QUOTEThe contents of the text output from R may be very sophisticated,
but the
Roger,
that's beautiful -- thanks!
Andrew
On Sat, February 11, 2006 5:58 am, Roger Bivand said:
On Sat, 11 Feb 2006, Andrew Robinson wrote:
Dear R community,
I'm trying to figure out if there is any way to split the printed output
of some commands, for example summary.lme, so that I can
On Fri, 10 Feb 2006 20:29:40 +0100 Angelo Secchi wrote:
Hi,
On Fri, 10 Feb 2006 10:56:32 -0800
Berton Gunter [EMAIL PROTECTED] wrote:
1. ?options (Look for scipen)
I'll check it
2. std error = Value/t value
3. Why do you care with these values?! Seems silly. It's
At 17:10 09.02.2006 -0600, you wrote:
Tom Backer Johnsen wrote:
There has been an incredible number of responses in a short time, with a
number of different suggestions. With hindsight, I must admit I have not
been quite clear, so additional (somewhat lengthy) explanation is needed.
I want to
ListeRs,
Within the last two months, I thought I saw mention of an R function
that would create a new data frame composed of duplicates or multiple
copies of rows of an input data frame given one or several columns of
values indicating how many times each row should be copied. As a simple
Eric Archer wrote:
ListeRs,
Within the last two months, I thought I saw mention of an R function
that would create a new data frame composed of duplicates or multiple
copies of rows of an input data frame given one or several columns of
values indicating how many times each row should
I'm using classification trees for the first time.
I understand the difference between these 2 packages, but I'm having a
bit of trouble interpreting the results.
I have 3 different response variables, but I'll only use 1 in this
discussion.
I first ran Tree. I was happy with the results, 6
Dear R helper,
I mange to transform uniform sequences to mixture
normal distributions using the following cods:
you forgot to mention the important fact that you are working
with package nor1mix (of which I am the maintainer which you
could have seen from library(help = nor1mix) or
Luis Ridao Cruz wrote:
R-help,
I'm using the following code:
plot(1985:2005, data, type = h, lwd = 10)
which should plot histogram like figure.
The bar tops are square and I was wondering if
it is intended to be so or not (I can of course
use barplot instead).
For the line endings
Martin == Martin Maechler [EMAIL PROTECTED]
on Fri, 10 Feb 2006 22:20:11 +0100 writes:
Dear R helper, I mange to transform uniform sequences to
mixture normal distributions using the following cods:
Martin you forgot to mention the important fact that you
Martin are
Division of Research at Kaiser Permanente, Oakland, CA
Seeking a Statistician/Senior SAS and R Programmer
Schedule: Full-time regular, 40 hours/week, Monday through Friday
Send Resume to: [EMAIL PROTECTED]
Position Summary:
Under general supervision, this position exists to assist the
Hello!
I would like to use lmer() to fit data, which are some estimates and
their standard errors i.e kind of a meta analysis. I wonder if weights
argument is the right one to use to include uncertainty (standard
errors) of data into the model. I would like to use lmer(), since I
would like
Given that this may very well be the most common use of the
R2HTML package I wonder if the R2HTML package developer would
be interested in providing an HTML2clip convenience wrapper
as part of the R2HTML package like this:
HTML2clip - function(x, file. = file(clipboard, w), append = FALSE,
On 2/10/06, Tom Backer Johnsen [EMAIL PROTECTED] wrote:
At 17:10 09.02.2006 -0600, you wrote:
Tom Backer Johnsen wrote:
There has been an incredible number of responses in a short time, with a
number of different suggestions. With hindsight, I must admit I have not
been quite clear, so
Hi all,
I'm new to the list...pretty new at learning to code in R...
Is there a way to relate 2 different arrays in R?
Hypothetical example:
data1
IDz
1100
2250
375
412
589
data2
IDz
1
1
1
1
2
3
4
3
4
5
5
5
etc.
Goal is to fill column z in data2 with
thanks, it works.
best, Alex
On Fri, 10 Feb 2006, Prof Brian Ripley wrote:
Most likely this is FAQ 7.22.
althogh how we are supposed to know what plotHeatTrellis does is beyond
me.
On Fri, 10 Feb 2006, [EMAIL PROTECTED] wrote:
Hello:
I have problems generating big JPG
Dear Jeff,
Assuming that the column named z in the matrix data2 already exists and
has arbitrary content (such as 0's or NA's), how about the following?
data2[,z] - data1[data2[,ID], z]
I hope this helps,
John
John Fox
Department of Sociology
McMaster
The following appears to be an answer to your question, though I'd be
pleased to receive critiques from others. Since your example is NOT
self contained, I modified an example in the glmmPQL help file:
(fit - glmmPQL(y ~ factor(week)-1+trt, random = ~ 1 | ID,
+
I use Mac. I installed R with the download from R-project, so R is in
/usr/bin/R.
My TclTk library is installed via fink (although I don't remember I
intentionally
installed it) so libtk8.4.dylib is in /sw/lib.
I tried to install the HDF5 package from within R:
install.packages('hdf5')
and
thanks, Andy and John, for suggestions.
worked great.
jeff
**
Jeffrey Moore, Ph.D.
Postdoctoral Research Scientist
Duke Center for Marine Conservation
Duke University Marine Laboratory
135 Duke Marine Lab Road
Beaufort, NC 28516
Phone: (252) 504-7653
Fax:
Hi,
I have a set of data in x,y coordinates across the range of -5 to 5 in each
dimension. I would like to obtain the frequency distribution of the
different points, and then graph them so you can see which of the points are
the most frequently occurring.
This would seem to be easy in Matlab,
78 matches
Mail list logo