On 7/12/07, Pete Kazmier [EMAIL PROTECTED] wrote:
hadley wickham [EMAIL PROTECTED] writes:
On 7/12/07, Pete Kazmier [EMAIL PROTECTED] wrote:
Is there a way in ggplot to make a histogram with the left-hand y-axis
label as frequency, and a right-hand y-axis label as percentage?
Not
Dear Friends,
I have done the analysis of SELDI
data using the PROcess package. And now, I would like to confirm
the analysis using CaMassClass package. I have downloaded and
tried to analyze the Test samples that provided, it shows
the following error.
directory = system.file(Test,
Dear R users,
I am trying to connect an Access database.
The code I use is :
channel-odbcConnectAccess(Base)
table.db-sqlFetch(Base,table,colnames=T)
When I look at the table in R, the table imported lacks the first 41 lines, and
the column names are the values of the 42nd record.
Any hints ?
Dear Professor Schaefer
I am experiencing a technical difficulty with your mix package.
I would appreciate it if you could help me with this problem.
When I run the following code, R 2.5.1 and R 2.6.0 crashes.
It's been tested on at least 2 windows machine and it is consistent.
Execution code
I did it a long time ago, so I do not remember why I have to use
exim4... sorry. mail did not work alone.. that is why I use exim4.
Perhaps a simpler solution exists.
Cheers
(Ted Harding) wrote:
On 12-Jul-07 16:10:46, Stéphane Dray wrote:
Here is a small function that I used on Debian.
Hi
I have dataframe which contain 5 columns and 1000 records. I want standard
each cell.
I want range each column between 0 and 1 . I think i must use loop?
could you help me?
-
Moody friends. Drama queens. Your life? Nope! - their life, your story.
This is an ODBC issue, not an R issue. Use a command-line ODBC client to
find a suitable SQL invocation to retrieve the data, and run that via
sqlQuery.
(It might be a quoting issue: Access is not using standard SQL.)
On Thu, 12 Jul 2007, Patrick Gonin wrote:
Dear R users,
I am trying to
Hi R help group,
Any budy give me tutorial/paper for eigen value calculation and PCA.
Thanking you.
--
Nitish Kumar Mishra
Junior Research Fellow
BIC, IMTECH, Chandigarh, India
E-Mail Address:
[EMAIL PROTECTED]
[EMAIL PROTECTED]
__
Gustaf Rydevik wrote:
Hi all,
I've run into a slightly illogical (to me) behaviour with the $
subsetting function.
consider:
Test
A B
1 1 Q
2 2 R
Test$A
[1] 1 2
vector-A
Test$vector
NULL
Test$A
[1] 1 2
Test[,vector]
[1] 1 2
Hi, I have got four correlation matrix. They are the same set of variables
under different conditions. Is there a way to test whether the correlation
matrix are significently different among each other? Could
anyone give me some advice?
--
View this message in context:
Hi,
Again, screen was described in the April issue of R News.
Augmenting R with Unix Tools by Andrew Robinson.
@ARTICLE{Rnews:Robinson:2007,
AUTHOR = {Andrew Robinson},
TITLE = {Augmenting {R} with {Unix} Tools},
JOURNAL = {R News},
YEAR = 2007,
VOLUME = 7,
NUMBER = 1,
PAGES =
On 7/12/07, JVVerkuilen [EMAIL PROTECTED] wrote:
Hello all,
I am using lme4 to fit some mixed logistic regressions. I need to
impose an identification constraint of the following form:
(1sig12)
(sig12 sig22)
and have not figured out how to do it, i.e., sig11 = 1 but the rest of
Hi all,
How does one create a matrix of values in a loop. For example, I have 46*70
matrix with columns giving the values for some blood tests performed upon
the 46 patients(rows). I want to perform a simple test of normality
(shapiro-wilk test) upon each of the biochemical tests for the 46
you could also have a look at function posdefify() from package
`sfsmisc'.
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, Belgium
Tel: +32/(0)16/336899
Fax:
Thanks.
I am working (making some progress) on two installations (different sun's).
I am working to get the system admins to install the right libraries and
compilers to get the job done.
Cheers,
Dan
=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=-=
Daniel A. Powers, Ph.D.
Department of Sociology
You can control the panel sequence with subscripting and transpose.
Here are several examples. I think tmp.tr3 is the one you asked for.
library(lattice)
tmp - data.frame(x=rnorm(24), y=rnorm(24), a=rep(letters[1:6],4),
b=rep(LETTERS[1:4],each=6))
tmp.tr - xyplot(y ~ x | a*b, data=tmp)
tmp.tr
spime wrote:
Suppose I have
Training data: my.train
Testing data: my.test
The bootstrap does not need split samples.
I want to calculate bootstrap error rate for logistic model. My wrapper
function for prediction
pred.glm - function(object, newdata) {
ret -
Suppose I have
Training data: my.train
Testing data: my.test
I want to calculate bootstrap error rate for logistic model. My wrapper
function for prediction
pred.glm - function(object, newdata) {
ret - as.factor(ifelse(predict.glm(object, newdata,
type='response') 0.4, 0, 1))
Balazs Torma wrote:
Thank you all for your answers!
The problem is that I don't know the length of the list in advance!
And hoped for a convinience structure which reallocates once the
preallocated list (or matrix) becomes full.
That's not massively hard to do yourself, is it?
As in
Hello,
I'm using the weighted nonlinear regression proposed in nls help.
I would like to know which is proper reference to the Weighted version of
Michaelis-Menten model. Which is The White book cited in nls help?
And, by-the-way, many thanks to Martin Maechler for this function.
Best
Dear all,
Has someone implemented in R (or any other language)
Knol DL, ten Berge JMF. Least-squares approximation of an improper correlation
matrix by a proper one. Psychometrika, 1989, 54, 53-61.
or any other similar algorithm?
Best regards
Jens Oehlschlägel
Background:
I want to
On Fri, 13 Jul 2007, Philippe Grosjean wrote:
If all the data coming from your iterations are numeric (as in your toy
example), why not to use a matrix with one row per iteration? Also, do
preallocate the matrix and do not add row or column names before the end
of the calculation. Something
Hi,
Using library(lattice), is there any way to tell xyplot to plot
panels top to bottom, then left to right (i.e. panels are appended
vertically, then horizontally). as.table changes the plot direction
from left-to-right then top-to-bottom, to right-to-left then bottom-
to-top, but that's
You may find screen useful, see e.g.
http://tolstoy.newcastle.edu.au/R/e2/help/07/02/10824.html
for a short description how to use it with R (reattaching works via
screen -r). (Maybe your sysadmin has to install screen first, it's a
unix tool.)
KR
Martin
[EMAIL PROTECTED] wrote:
Try using the unix command screen
On 13/07/07, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote:
Goodmorning everybody,
I need to run an R program via SSH. Usually I open R, I run the
program and I stay logged-in, waiting for the output. As a matter of
fact, if I close the connection with SSH I
another solution:
l - vector(mode = list, length = 10^5)
system.time(for(i in (1:10^5)) l[[i]] - c(i,i+1,i))
On my system this version is even slightly faster than the matrix version ...
Best,
Matthias
- original message
Subject: Re: [R] filling a list faster
Sent: Fri, 13 Jul
You don't need to do a loop. You can do it with an
apply statement I think. See ?apply
If the matrix is mat to do the test by column try:
apply(mat, 2, shapiro.test)
--- Tavpritesh Sethi [EMAIL PROTECTED] wrote:
Hi all,
How does one create a matrix of values in a loop.
For example, I have
See this:
http://tolstoy.newcastle.edu.au/R/help/06/08/32543.html
and also the other messages in that thread.
On 7/13/07, Yan Wong [EMAIL PROTECTED] wrote:
Hi,
Using library(lattice), is there any way to tell xyplot to plot
panels top to bottom, then left to right (i.e. panels are appended
Hi Karl,
There's no official way to do it, but you can hack the colour
gradient scale to do what you want:
x=-10:10
y=-10:10
dat=expand.grid(x=x,y=y)
dat$z=dat$x^2+dat$y^2-100
# Create a modified scale
gr - scale_fill_gradient2()$clone()
gr$breaks - function(.) seq(-100, 100, by=10)
I have a model with 5 parameters that I am optimising where the (best)
value of the objective function is negative. I would like to use the
Hessian matrix (from genoud and/or optim functions) to construct the
covariance and correlation matrices.
This is the code that I am using:
est -
If all the data coming from your iterations are numeric (as in your toy
example), why not to use a matrix with one row per iteration? Also, do
preallocate the matrix and do not add row or column names before the end
of the calculation. Something like:
m - matrix(rep(NA, 3*10^5), ncol = 3)
Thank you all for your answers!
The problem is that I don't know the length of the list in advance!
And hoped for a convinience structure which reallocates once the
preallocated list (or matrix) becomes full.
By the way, what is your very fast machine, that is actually four
times faster
hello,
first I create a list:
l - list(1-c(1,2,3))
then I run the following cycle, it takes over a minute(!) to
complete on a very fast mashine:
for(i in (1:10^5)) l[[length(l)+1]] - c(i,i+1,i)
How can I fill a list faster? (This is just a demo test, the elements
of the list are
On Fri, 13 Jul 2007, David Barron wrote:
Try having a look at the scale and sweep functions.
sweep applies to arrays, not data frames, and scale converts to a matrix.
For a data frame
df2 - df1
df2[] - lapply(df1, function(x) {r - range(x, na.rm=TRUE);
Gustaf Rydevik gustaf.rydevik at gmail.com writes:
Is there a reason for the $ operator not evaluating
the vector before executing?
Yes. Consider:
x - list(a = 1, b = 'n')
a - 'b'
x[[a]] == 'n'
[1] TRUE
x[['a']] == 1
[1] TRUE
x$a == ???
If the result is 'n', then all bare $ references
Hi all,
I've run into a slightly illogical (to me) behaviour with the $
subsetting function.
consider:
Test
A B
1 1 Q
2 2 R
Test$A
[1] 1 2
vector-A
Test$vector
NULL
Test$A
[1] 1 2
Test[,vector]
[1] 1 2
Is there a reason for the $ operator not evaluating the vector before executing?
On 12-Jul-07 23:47:39, Masanao Yajima wrote:
Dear Professor Schaefer
I am experiencing a technical difficulty with your mix package.
I would appreciate it if you could help me with this problem.
When I run the following code, R 2.5.1 and R 2.6.0 crashes.
It's been tested on at least 2
Try having a look at the scale and sweep functions.
David
On 13/07/07, Amir_17 [EMAIL PROTECTED] wrote:
Hi
I have dataframe which contain 5 columns and 1000 records. I want standard
each cell.
I want range each column between 0 and 1 . I think i must use loop?
could you help me?
How do I make Lithuanian characters display correctly in R graphics?
Instead of the special characters for Lithuanian language I get question
marks...
I use Ubuntu Feisty, the locale is utf-8 ...
Do I need to specify somewhere the locale for R, or - default font for the
graphics?
--
Donatas
one way is the following (maybe there're better):
pats - do.call(paste, c(as.data.frame(M), sep = \r))
pats - factor(pats, levels = unique(pats))
cbind(unique(M), Freq = as.vector(table(pats)))
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School
Thank you Pr Ripley, it worked well with a simple sqlQuery rather than a
sqlFetch command:
I paste it here in case some other users encounter a similar problem:
table.db-sqlQuery(channel,paste(select * from table))
Dr Patrick Gonin
Institut Gustave Roussy
94800 Villejuif
France
Prof Brian
Dear r-help,
I would like to use the 'ppm' function of the 'spatstat' package to
fit a Strauss inhibition model. I understand that I can specify a
parametric model for the background trend, but how would I specify a
trend which is estimated using a Kernel density smoother?
In particular, I
On Fri, 13 Jul 2007, Donatas G. wrote:
How do I make Lithuanian characters display correctly in R graphics?
Instead of the special characters for Lithuanian language I get question
marks...
I use Ubuntu Feisty, the locale is utf-8 ...
Do I need to specify somewhere the locale for R, or -
Please could you write an example of this command:
table.db-sqlQuery(channel,paste(select * from table))
With real names of .mdb file and table?
Thanks
Michael
- Original Message -
From: Patrick Gonin [EMAIL PROTECTED]
To: Prof Brian Ripley [EMAIL PROTECTED]
Cc: r-help@stat.math.ethz.ch
Hi All,
convhulln {geometry} computes the convex hull of a set of points in n-
dimensions via .call, returning the hull itself, and also
unavoidably generates a diagnostic report on an Rterm console. See
the example below.
I need to access the results of the diagnostic report (specifically
Donatas G. wrote:
How do I make Lithuanian characters display correctly in R graphics?
Instead of the special characters for Lithuanian language I get question
marks...
I use Ubuntu Feisty, the locale is utf-8 ...
Do I need to specify somewhere the locale for R, or - default font for the
Goodmorning everybody,
I need to run an R program via SSH. Usually I open R, I run the
program and I stay logged-in, waiting for the output. As a matter of
fact, if I close the connection with SSH I loose the calculations and
the output of my R program. What command I have to use in order to
Here is what works for me:
base.mdb being an Access database in the R default data path, Piece being a
table in base:
channel - odbcConnectAccess(base)
piece.db-sqlQuery(channel,paste(select * from Piece))
Patrick
Michal Kneifl [EMAIL PROTECTED] wrote: Please could you write an example of
It all depends on what you want to do. In your example, it is faster
to first fill in a matrix and then convert the matrix to a list. The
problem with filling in the list is that you are dynamically
allocating space for each iteration which is probably taking at least
an order of magnitude more
At 17:17 12/07/2007, kdestler wrote:
Hi
I recently imported data from r into Stata. I then ran the linear
regression model I've been working on, only to discover that the results are
somewhat (though not dramatically different). the standard errors vary more
between the two programs than do the
Actually if you are really interested in the list, then just do the
lapply and compute your data; it seems to be even faster than the
matrix:
system.time(l.1 - lapply(1:10^5, function(i) c(i, i+1, i)))
user system elapsed
0.500.000.61
l.1[1:4]
[[1]]
[1] 1 2 1
[[2]]
[1] 2 3 2
Hi
I imported a lot of text-files with R and saved them by using:
save(list=ls(),file=Grunddaten.Rdata)
After that I wanted to check my saved data and wrote
source(Grunddaten.Rdata)
This gives me an error:
Error in parse(file, n = -1, NULL, ?) : Syntax error at 1:
Can
A seemingly simple problem has me stumped. Is it possible to choose the
number of colour breaks for a gradient scale in the current version of
ggplot2?
Here is a simple example:
-
x=-10:10
y=-10:10
dat=expand.grid(x=x,y=y)
dat$z=dat$x^2+dat$y^2-100
Prof Brian Ripley wrote:
On Fri, 13 Jul 2007, Peter Dalgaard wrote:
The correct incantation seems to be
postscript(font=URWHelvetica, encoding=ISOLatin7)
plot(0,main=tolower(\u104\u116\u0118\u012e\u0172\u016a\u010c\u0160\u017d))
dev.off()
The encoding should happen automagically
I need help regarding to the following problem:
Consider this matrix:
M - matrix(c(1,2, 4,3, 1, 2),3, 2, byrow=TRUE)
M
[,1] [,2]
[1,]12
[2,]43
[3,]12
I would like to have a matrix which counts the identical rows and places the
counts into its third column. I
On Fri, 13 Jul 2007, Peter Dalgaard wrote:
Donatas G. wrote:
How do I make Lithuanian characters display correctly in R graphics?
Instead of the special characters for Lithuanian language I get question
marks...
I use Ubuntu Feisty, the locale is utf-8 ...
Do I need to specify somewhere
DR == Dimitris Rizopoulos [EMAIL PROTECTED]
on Fri, 13 Jul 2007 14:43:08 +0200 writes:
DR you could also have a look at function posdefify() from
DR package `sfsmisc'.
DR I hope it helps.
Yes, thanks, Dimitris; note that my posdefify() function uses a
pretty arbitrary fudge
Martin,
I sent you the Matlab code for this which I had obtained from Nich Higham.
Cheng, Sheung Hun and Higham, Nick (1998)
A Modified Cholesky Algorithm Based on a Symmetric Indefinite
Factorization;
\emph{SIAM J. Matrix Anal.\ Appl.}, \bold{19}, 1097--1110.
Do you remember?
Ravi.
Caskenette, Amanda wrote:
I have a model with 5 parameters that I am optimising where the (best)
value of the objective function is negative. I would like to use the
Hessian matrix (from genoud and/or optim functions) to construct the
covariance and correlation matrices.
This is the code
Thank you very much for the quick response,
Yes that is correct, some values for V are negative.
I wish that looking at the H, eigen(H), V, s, s%o%s made me wiser, however I
really have no idea what I am looking for.
The genoud function does not converge if left running for a week so I have a
Dear All,
I have a data frame with the columns Mass and Intensity (this is mass
spectrometry stuff). Each of the mass values gives rise to a mass window of
5 ppm around the individual mass (from mass - mass/1E6*5 to mass +
mass/1E5*5). I need to filter the array such that in case these mass
On Fri, 13 Jul 2007, Peter Dalgaard wrote:
Prof Brian Ripley wrote:
On Fri, 13 Jul 2007, Peter Dalgaard wrote:
The correct incantation seems to be
postscript(font=URWHelvetica, encoding=ISOLatin7)
plot(0,main=tolower(\u104\u116\u0118\u012e\u0172\u016a\u010c\u0160\u017d))
dev.off()
The
On 7/13/07, Richard M. Heiberger [EMAIL PROTECTED] wrote:
You can control the panel sequence with subscripting and transpose.
Here are several examples. I think tmp.tr3 is the one you asked for.
library(lattice)
tmp - data.frame(x=rnorm(24), y=rnorm(24), a=rep(letters[1:6],4),
[EMAIL PROTECTED] wrote:
Hello, R experts,
Sorry for asking this question again again since I really want a help!
I have a two-factor experiment data and like to calculate estimates of
interation contrasts say factor A has levels of a1, a2, and B has
levels of b1, b2, b3, b4, and b5 with 3
See help(legend) and help(identify).
Ajay Singh wrote:
Hi,
I have two problems in R.
1. I need 10 cdfs on a graph, the graph needs to have legend. Can you let
me know how to get legend on the graph?
2. In ecdf plot, I need to know the x and y co-ordinates. I have to get
Sorry, this sounds like a fairly basic question that can be resolved by
which() and possible ifelse(). There is no details in your email.
I am afraid you have to learn the basics of R or ask question with more
details (e.g. example data).
Or ask someone locally.
Regards, Adai
Johannes
I think we need a bit more information and perhaps a
small example data set to see what you want.
I am not familiar with term mass window. Is this a
confidence interval around the mass value?
--- Johannes Graumann [EMAIL PROTECTED]
wrote:
Dear All,
I have a data frame with the columns
--- Renger van Nieuwkoop [EMAIL PROTECTED] wrote:
Hi
I imported a lot of text-files with R and saved them
by using:
save(list=ls(),file=Grunddaten.Rdata)
After that I wanted to check my saved data and wrote
source(Grunddaten.Rdata)
This gives me an error:
Error in
Hello
I am running a simulation study to assess the type I error rate of a likelihood
ratio test. This requires me to minimize over restricted spaces under the null
and under the alternative hypothesis. My problem is that occassionally I run
into an error message that reads: Error in
Hi,
I want to know how to increase the limit size of R object. For example, i
want to store a matrix (1 rows and 8000 columns)into an object x.
However, it gives me an error message 'can't allocate a vector of size
396.7 Mb'. But there is no problem to store a matrix of 5000*3000 matrix
into
Based on the feedback received, I did the following:
a) moved my lib sub-directory from the existing installed R version to
c:\myRLib
b) installed the updated R version
c) created .Renviron file in the home directory (C:\R-2.5.1) with the line
R_LIBS=c:/myRLib
d) used .libPaths() command to
Dear list,
I have a dataset consists of duplicated sequences within day for each patient
(see below data) and I want to reshape the data with patient as time variable.
However the reshape function only takes the first sequence of the replicates
and ignores the second. How can I 1) average the
On 13 Jul 2007, at 17:11, [EMAIL PROTECTED] wrote:
Another high level option is to change the rule determining how
packets are chosen for a given panel in the layout.
print(tmp.tr3,
packet.panel = function(layout, row, column, ...) {
layout - layout[c(2, 1, 3)]
Hi Tom,
I have a dataset consists of duplicated sequences within day for each
patient (see below data) and I want to reshape the data with patient as time
variable. However the reshape function only takes the first sequence of the
replicates and ignores the second. How can I 1) average
Dimitris,
Thanks a lot for the quick response with the pointer to posdefify. Using its
logic as an afterburner to the algorithm of Higham seems to work.
Martin,
Jens, could you make your code (mentioned below) available to the community,
or even donate to be included as a new method of
Hi all.
How do I extract date from fcs format file with R. I.e I'd like
make statistical analysis using R-program, but I don't know if there
are R-packets for fcs format file, and using examples.
Thanks.
Pta: In Linux SO exist any program that transform from fcs format to
ASCII text file?
Hi:
We're trying to install the Hmisc package on a Solaris 9 machine.
Here's what we get:
R CMD INSTALL /usr/local/srctree/Hmisc_3.4-2.tar.gz
* Installing to library '/usr/local/lib/R/library'
* Installing *source* package 'Hmisc' ...
** libs
g95 -fPIC -g -O2 -c cidxcn.f -o cidxcn.o
g95
This sounds like a solution I've been looking for. With this setup now
in place, when you download and install some new packages, where will
they be put? Into C:\myRlib ?
Thanks.
--Chris
Christopher W. Ryan, MD
SUNY Upstate Medical University Clinical Campus at Binghamton
40 Arch Street,
Hello, Chuck,
Thank you very much for your help! But the contrasts I want to do
simutaneously is
contrasts(B)
[,1] [,2] [,3] [,4]
b1 -4 -3 -2 -1
b21 -3 -2 -1
b312 -2 -1
b4123 -1
b51234
Could you please show me how to calculate
there is rflowcyt (older and on its way out) and flowCore, which
contains a newer tool set
available at
www.bioconductor.org
Horacio Castellini wrote:
Hi all.
How do I extract date from fcs format file with R. I.e I'd like
make statistical analysis using R-program, but I don't know if there
This is not related to R but I was hoping that someone could help me. I
am reading the Understanding the Metropolis Hastings Algorithm
paper from the American Statistician by Chip and Greenberg, 1995, Vol
49, No 4. Right at the beginning they explain the algorithm for basic
acceptance rejection
On Friday 13 July 2007 17:56:45 Prof Brian Ripley wrote:
On Fri, 13 Jul 2007, Peter Dalgaard wrote:
Prof Brian Ripley wrote:
On Fri, 13 Jul 2007, Peter Dalgaard wrote:
The correct incantation seems to be
postscript(font=URWHelvetica, encoding=ISOLatin7)
-- Forwarded message --
From: John C Frain [EMAIL PROTECTED]
Date: 13-Jul-2007 22:30
Subject: Re: [R] THANK YOU: Updating R version
To: Christopher W. Ryan [EMAIL PROTECTED]
When I update R the following has worked for me (Windows XP)
1. Install the new version to a new
Hello everybody
I have done a generalized linear model (glm-function) with a binary
response variable (presence/absence) and two explanatory variables
(including the interaction term (model-glm(x~factor1*factor2,
family=quasibinomial)). The interaction is what we are mainly
interested in
Hello All,
I wanted to do many plots (in my case, wanted to get 6 histograms) on the
same figure. Is there a method in R that analogous to 'subplot' in MATLAB?
Any help will be very much appreciated.
thanks,
Suman
[[alternative HTML version deleted]]
Jens,
An alternative solution to the improper matrix
problem is to do a principal factor solution
rather than a maximum likelihood factor analysis
solution. In the following discussion, I am
using the factor.pa (principal factor) function
from my psych package.
Using your test set
Could someone show me how to get a blue line in this plot?
ggplot(movies, aes(x=rating)) + stat_qq(geom=line,
quantiles=seq(0,1,0.005), distribution=qunif)
I've tried many permutations but cannot seem to find the right
combination. I've tried these flavors:
ggplot(movies, aes(x=rating))
On 7/13/07, suman Duvvuru [EMAIL PROTECTED] wrote:
Hello All,
I wanted to do many plots (in my case, wanted to get 6 histograms) on the
same figure. Is there a method in R that analogous to 'subplot' in MATLAB?
Here are a few possibilities:
data(singer, package = lattice)
## using
On Fri, 13 Jul 2007, Leeds, Mark (IED) wrote:
This is not related to R but I was hoping that someone could help me. I
am reading the Understanding the Metropolis Hastings Algorithm
paper from the American Statistician by Chip and Greenberg, 1995, Vol
49, No 4. Right at the beginning they
On Fri, 13 Jul 2007, simone wrote:
Hello everybody
I have done a generalized linear model (glm-function) with a binary
response variable (presence/absence) and two explanatory variables
(including the interaction term (model-glm(x~factor1*factor2,
family=quasibinomial)). The interaction is
Let's say I have a list called the_list consisting of three components:
the_list$component_1
the_list$component_2
the_list$component_3
Now, I want to access it using a variable called comp.
comp - component_1
I'm looking for some function that let's me do this:
unknown_function(the_list, comp)
I feel like I'm not firing at all cylinders with the R language,
despite reading over the R Language Definition.
Has anyone seen something like an R for Rubyists guide (i.e.
teaching the R language for those who are more familar with Ruby)?
Alternately, a book with lots of examples about how
Hi All,
I want to upload J K Lindsey's repeated in a LINUX OS..
I had tried the command..
[EMAIL PROTECTED] Desktop] # R CMD INSTALL repeated.gz
WARNING: invalid package 'repeated'
*Installing to library
Do you really have repeated.gz? Lindsey usually packs his packages with
the .tgz extension, e.g. at
http://popgen.unimaas.nl/~jlindsey/rcode.html
You need to get repeated.tgz and rmutil.tgz and use
tar zxf rmutil.tgz
tar zxf repeated.tgz
R CMD INSTALL rmutil repeated
Optional (but highly
Dear useRs,
I am trying to model a time series with a transfer function. I think
it can be put into the ARMA framework, and estimated with the 'arima'
function (and others have made similar comments on this list). I have
tried to do that, but the results have so far been disappointing.
Maybe I am
Hi,
is there an R function to convert unixtime to one of the R time
formats (using chron or POSIXct)?
Example data:
unixtime year month day hour
1183377301 2007 7 2 13
I would like to only use the first column.
If such a function does not exist: What would
On Sat, 14 Jul 2007, Patrick Drechsler wrote:
Hi,
is there an R function to convert unixtime to one of the R time
formats (using chron or POSIXct)?
Example data:
unixtime year month day hour
1183377301 2007 7 2 13
I would like to only use the first
when producing boxplot from bwplot, I have five groups: Nitrogen, Duration,
Pressure, A, Z. I wish the graphical display is according to the original
order. But the R-function bwplot seems to automatically adjust the groups
according to the alphabetic oder and thus creat a graph for A, Duration,
98 matches
Mail list logo