Dear all
as usual I am again lost in virtues of regular expressions.
I have such character vector named vzor:
[365] "61A" "62C/27" "65A/27" "66C/29" "69A/29" "70C/31"
"73A/31" "74C/33" "77A/33" "81A/35" "82C/37" "85A/37" "86C/39"
[378] "89A/39" "90C/41" "93A/41"
gsub("^.*([[:alpha:]]).*$", "\\1", vzor)
..<°}))><
) ) ) ) )
( ( ( ( (Prof. Philippe Grosjean
) ) ) ) )
( ( ( ( (Numerical Ecology of Aquatic Systems
) ) ) ) ) Mons-Hainaut University, Belgium
( ( ( ( (
..
On Wed, 18 Apr 2007, Atte Tenkanen wrote:
> I often get a following error with R
>
> *** caught segfault ***
> address 0x78807e00, cause 'memory not mapped'
>
> Possible actions:
> 1: abort (with core dump)
> 2: normal R exit
> 3: exit R without saving workspace
> 4: exit R saving workspace
>>
> S
On Wed, 18 Apr 2007, Ulrik Stervbo wrote:
> As I har problems installing the Cairo package, I went for Henriks solution
> - and it works almost perfect. I would like to have been able to generate
> transparent png.
You cannot do transparency via postscript.
I would suggest using pdf() and conver
> What is the optimum range to look for a value of lambda while doing ridge
> regression. Can/ should lambda be greater than 1 ?
-- I think it's data dependent, but lambda can certainly be greater than one.
For many ridge regression problems you can choose lambda `objectively' by
generalized cro
Hi all,
I used to connect internet via a proxy. Before update packages I wrote in R
Sys.putenv("http_proxy"="http://proxy3.redegov.sp.gov.br:80/";)
Nevertheless the way the connection is done has changed. For example, in the
browser the proxy is not indicated and I have to give an username and
Hi,
Help me, how I can Install made4(micoarray analysis tool) in R using linux
OS.
thanking you.
--
Nitish Kumar Mishra
Junior Research Fellow
BIC, IMTECH, Chandigarh, India
E-Mail Address:
[EMAIL PROTECTED]
[EMAIL PROTECTED]
__
R-help@stat.math.ethz.
Another way you can do it, if the data has the pattern shown in your
sample, it to select all the lines that start with a numeric:
> input <- " FILE-CONTENT ##
+ EXAM NUM:2
+ -
+ EXAM #1
+ ASTIG:-2.4D
+ AXIS:4.8
+ START OF HEIGHT DATA
+ 0 0.0 0.
+ 0 0.
On Tue, 17 Apr 2007, Antonio Olinto wrote:
> Hi all,
>
> I used to connect internet via a proxy. Before update packages I wrote in R
>
> Sys.putenv("http_proxy"="http://proxy3.redegov.sp.gov.br:80/";)
>
> Nevertheless the way the connection is done has changed. For example, in the
> browser the pr
Sorry, but this works under all the circumstances I tried on my Vista
system, so there is nothing I can do to debug it.
On Tue, 17 Apr 2007, Prof Brian Ripley wrote:
I suspect tcl's own version of 'access', but can you please confirm that this
still happens under 'Run as Administrator', assumi
> On Wed, 18 Apr 2007, Atte Tenkanen wrote:
>
> > I often get a following error with R
> >
> > *** caught segfault ***
> > address 0x78807e00, cause 'memory not mapped'
> >
> > Possible actions:
> > 1: abort (with core dump)
> > 2: normal R exit
> > 3: exit R without saving workspace
> > 4: exit
A backreference is contained in parentheses and there are no parentheses
in your regular expression, hence the error message.
Its probably easiest just to remove all non-letters:
x <- "45x53yy66"
gsub("[^[:alpha:]]", "", x) # xyy
On 4/18/07, Petr PIKAL <[EMAIL PROTECTED]> wrote:
> Dear all
>
> a
Now problems with plot-command...
I try to plot Julia set and the algorithm works, if the points are drawn during
the loop. But if I want to save the values first to the matrix and then
afterwards plot them at once, the picture is distorted. What's wrong with the
plot-command? I think the Point
Thank you all for your working solutions
Petr Pikal
[EMAIL PROTECTED]
[EMAIL PROTECTED] napsal dne 18.04.2007 13:26:36:
> A backreference is contained in parentheses and there are no parentheses
> in your regular expression, hence the error message.
>
> Its probably easiest just to remove all
Thomas Lumley wrote:
> On Tue, 17 Apr 2007, Frank E Harrell Jr wrote:
>
>> The points that Thomas and Brian have made are certainly correct, if
>> one is truly interested in testing for differences in medians or
>> means. But the Wilcoxon test provides a valid test of x > y more
>> generally.
Dear group members,
I want to compare response variables ("logAUC") of two groups (treatment
"Test", "Reference") of a subset ("period == 1") in dataframe "resp"
(below):
sequence subject period treatment AUC logAUC
1RT 1 1 Reference 44.1 3.786460
2RT 1
On Apr 18, 2007, at 8:46 AM, Helmut Schütz wrote:
> Dear group members,
>
> I want to compare response variables ("logAUC") of two groups
> (treatment
> "Test", "Reference") of a subset ("period == 1") in dataframe "resp"
> (below):
[ snip ]
> The formula method of t.test
>
>> result <- t.test
On 4/18/07, Helmut Schütz <[EMAIL PROTECTED]> wrote:
> Dear group members,
>
> I want to compare response variables ("logAUC") of two groups (treatment
> "Test", "Reference") of a subset ("period == 1") in dataframe "resp"
> (below):
>
>sequence subject period treatment AUC logAUC
> 1
Dear list
I'm plotting ( boxplot() and plot() ) some data for a publication.
The editor would like the text labels on the plots in a larger font.
I'm doing something like this:
jpeg(
filename = "D:/Martin/Work/CleanPath/RAF1%03d.jpg",
width = 1000,
height = 600,
take a look at
?relevel()
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
Web: http://med.kuleuven.be/biostat/
http://w
I updated R to the last 2.4.1 version and unfortunately I can not
load languageR any longer.
In R-2.4.1, LanguageR requires acepack, but Hmisc doesn't work when
acepack is loaded.
> library(languageR)
Loading required package: Design
Loading required package: Hmisc
take a look at ?par() and especially at the 'cex.lab', 'cex.main', and
'cex.axis' arguments. For instance, you could use something like:
plot(
dafExpo$grn07,dafExpo$grn09,
xlim=c(0,numAxeMax),
ylim=c(0,numAxeMax),
pch=1,
xlab="Green at 7 AM : NOx [(µg/m3)*hour]",
ylab="Green at 9 AM : NOx [(µg/m3
Prof Brian Ripley wrote:
> Sorry, but this works under all the circumstances I tried on my Vista
> system, so there is nothing I can do to debug it.
>
You (i.e. Sofia) could do some investigation yourself. It may prove
informative if you search for init.tcl and check whether it is readable
(for
Dear Charilaos!
Charilaos Skiadas wrote:
>>
>> Do you know a more elegant way than the clumsy one I have tried?
>>
>>> as.numeric(exp(result$estimate[2]-result$estimate[1]))
>>> as.numeric(exp(-result$conf.int[2]))
>>> as.numeric(exp(-result$conf.int[1]))
>
> First off, those three could probably
Prof Brian Ripley wrote:
> On Wed, 18 Apr 2007, Ulrik Stervbo wrote:
>
>> As I har problems installing the Cairo package, I went for Henriks solution
>> - and it works almost perfect. I would like to have been able to generate
>> transparent png.
>
> You cannot do transparency via postscript.
>
Dear R-Experts,
how can I divide the number 0.285 with 2. I need a function.
Result: 0.285 / 2 = 0.1425
Thanks, Corinna
__
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-
Schmitt, Corinna wrote:
> Dear R-Experts,
>
> how can I divide the number 0.285 with 2. I need a function.
> Result: 0.285 / 2 = 0.1425
Just get the / operator:
> divide = get("/")
>
> divide(0.285,2)
[1] 0.1425
Is that what you want?
Barry
__
Thanks for the help! I found out that it was a problem with the search path
for R. Reinstalling R seems to have solved the problem. Sorry for bothering
you before having tested that.
Sofia
-Ursprungligt meddelande-
Från: Peter Dalgaard [mailto:[EMAIL PROTECTED]
Skickat: den 18 april 2007
Well, i think
half.of.0.285 <- function() {
0.1425
}
would do the trink.
Gabor
On Wed, Apr 18, 2007 at 04:42:49PM +0200, Schmitt, Corinna wrote:
> Dear R-Experts,
>
> how can I divide the number 0.285 with 2. I need a function.
> Result: 0.285 / 2 = 0.1425
>
> Thanks, Corinna
>
> _
Dear R-experts,
It is a quite stupid question but please help me. I am very confuced. I
am able to import normal txt ant mat-files to R but unable to import
.xls-file
I do not understand the online help. Can please anyone send me the
corresponding command lines? The .xls-file is attached. In my f
Hi,
is there a way of telling Emacs + ESS to show words that are already
a function in R (such as 'length') is a different colour/font?
Best,
Federico
--
Federico C. F. Calboli
Department of Epidemiology and Public Health
Imperial College, St. Mary's Campus
Norfolk Place, London W2 1PG
Tel +
Hello,
I have got a numeric matrix with three colums of a few measurments. My x and
y coordinates are rising numbers. The z coordinate is the measurment. In
my matrix if have got ~6000 sorted values for one measurments multiplied
with the number of measurments. An example in the end.
My aim is to
Has anyone proposed using a bootstrap for Pedro's problem?
What about taking a boostrap sample from x, a boostrap sample from y, take
the difference in the medians for these two bootstrap samples, repeat the
process 1,000 times and calculate the 95th percentiles of the 1,000
computed differences?
Dear all R gurus,
I have following syntax:
y = c(1:10)
chippy <- function(x)
{
y[5] = x
sin(cos(t(y)%*%y)*exp(-t(y)%*%y/2))
}
curve(chippy, 1, 20, n=200)
But I am getting error while executing :
Error in xy.coords(x, y, xlabel, ylabel, log) :
Hi all,
I am trying to fit an ANOVA model in R using the aov/lm commands. I have a
set of observational (i.e. no fixed experimental effects) data, in which I
have identified high and low clusters of the response variable. The design
is unbalanced, with 773 high cluster observations, and 523 low cl
Your chipply function is not vectorized. See ?curve and try:
curve(Vectorize(chippy)(x), 1, 20, n=200)
On 4/18/07, Ron Michael <[EMAIL PROTECTED]> wrote:
>
> Dear all R gurus,
>
>
>
> I have following syntax:
>
>
>
> y = c(1:10)
> chippy <- function(x)
>{
> y[5] = x
> sin
[EMAIL PROTECTED] wrote:
> Has anyone proposed using a bootstrap for Pedro's problem?
>
> What about taking a boostrap sample from x, a boostrap sample from y, take
> the difference in the medians for these two bootstrap samples, repeat the
> process 1,000 times and calculate the 95th percentiles
Corinna Schmitt wrote:
>
> It is a quite stupid question but please help me. I am very
> confuced. I am able to import normal txt ant mat-files to R but
> unable to import .xls-file
>
I've tried two ways to import excel files, but none of them
seems perfect.
Method 1:
This method uses library
2007/4/18, Schmitt, Corinna <[EMAIL PROTECTED]>:
> It is a quite stupid question but please help me. I am very confuced. I
> am able to import normal txt ant mat-files to R but unable to import
> .xls-file
Searching for Excel on e.g. http://www.r-project.org/search.html,
http://tolstoy.newcastle.e
There is also a read.xls command in package gdata, it seems that it uses
a perl script called 'xls2csv'. I've have no idea how good this is,
never tried it.
Btw, xlsReadWrite is Windows-only, so you can use it only if
you use windows.
Gabor
ps. Corinna, to be honest, i've no idea what kind onli
Hi R users,
Does anybody knows for the following erro after running
geweke.diag(MCMC.sampled, frac1=0.1, frac2=0.5)
Erro em glm.fit(x = X, y = Y, weights = weights, start = start, etastart =
etastart, :
laço interno 1; não é possÃvel corrigir o tamanho do passo
Além disso: W
I use gdata and it works quite well for me. It's as easy as
install.packages(gdata)
library(gdata)
data = read.xls("mydata.xls",sheet=1)
[read.xls() can take other arguments]
It requires concurrent installation of Perl, but installing Perl is also
simple. For Windows, you can get it here:
http:
To avoid complications, save your file as comma separated and use one
of the instructions for reading delimited files. If you are using a
comma as a decimal point you are probably using ; as a separator. If
this is so use read.csv2. Please see the help files for read.table.
Best Regards
John
...is this what you're looking for?
donedat <- subset(data,ID < 6000 | ID >= 7000)
findat <- donedat[-unique(rapply(donedat,function(x)
which( x < 0 ))),,drop=FALSE]
the second line looks through each column, and finds the indices of negative
values - rapply() ret
Gabor Csardi wrote:
>
> There is also a read.xls command in package gdata, it seems that it uses
> a perl script called 'xls2csv'. I've have no idea how good this is,
> never tried it.
>
> Btw, xlsReadWrite is Windows-only, so you can use it only if
> you use windows.
>
Ok, but who would be insa
John C Frain wrote:
>
> To avoid complications, save your file as comma separated and use one
> of the instructions for reading delimited files. If you are using a
> comma as a decimal point you are probably using ; as a separator. If
> this is so use read.csv2. Please see the help files for rea
Hello Everybody,
Install the package: install.packages("xlsReadWrite")
Load it: library(xlsReadWrite)
testfile = read.xls("TesFile.xls")
Have Fun!
Kind Regards,
Soare Marcian-Alin
PS: If dont works, then install also the package "xtable", but it should
work without installing it!
2007/4/18, J
On Wed, Apr 18, 2007 at 03:51:35PM -0200, Alberto Monteiro wrote:
> Gabor Csardi wrote:
> >
> > There is also a read.xls command in package gdata, it seems that it uses
> > a perl script called 'xls2csv'. I've have no idea how good this is,
> > never tried it.
> >
> > Btw, xlsReadWrite is Windows-
For testing, the permutation test may be prefered to the bootstrap
(though the bootstrap could be used for a confidence interval).
I remember in grad school doing a project on comparing the efficiency of
a permutation test on medians compared to the MannWhitney test, but I
don't remember the speci
On Wed, 18 Apr 2007, Greg Snow wrote:
> For testing, the permutation test may be prefered to the bootstrap
> (though the bootstrap could be used for a confidence interval).
>
> I remember in grad school doing a project on comparing the efficiency of
> a permutation test on medians compared to the
Dear all,
Sorry to bother you, but I didn't find the solution in the R-help archive.
I would like to change the thickness of stripes inside the barplot. Is there any
solution
to do that when using the barplot() fuction and the density option ?
Or is there any other function ?
Thanks very much fo
same problem here.
last time I had a similar one when I did "library(MASS)", I solved
that by re-installation of R 2.4.1. However, this time it does not
work.
On 4/18/07, Lara Tagliapietra <[EMAIL PROTECTED]> wrote:
> I updated R to the last 2.4.1 version and unfortunately I can not
> load langua
Here is a workable version for the Julia set. I put it also to
http://fractalswithr.blogspot.com/
Atte
> Now problems with plot-command...
>
> I try to plot Julia set and the algorithm works, if the points are
> drawn during the loop. But if I want to save the values first to
> the matrix a
Hello,
Same problem on a MacBook Pro (intel) with RODBC, nortest and gplots.
Best regards
Le 18-avr.-07 à 21:09, Weiwei Shi a écrit :
> same problem here.
>
> last time I had a similar one when I did "library(MASS)", I solved
> that by re-installation of R 2.4.1. However, this time it does not
Dear All:
Pleas help me to increase the memory in R.
I am trying to make euclidean distance matrix.
The number of low in data is 500,000. Therefore, the dimension of euclidean
distance matrix is 500,000*500,000.
When I run the data in R. R could not make distance matrix because of memory
al
You would need 2TB (2,000,000,000,000) to store a single copy of your
data. You probably need to rescale your problem. Even if you had the
memory, the computation would take a very long time.
On 4/18/07, Hong Su An <[EMAIL PROTECTED]> wrote:
> Dear All:
>
> Pleas help me to increase the memory i
In their paper, "Lexical Scope and Statistical Computing", the authors (
Gentleman and Ihaka ) go to great length explaining why R's use of
lexical scoping creates advantages when doing statistical computations.
If anyone has or is familiar with this paper, could they provide the
main program code
Hi Vadim,
Estimates of the "alpha_0" terms in MASS are the $yb component
of the object returned by ppr(). As I understand it,
the original PPR algorithm assumes the response
variable(s) are centered, so the 'alpha_0' term in
MASS is just the mean of the response if the user
does not center the r
Hi I just installed the gmodels package
and the installation was successful but when I was trying to load the
library I got an error (see below).
Interesting, yesterday I wrote to the maintainer of RSQLite apckage
because I got the same error.
Does somebody knows what is going on ??
thanks,
Ma
If I print a sting I get an initial [1]:
> xx="a"
> xx
[1] "a"
How do I get it to print just
"a"
with no [1]?
I tried looking this up, but I don't know what the initial [1] is called.
Steve
__
R-help@stat.math.ethz.ch mailing list
https://stat.e
> x="a\n"
> cat(x)
a
On Apr 18, 2007, at 5:31 PM, steve wrote:
> If I print a sting I get an initial [1]:
>
>> xx="a"
>> xx
> [1] "a"
>
>
> How do I get it to print just
>
> "a"
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/li
-- Forwarded message --
From: John C Frain <[EMAIL PROTECTED]>
Date: 18-Apr-2007 22:35
Subject: Re: [R] importing excel-file
To: Alberto Monteiro <[EMAIL PROTECTED]>
One additional suggestion would be to use gretl. Gretl will read
excel files with an option to possibly ignore the
Hello Steve,
You can print strings in R with the method cat()
\n . new line
\t . tabulator
Try:
name <- c("Steve")
age=22
cat("\tHello my name is", name ,"and I am", age ,"years old.\n")
Have Fun!
Kind Regards,
Soare Marcian-Alin
2007/4/18, steve <[EMAIL PROTECTED]>:
>
> If I print a
is there no package/function in R to calculate the conditional power or the
bayesian predictive power for trials with binary endpoints? Thanks
--
View this message in context:
http://www.nabble.com/Conditional-power%2C-predictive-power-tf3603396.html#a10066991
Sent from the R help mailing list a
As part of carrying out a complicated maximum likelihood estimation, I
am trying to learn to program likelihoods in R. I started with a simple
probit model but am unable to get the code to work. Any help or
suggestions are most welcome. I give my code below:
ml
not sure just what you want, but here are some snippets
newton <-
function(lfun, est, tol = 1e-7, niter = 500) {
cscore <- lfun$score(est)
if (abs(cscore) < tol)
return(est)
for (i in 1:niter) {
new <- est - cscore / lfun$d2(est)
cscore <- lfun$score(new)
On Apr 18, 2007, at 1:58PM , Gabor Csardi wrote:
> On Wed, Apr 18, 2007 at 03:51:35PM -0200, Alberto Monteiro wrote:
>> Gabor Csardi wrote:
>> Ok, but who would be insane enough to use Excel in Linux or Mac? :-)
The original reason for read.xls() was to allow a web application
running computat
Dear Happy R-users & experts,
I am in need of advice,
While working with spatial data (x & y coordinates of seed locations) I have
come accross the problem that I need to convert my point data into a matrix or
grid system. I then need to count how often a point falls into a certain
position in
I have received a number of reports of problems with recent unversal
Mac packages from CRAN when used with R 2.4.1. Has something in the
build script changed?
-G
On Apr 18, 2007, at 4:49PM , Mayte Suarez-Farinas wrote:
> Hi I just installed the gmodels package
> and the installation was su
If I have a data frame X that looks like this:
A B
- -
1 2
1 3
1 4
2 3
2 1
2 1
3 2
3 1
3 3
and I want to make another column which has the rank of B computed
separately for each value of A.
I.e. something like:
A B C
- - -
1 2 1
1 3 2
1 4 3
2 3 3
2 1 1
2 1 2
3 2 2
3 1 1
3 3 3
by(X, X[,1], func
The issue Greg mentions, that "most scientists store their
experimental data in MS-Excel spreadsheets"
is the motivation for one of the sessions at the Interface
2007 conference
http://sbm.temple.edu/interface07/index.html
in Philadelphia, May 23-26, 2007
Erich and Thomas designed the RExcel int
On Wed, 18 Apr 2007, Marco Visser wrote:
> Dear Happy R-users & experts,
>
> I am in need of advice,
> While working with spatial data (x & y coordinates of seed locations) I have
> come accross the problem that I need to convert my point data into a matrix
> or grid system. I then need to count
Consider the code:
x <- seq(0,1,0.2)
y <- seq(0,1,0.01)
cbind(match(y,x),y)
which, surprisingly, doesn't show a match at 0.6! (It gives correct
matches at 0, 0.2, 0.4, 0.8 and 1, though)
In addition,
x[4]==y[61]
yields FALSE. (but x[5]==y[81], the one for 0.8, yields TRUE)
Is this a consequ
On 4/18/07, Bernhard Klingenberg <[EMAIL PROTECTED]> wrote:
> Consider the code:
>
> x <- seq(0,1,0.2)
> y <- seq(0,1,0.01)
> cbind(match(y,x),y)
>
> which, surprisingly, doesn't show a match at 0.6! (It gives correct
> matches at 0, 0.2, 0.4, 0.8 and 1, though)
>
> In addition,
>
> x[4]==y[61]
>
>
Deepankar,
Some general advice from a non-expert:
Write your likelihoods without a for loop. This is important because
the likelihood is evaluated multiple times in the maximization process
and you don't want to be looping looping looping ...
Always try multiple starting values
Sometimes it hel
Does this do what you want?
> x <- "A B
+ 1 2
+ 1 3
+ 1 4
+ 2 3
+ 2 1
+ 2 1
+ 3 2
+ 3 1
+ 3 3"
> x <- read.table(textConnection(x), header=TRUE)
> x$C <- ave(x$B, x$A, FUN=rank)
> x
A B C
1 1 2 1.0
2 1 3 2.0
3 1 4 3.0
4 2 3 3.0
5 2 1 1.5
6 2 1 1.5
7 3 2 2.0
8 3 1 1.0
9 3 3 3.0
On 4/18/07, Lu
Dear R-helpers,
I am looking for a segmentation package that gives some "tree identifier"
as output for every observation in the data set (my response variable is
binary). I have skimmed through "rpart", "ada" and "adabag": The output
"trees" gives you the formula, but I have to run several tho
Hello I am using the "for (i...)" and a sink() into a file. But the output I
am having is not arranged in either a vector or any other good structure. I
would like to have the output in a file directly as a vector so that I do
not have to edit the [1] and [6] etc and that the values are comma
sepa
78 matches
Mail list logo