On Thu, 31 Aug 2006, Petr Pikal wrote:
> Hi
>
> or you can change
>
> Rprofile.site in etc directory to whatever startup commands you want
> to execute to have the same profile in all sessions.
Please see the warning about this in ?Startup: it is not the same thing as
putting .Rprofile in you
Hello
I have 2 questions:
1. How I can do the test of hartley for homogeneity of variances?
2. How I can do the test of Box for homogeneity of variances?
Thanks in advance
--
Marco Vinisio Martinez
Profesor Departamento de Matematicas y Departamento de Biologia
Pontificia Universidad Javeriana
Hi All,
I'm trying to download data using following code.
require(UsingR) ## This is the R-package you need
to run command yahoo.get.hist.quote
#Initialize empty table
closing <-NULL
#Downalod consituents since I don't have it on my comp
download.file("http://www2.standar
Hi
or you can change
Rprofile.site in etc directory to whatever startup commands you want
to execute to have the same profile in all sessions.
HTH
Petr
On 30 Aug 2006 at 15:11, BBands wrote:
Date sent: Wed, 30 Aug 2006 15:11:35 -0700
From: BBands <[EMAIL PROT
Carlo Trimarchi a écrit :
> Hi,
> how can I read, using for example read.table() or scan(), just one
> column from a text file that has more columns without any header?
> Thanks, bye.
afaik, you have to read all the table
and then you select the column you want.
eg read.table(blabla)[3] to get th
'merge' is the key here. You say you want to merge, but it seems did not
try merge()
> (res <- merge(cvsFrame, downloadsFrame, by="row.names", all=TRUE))
Row.names cvsactions downloads
1 2002-11-15 412
2 2002-12-15 9 8
3 2003-01-15 5NA
4 2003
If you are converting them to 'its' anyways then after the
conversion to 'its' use the 'its' union command.
On 8/31/06, James Howison <[EMAIL PROTECTED]> wrote:
> Hi,
>
> I've examined the archives and found quite a few questions on
> concatenating dataframes, but none that really addressed my iss
Hi,
I've examined the archives and found quite a few questions on
concatenating dataframes, but none that really addressed my issue,
I'm afraid. I've also examined the cbind and rbind documentation but
nonetheless here I am writing to r-help ;)
This is what I have (the row names are dates
In thinking about this a bit more we can use
panel.superpose/panel.groups to shorten it:
# define data -- df
# note that your val2 and val3 lines had a syntax
# so we have commented them out and
# replaced them as shown.
n <- 18
x1 <- seq(1,n)
val1 <- -2*x1+50
# val2 <- (-2*(x1-8)2)+100
val2 <- (
Dear useRs,
A few days ago, the initial version 0.1.0 of a new package 'random' was
uploaded to CRAN.
The random packages provides convenient access to the non-deterministic
random numbers provided by the random.org site created by Mads Haahr
(http://www.random.org).
While certain hardware an
Or maybe this is what you are looking for where pnl below was
created by modifying source to the panel.plot.default in the zoo
package (there might be a simpler way):
pnl <- function (x, y, subscripts, groups, col, pch, type, ...) {
for (g in levels(groups)) {
idx <- g == groups[subsc
Hi all
I'm getting a NaN returned on using dffits, as explained
below. To me, there seems no obvious (or non-obvious reason
for that matter) reason why a NaN appears.
Before I start digging further, can anyone see why dffits
might be failing? Is there a problem with the data?
Consider:
#
Rolf
> The last invocation indicated that the name of the file was *really*
> ``.Rprofile.txt'' --- although I'd tried to save it as (simply)
> ``.Rprofile''. Is that the problem? If so, how can I persuade
> Windoze NOT to stick that damned .txt tag on the end?
The easiest way is to use a smar
To handle conditioning on survey we provide a panel function
that subsets col and pch:
# define test data - df
# note that your val2 and val3 lines had a syntax
# so we have commented them out and
# replaced them as shown.
n <- 18
x1 <- seq(1,n)
val1 <- -2*x1+50
# val2 <- (-2*(x1-8)2)+100
val2 <-
you can read them all in and delete the ones you don't want. Or check
out 'what' on 'scan' or colClasses on 'read.table'
On 8/30/06, Carlo Trimarchi <[EMAIL PROTECTED]> wrote:
> Hi,
> how can I read, using for example read.table() or scan(), just one
> column from a text file that has more column
On 8/30/2006 5:49 PM, Rolf Turner wrote:
> I am (for my sins) having to do some work using R under Windoze. I
> wanted to set up a .Rprofile to control my set-up. The docs on
> .Rprofile say that it can/should be placed in ``the user's home
> directory''. ``An Introduction to R'' observes lucidl
> -Original Message-
> From: [EMAIL PROTECTED] [mailto:r-help-
> [EMAIL PROTECTED] On Behalf Of Rolf Turner
> Sent: Wednesday, August 30, 2006 2:50 PM
> To: r-help@stat.math.ethz.ch
> Subject: [R] .Rprofile under Windoze.
>
> I am (for my sins) having to do some work using R under Windoze.
Hi,
how can I read, using for example read.table() or scan(), just one
column from a text file that has more columns without any header?
Thanks, bye.
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read th
And... If you have/use shortcuts to R, you may also save an
".Rprofile" to whatever directory you name in the "Start in:" field of
the shortcut. This allows one to have many profiles.
jab
--
John Bollinger, CFA, CMT
www.BollingerBands.com
If you advance far enough, you arrive at the begin
> In fact, my problem is to fit the data for every level of the f2 factor,
> showing the levels of the f1 factor and that for several surveys .
> Here's an example closer to my actual data :
Then maybe you want:
qplot(x, y, . ~ surveys, data=df, type="line", colour=f1, id=f2, size=f2)
(which does
On 8/30/06, Rolf Turner <[EMAIL PROTECTED]> wrote:
> I am (for my sins) having to do some work using R under Windoze. I
> wanted to set up a .Rprofile to control my set-up. The docs on
> .Rprofile say that it can/should be placed in ``the user's home
> directory''. ``An Introduction to R'' obser
Under Windows mine is located here
C:\Program Files\R\R-2.3.1\library\base\R
The file name, however is not .Rprofile, but rather Rprofile
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: [EMA
Gabor Grothendieck a écrit :
>Note that before entering this you need:
>
>library(lattice)
>library(grid) # to access the viewport function
>
>On 8/29/06, Gabor Grothendieck <[EMAIL PROTECTED]> wrote:
>
>
>>Try this:
>>
>>xyplot(val ~ x, data = df, type = "p",
>> col = as.numeric(df$f1), p
Gabor Grothendieck a écrit :
> Try this:
>
> xyplot(val ~ x, data = df, type = "p",
> col = as.numeric(df$f1), pch = as.numeric(df$f2))
>
> key1 <- list(border = TRUE, colums = 2, text = list(levels(df$f1)),
> points = list(pch = 1:nlevels(df$f1))
> )
>
> key2 <- list(border = TRUE, colums
I am (for my sins) having to do some work using R under Windoze. I
wanted to set up a .Rprofile to control my set-up. The docs on
.Rprofile say that it can/should be placed in ``the user's home
directory''. ``An Introduction to R'' observes lucidly that this
concept needs to be clarified under W
hadley wickham a écrit :
>> I would like to use the lattice library to show several groups on
>> the same graph. Here's my example :
>>
>> ## the data
>> f1 <- factor(c("mod1","mod2","mod3"),levels=c("mod1","mod2","mod3"))
>> f1 <- rep(f1,3)
>> f2 <- factor(rep(c("g1","g2","g3"),each=3),level
I tried to run gdb in linux with emacs
But could not even run a simple example in the writing extensions
tutorial.
The execution history is as follows.
Gdb worked fine for other debugging such as C++ codes.
Thanks
Han
I started R at echo of emacs by typing
(also tried other methods mentioned
There are functions to do weighted summary statistics in the Hmisc
package (wtd.quantile, ...).
For more complicated analyses (but not plots yet) the biglm package has
a bigglm function that expects the data in chunks, you could write a
function that expand parts of the dataset at a time.
Hope th
You may want to look at the logspline package, it uses a different
technique than density does, but it estimates densities and allows you
to tell the routine that there is a minimum value and that the density
does not extend beyond there.
Hope this helps,
--
Gregory (Greg) L. Snow Ph.D.
Statis
One option is to use VNC along with vncrec to do the recording (see the
website: http://www.sodan.org/~penny/vncrec/). I think there are some
other recorders also available for vnc, so you might try a google
search.
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthca
Dear R helpers
I would like to know if the "Ranking and Selection" statistical
procedure has been implemented in R. I made a quick search in the R
packages list but I could not find it.
Thanks in advance
Prasanna
__
R-help@stat.math.ethz.ch mailing lis
Dear all,
To Gabor Grothendieck, (again) thanks you very much for your help.
Now, I can play around with lattice package.
Best, Muhammad Subianto
#Gabor
#reduce the data to a frequency matrix and
#then plot it using classic and then lattice graphics:
zm <- as.matrix(rowsum(z1[-9], z1[,9]))
ba
> I have a table with following collumns: State, SamplePlot, Species and
BodySize. I sampled bird species at
> 34 SamplePlots and 5 States (regions) monthly during two years. On each bird
record I measured bodysize
> and identified the species. So I have many records of each species (about 150
sp
Dear R-friends,
I have a table data structured by group, subgroups, records and attributes.
For each group and subgroup I have differente number of records (more than
200). I need bootstrap 100 records for each group/subgroup combinations and
repeat it a big number of times.
Could so
Dear all,
During my pre-R era I tried (yes, tried) to understand mixed models by
working through the 'rat example' in Sokal and Rohlfs Biometry (2000)
3ed p 288-292. The same example was later used by Crawley (2002) in his
Statistical Computing p 363-373 and I have seen the same data being used
On Wed, 30 Aug 2006, Rick Bischoff wrote:
>>> Unfortunately, it seems that most(all?) of R's graphics and summary
>>> statistic functions don't take a weight or frequency argument.
>>> (Fortunately the models do...)
>>
>> I have been been meaning to add this functionality to my graphics
>> package
>> Unfortunately, it seems that most(all?) of R's graphics and summary
>> statistic functions don't take a weight or frequency argument.
>> (Fortunately the models do...)
>
> I have been been meaning to add this functionality to my graphics
> package ggplot (http://had.co.nz/ggplot), but unfortunat
maybe something like this could be of help:
max.score <- c(3,4,3) # max score for each item
all.pats <- as.matrix(expand.grid(lapply(max.score, ":", 1)))
all.pats[rowSums(all.pats) == 5, ]
Best,
Dimitris
Dimitris Rizopoulos
Ph.D. Student
Biostatistical Centre
School of Public Health
Catho
OK,
With this extra detail, a
third solution follows, which may be closer in spirit
to your application.
It may or may not be faster than the other two,
depending on the exact parameters used:
library(partitions)
1+blockparts(n=15,y=c(3,4,2,5,6)-1,include.fewer=T)
(720 distinct solutions)
rks
In each case, look around (help.search,
RSiteSearch) to see if you can find a function
that handles weights. For the case you mention,
medians, it can be done via quantile regression:
x <- w <- 1:5
library(quantreg)
coef(rq(x ~ 1, weight = w))
On 8/30/06, Rick Bischoff <[
"Romain Lorrilliere" <[EMAIL PROTECTED]> wrote in message
news:[EMAIL PROTECTED]
Hi,
do you know, a method to convert an decimal value (integer) to the
corresponding hexadecimal value ?
Starting in R 2.1.0, sprintf can be used:
> x <- c(0, 65535, 65536, 305419896, 2^31-1)
> y <- sprintf("0x%X",
An other way would be:
a <- 123
class(a) <- "hexmode"
a
[1] "7b"
On 30 Aug 2006, at 16:26, mel wrote:
> ?sprintf
>
> ex :
>> sprintf('%X',10)
> [1] "A"
>
> hih
>
> __
> R-help@stat.math.ethz.ch mailing list
> https://stat.ethz.ch/mailman/listinfo/r-hel
Hi Duncan
Here is a bit more detail, this is a bit tough to explain, sorry for not
being clear. Ordering is not important because the vector I am creating
is used as a sufficient statistic in an optimization routine to get some
MLEs. So, any combination of the vector that sums to X is OK. But, the
The data sets I am working with all have a weight variable--e.g.,
each row doesn't mean 1 observation.
With that in mind, nearly all of the graphs and summary statistics
are incorrect for my data, because they don't take into account the
weight.
For example "median" is incorrect, as th
how about learning to use help.search? (In may you already asked a
similiar question...)
help.search("hexadecimal")
which would lead you to format.hexmode
?format.hexmode
Romain Lorrilliere schrieb:
> Hi,
>
> do you know, a method to convert an decimal value (integer) to the
> corresponding
Hi,
Can anyone tell me which function in R gives the hessian matrix of the
log-likelihood of a nonlinear mixed model? fdHess is for scarlar function
only.
Thanks in advance!
Hongmei
[[alternative HTML version deleted]]
__
R-help@stat.math.eth
?sprintf
ex :
> sprintf('%X',10)
[1] "A"
hih
__
R-help@stat.math.ethz.ch 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-contai
use sprintf
> sprintf("%x",123)
[1] "7b"
>
On 9/30/06, Romain Lorrilliere <[EMAIL PROTECTED]> wrote:
> Hi,
>
> do you know, a method to convert an decimal value (integer) to the
> corresponding hexadecimal value ?
>
> thinks for help.
>
> Romain
>
> --
>
> Lorrillière Romain
>
>
>
> UMR 8079 La
Try:
sprintf("%x", 109)
On 9/30/06, Romain Lorrilliere <[EMAIL PROTECTED]> wrote:
> Hi,
>
> do you know, a method to convert an decimal value (integer) to the
> corresponding hexadecimal value ?
>
> thinks for help.
>
> Romain
>
> --
>
> Lorrillière Romain
>
>
>
> UMR 8079 Laboratoire Ecologie,
Try this. First we reduce the data to a frequency matrix and
then plot it using classic and then lattice graphics:
zm <- as.matrix(rowsum(z1[-9], z1[,9]))
barplot(zm, beside = TRUE, col = grey.colors(2))
legend("topleft", legend = levels(z1[,9]), fill = grey.colors(2))
library(lattice)
barchart
I think I've got it now.
If I understand your question, try:
> x <- do.call("cbind",lapply(5:7,function(i){restrictedparts(i,
5,include.zero=FALSE)}))
> acceptable <- function(x){all(x<=c(3,4,5,2,6))}
> x[,apply(x,2,acceptable)]
[,1] [,2] [,3] [,4]
[1,]1232
[2,]1
On 8/30/2006 9:49 AM, Doran, Harold wrote:
> Dear list
>
> Suppose I have the following vector:
>
> x <- c(3,4,2,5,6)
>
> Obviously, this sums to 20. Now, I want to have a second vector, call it
> x2, that sums to x where 5 <= x <= 20, but there are constraints.
>
> 1) The new vector must be s
Hi,
do you know, a method to convert an decimal value (integer) to the
corresponding hexadecimal value ?
thinks for help.
Romain
--
Lorrillière Romain
UMR 8079 Laboratoire Ecologie, Systématique et Evolution
Bât. 362
Université Paris-Sud
91405 Orsay cedex
France
tel : 01 69 15 56
Dear R-list,
I'm trying to estimate the relative importance of 6 environmental variables
in determining clam yield. To estimate clam yield a previous work used the
function Yield = (A^a*B^b*C^c...)^1/(a+b+c+...) where A,B,C... are the
values of the environmental variables and the weights a,b,c... h
Dear Harold
package "partitions" does almost this:
> library(partitions)
> x <- 1+restrictedparts(15,5)
> x[,1:10]
[,1] [,2] [,3] [,4] [,5] [,6] [,7]
[1,] 16 15 14 13 12 11 10
[2,]1234567
[3,]1111111
[4,]11
Dear list
Suppose I have the following vector:
x <- c(3,4,2,5,6)
Obviously, this sums to 20. Now, I want to have a second vector, call it
x2, that sums to x where 5 <= x <= 20, but there are constraints.
1) The new vector must be same length as x
2) No element of the new vector can be 0
3) Ele
Andras,
At this point you need to write your own function to take the
posterior density sample (stored in a coda mcmc object) and
covariates of interest to get a Monte Carlo estimate of these
probabilities.
Best,
ADM
On Aug 30, 2006, at 5:10 AM, Andras Treszl wrote:
> Hi,
>
> I am using
Chris Barker ucdavis.edu> writes:
>
However, after reading the packages' documentation and
> searching the mailing lists and other sources, it seems that the
> functions in geoRglm do not acknowledge the existence of multiple
> realisations per site. I see that the likfit() function in geoR
"COMTE Guillaume" <[EMAIL PROTECTED]> writes:
> Hi all,
>
> I'm getting confused by date handling.
>
>
>
> I wish to read a date from a file wich is a number of seconds since 1970
> (POSIXct).
>
> Then i wish to convert this date to a human readable form (POSIXlt)
>
>
>
> By example :
>
Uwe Ligges wrote:
>
>
> Rainer M Krug wrote:
SNIP
>> My scenario: I simulated densities of a plants originating from an sseed
>> source at distance zero. Therefore the density of the plants will be
>> highest close to zero. Is there anything I can do to have this pattern?
>> If I use 'from' or
Hi all,
I'm getting confused by date handling.
I wish to read a date from a file wich is a number of seconds since 1970
(POSIXct).
Then i wish to convert this date to a human readable form (POSIXlt)
By example :
ctDate<-1132963200 #"Wed Aug 30 14:24:37 2006"
is(ctDate)
[1] "numeric"
I've submitted a new, and substantially enhanced, version (1.2-0) of the
Rcmdr package to CRAN. Some highlights (from the CHANGES) file:
o Added ability to import from Excel, Access or dBase files (contributed
by Samir Messad, Renaud Lancelot and Matthieu Lesnoff).
o Added ability to read
Dear all,
Many Thanks to Jacques VESLOT and Jim Lemon for their helps.
Best, Muhammad Subianto
#Jacques VESLOT
barplot(t(sapply(split(z1[,1:8], z1$V9),colSums)), beside=T)
#Jim Lemon
barplot(sapply(z1[1:8],by,z1[9],sum),beside=TRUE)
On this day 30/08/2006 11:43, Muhammad Subianto wrote:
> D
On 08/30/06 11:55, Romain Francois wrote:
> Robert Mcfadden a ?crit :
> > May be it's not a bug, but I tried to search for the package rpanel and I
> > was not find. At the r-project's site that package is available.
> > How to explain it?
> > Rob
> >
> Hi Rob,
>
> If it's not there : http://finzi
Romain Francois wrote:
> Robert Mcfadden a écrit :
>> May be it's not a bug, but I tried to search for the package rpanel and I
>> was not find. At the r-project's site that package is available.
>> How to explain it?
>> Rob
>>
> Hi Rob,
>
> If it's not there : http://finzi.psych.upenn.edu
Rainer M Krug wrote:
> Hi
>
> the function density() does normally integrate to one - I've checked it
> and it works and I also read the previous threads.
> But I realised that it does not integrate to one if I use from, to or cut.
>
> My scenario: I simulated densities of a plants originating
stat stat wrote:
> Dear all R users,
>
> Suppose,
>
>
> Dear all R users,
>
> Suppose,
>
>
> pauto.cor = pacf(lh, plot=F)
> max.lag = max(pauto.cor$lag)
> min.lag = min(pauto.cor$lag)
> centre = (max.lag - min.lag)/2
> pauto.cor = pauto.cor$acf
> pauto.cor = pauto.cor[-1]
>
> par(mar=c(3,0
Robert Mcfadden a écrit :
> May be it's not a bug, but I tried to search for the package rpanel and I
> was not find. At the r-project's site that package is available.
> How to explain it?
> Rob
>
Hi Rob,
If it's not there : http://finzi.psych.upenn.edu/R/library/
it's not on the extension
May be it's not a bug, but I tried to search for the package rpanel and I
was not find. At the r-project's site that package is available.
How to explain it?
Rob
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE
Muhammad Subianto wrote:
> ...
> I have tried to learn ?xtabs ?table and ?ftable but I can't figure out.
> I need a barplot for all variables and the result maybe like
>
> | | | |
> | | | | | || | |
> |pos|neg| |pos|neg||pos|neg
On 8/30/06, Mark Lyman <[EMAIL PROTECTED]> wrote:
> > I have a problem with apply function. I have to two matrices of dimension of
> > one column but n rows. I have to check whether one matrix is greater than
> > other by going thru each row (ie) using if condition to check one matrix
> > with anot
Hi
the function density() does normally integrate to one - I've checked it
and it works and I also read the previous threads.
But I realised that it does not integrate to one if I use from, to or cut.
My scenario: I simulated densities of a plants originating from an sseed
source at distance zero
Hi,
I am using MCMCpack and the MCMClogit function to create logistic
regression models in a medical (adverse event) study. My question is,
is there a way where I can directly create the estimated probabilities
of the adverse outcome, and the confidence interval for the
estimated probabilities? Or
barplot(t(sapply(split(z1[,1:8], z1$V9),colSums)), beside=T)
---
Jacques VESLOT
CNRS UMR 8090
I.B.L (2ème étage)
1 rue du Professeur Calmette
B.P. 245
59019 Lille Cedex
Tel : 33 (0)3.20.87.10.44
Fax : 33 (0)3.20.87.10.31
http://www-
That is not an R package.
It looks like the sources of the server described on
http://cran.r-project.org/contrib/extra/dcom/RSrv200.html
which tells you about the approriate mailing list to ask about it.
You probably want to use RSrv200.exe to install the server (and other
tools and examples).
Dear all,
I have a dataset. I want to make barplot from this data.
Zero1 <- "
V1 V2 V3 V4 V5 V6 V7 V8 V9
1 1 0 0 0 1 0 0 0 Positive
2 0 0 1 0 1 0 1 1 Negative
3 0 0 1 0 0 0 1 1 Positive
4 0 1 0 1 1 1 0 1 Negative
5 0 0 1 0 1 1 0 0 Positive
6 0
Dear all R users,
Suppose,
Dear all R users,
Suppose,
pauto.cor = pacf(lh, plot=F)
max.lag = max(pauto.cor$lag)
min.lag = min(pauto.cor$lag)
centre = (max.lag - min.lag)/2
pauto.cor = pauto.cor$acf
pauto.cor = pauto.cor[-1]
par(mar=c(3,0,1,1))
barplot(pauto.cor, axes=F,xlim=c(max(pauto.cor),
Hi,
I am using MCMCpack and the MCMClogit function to create logistic
regression models in a medical (adverse event) study. My question is,
is there a way where I can directly create the estimated probabilities
of the adverse outcome, moreover the confidence interval for the
estimated probabilitie
If you are using 'only' English then
str <- "dog"
strsplit(str,NULL)[[1]]
works perfectly and it is fast.
But if you also dealing with Unicode character have a look at
http://wiki.r-project.org/rwiki/doku.php?id=tips:data-
strings:decomposestring
Cheers,
Hans
> you can also use substring(
I am trying to install SrcStatConnectorSrv(2) and rcom from local zip
files.
I have successfully downloaded both files (from CRAN for rcom CRAN Other
for SrcStatConnectorSrv) and installed rcom. However, I get the following
error message when I try to install SrcStatConnectorSrv (the version date
Respected Sir/Madam,
I have a problem with apply function. I have to two matrices of dimension of
one column but n rows. I have to check whether one matrix is greater than
other by going thru each row (ie) using if condition to check one matrix
with another matrix.
I like to use apply() function
you can also use substring(), e.g.,
substring(x3, 1:nchar(x3), 1:nchar(x3))
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: +32/(0)16/337015
82 matches
Mail list logo
