On Fri, Nov 21, 2008 at 1:27 AM, A Ezhil <[EMAIL PROTECTED]> wrote:
> I have a matrix of size 1 x 50. I would like to calculate all possible
> pair-wise correlation coefficient (5x10^7 combinations) using cor(). How can
> I efficiently calcualte and save the result in a matrix?
>
You might al
I am not sure what the issue is here. Do you want to capture both stderr
and stdout (use 2>1 in the command with an sh-like shell), or is the
problem that you don't get immediate output?
The latter is a Perl issue: you need to circumvent output buffering.
See e.g
http://perl.plover.com/FAQs/Buf
Ben Zuckerberg schrieb:
Greetings,
I have several sets of oscillation data and would like to estimate the
parameters of a sine function to each set (and hopefully automate
this). A colleague provided an excel sheet that uses solver to minimize
the RSS after fitting the sine function to each
Hi,
I have an application in perl that prints some output to either stderr
or stdout.
Here's an example:
# tmp.pl
print STDERR "starting iterator\n";
for(my $i = 0; $i < 100; $i++) {
print $i . "\n";
}
# tmp.R
con <- pipe("perl tmp.pl")
r <- readLines(con, n = -1)
close(con)
However,
or rbinom(n, size = 1, prob)
Regards,
Yihui
--
Yihui Xie <[EMAIL PROTECTED]>
Phone: +86-(0)10-82509086 Fax: +86-(0)10-82509086
Mobile: +86-15810805877
Homepage: http://www.yihui.name
School of Statistics, Room 1037, Mingde Main Building,
Renmin University of China, Beijing, 100872, China
On Thu
Dear James,
Try this:
as.Date('20081121',"%Y%m%d")
[1] "2008-11-21"
HTH,
Jorge
On Thu, Nov 20, 2008 at 9:59 PM, Smith, James (EWLS) <
[EMAIL PROTECTED]> wrote:
> G'day all,
>
>
>
> I have many, very long time series' that come in to R in txt format and
> so the date column is structured:
>
>
Dear list,
First off, let me offer my apologies, I know this is a very basic question.
After amassing a large number of objects (from multiple projects) in one
working directory, I'd like to be able to start using different directories,
just for the sake of organization alone. But I have no i
G'day all,
I have many, very long time series' that come in to R in txt format and
so the date column is structured:
20081121
I would like them all to look like 2008-11-21 but cannot find any way to
add these characters to all dates in one column in R.
I have way too many files to do it man
Dear Anup,
Try this:
# Data
A <- c(1,2,3,4,5)
B <- rnorm(100)
# Results
t(apply(sapply(A,function(x) pnorm(x*B)),2,function(x)
c(Mean=mean(x),Var=var(x
HTH,
Jorge
On Thu, Nov 20, 2008 at 10:09 PM, Anup Menon <[EMAIL PROTECTED]>wrote:
> Dear Friends,
>
> I'm trying to see if there is so
Dear Friends,
I'm trying to see if there is some possibility that I can do the following
computations without a loop. I have attached a toy example below.
A <- c(1,2,3,4,5)
B <- rnorm(100)
store <- matrix(0,5,2)
for (i in 1:5)
{
store[i,1] <- mean(pnorm(A[i]*B))
store[i,2] <- var(pnorm(A[i]*
On Thu, 20 Nov 2008, A Ezhil wrote:
Dear All,
I have a matrix of size 1 x 50. I would like to calculate all
possible pair-wise correlation coefficient (5x10^7 combinations) using
cor(). How can I efficiently calcualte and save the result in a matrix?
By using cor()??
res <- c
Dear All,
I have a matrix of size 1 x 50. I would like to calculate all possible
pair-wise correlation coefficient (5x10^7 combinations) using cor(). How can I
efficiently calcualte and save the result in a matrix?
Thanks in advance.
Kind regards,
Ezhil
__
Hello,
Does anyone know whether a Windows version of the Ratings package will
be available?
Thank you, Paul Prew
CONFIDENTIALITY NOTICE: \ This e-mail communication an...{{dropped:11}}
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailm
yet another way:
> mapply(function(x,y)length(intersect(x,y)), q, w)
[1] 0 1 3 3
On Thu, Nov 20, 2008 at 7:15 PM, Stefan Th. Gries <[EMAIL PROTECTED]> wrote:
> Hi all
>
> I have two lists that have the same number of numeric vectors such as:
> q<-list(1, 2:3, 4:6, 7:10)
> w<-list(0, 1:2, 3:7, 8:
On 21/11/2008, at 10:13 AM, Steffy, Elizabeth A. wrote:
I got this error for this equation and i'm not sure what it means
or how to fix it:
Error in S[index] = S[index - 1] + (dSi - dSo - SC) * dt :
nothing to replace with
Does anyone know how to fix this?
No.
PLEASE do read th
try this:
> q<-list(1, 2:3, 4:6, 7:10)
> w<-list(0, 1:2, 3:7, 8:10)
> result <- sapply(seq(length(q)), function(.indx) length(intersect(q[[.indx]],
> w[[.indx]])))
> result
[1] 0 1 3 3
>
On Thu, Nov 20, 2008 at 7:15 PM, Stefan Th. Gries <[EMAIL PROTECTED]> wrote:
> Hi all
>
> I have two lists t
Reading in between the lines a little, maybe you want
lm(..., na.action = na.exclude)
That should return missing values for the influence statistics when
the predictor or responses is missing in the input.
Hadley
On Thu, Nov 20, 2008 at 4:19 PM, David Kaplan
<[EMAIL PROTECTED]> wrote:
> Let me
Hey,
I'm very curious how to let R identify the procedures which may get
performance improvement using GPU. While *apply and bootstrap
functions may be obviously recognizable, other tasks may not be so
obvious, I guess.
On Wed, Nov 19, 2008 at 11:30 AM, Peter Dalgaard
<[EMAIL PROTECTED]> wrote:
>
Hi David,
I inadvertently introduced a bug in ggplot in the last release. I
uploaded a fix to CRAN this morning and it should be available in the
near future. Sorry for the inconvenience.
Regards,
Hadley
On Thu, Nov 20, 2008 at 2:49 PM, David Hajage <[EMAIL PROTECTED]> wrote:
> Hello R users,
Dear All,
I have a matrix of size 1 x 50. I would like to calculate all possible
pair-wise correlation coefficient (5x10^7 combinations) using cor(). How can I
efficiently calcualte and save the result in a matrix?
Thanks in advance.
Kind regards,
Ezhil
__
Dear Niels
I am trying to construct some anthropometric standards for adults using LMS
and GAMLSS in R. I saw that you had a problem that I have with the package
lmsqreg of the package being reported as "not a valid package". I presume
you solved your problem, but I did not see how in the r-help l
Hello R users,
I have an error with package ggplot2 under linux (ubuntu 8.10), R 2.8.0 and
ggplot2 0.7, everything up to date :
> library(ggplot2)
Le chargement a nécessité le package : grid
Le chargement a nécessité le package : reshape
Le chargement a nécessité le package : plyr
Le chargement a
I got this error for this equation and i'm not sure what it means or how to fix
it:
Error in S[index] = S[index - 1] + (dSi - dSo - SC) * dt :
nothing to replace with
Does anyone know how to fix this?
__
R-help@r-project.org mailing list
https://sta
Hi all
I have two lists that have the same number of numeric vectors such as:
q<-list(1, 2:3, 4:6, 7:10)
w<-list(0, 1:2, 3:7, 8:10)
What I want to do is create a vector desired.result that looks like
this but I am thinking there must be some kind of non-loop /
"\\wapply" way to so ...
desired.re
Dear all,
A am trying to fit a generalized linear mixed effects model with a binomial
link function, my response data is binary, using the lme4 R package, for the
glmer model but with the cauchit link function (CDF of Cauchy distribution),
under the package this has not yet been coded and was won
Just a thought on this topic, I found Harminv quite powerful for this
sort of task. I wonder whether it could be wrapped into a R package
(it's GPL).
http://ab-initio.mit.edu/wiki/index.php/Harminv
On 20 Nov 2008, at 22:46, Prof Brian Ripley wrote:
See e.g.
http://finzi.psych.upenn.edu/R
See e.g.
http://finzi.psych.upenn.edu/R/Rhelp02a/archive/131024.html
RSiteSearch() produced this and similar relevant past postings.
On Thu, 20 Nov 2008, Ben Zuckerberg wrote:
Greetings,
I have several sets of oscillation data and would like to estimate the
parameters of a sine function to
2008/11/20 David Kaplan <[EMAIL PROTECTED]>:
> Let me try to be more specific.
>
> The x y coordinates are different because of NAs in the dataset. In this
> analysis, a set of hat values (a measure of influence in regression) is
> given for each observation. On the basis of the regression that w
Greetings,
I have several sets of oscillation data and would like to estimate the
parameters of a sine function to each set (and hopefully automate
this). A colleague provided an excel sheet that uses solver to minimize
the RSS after fitting the sine function to each data set, but this
cumbe
Is there anyway to label axes in 3D plots with mathematical expressions?
In the code below, I want to replace "delta_yrsed" with what "\Delta
\widehat{yrsed}" represents in TeX, but the [xyz]lab parameters of title3d
appear to only accept character strings.
require("rgl")
fn.delta.yrsed <- funct
Let me try to be more specific.
The x y coordinates are different because of NAs in the dataset. In
this analysis, a set of hat values (a measure of influence in
regression) is given for each observation. On the basis of the
regression that was run to get these hat values, the sample size wa
on 11/20/2008 02:28 PM Tul Gan wrote:
> Hi !
> �� I am new to R. Can somebody help me in reformatting�huge output
> files ,i.e, rearranging sets of columns in specific order.
> For example: I have data for three compunds 1, 2 and 3
> file1:
> ID CA1 CA3 CA2 MA2 MA1 MA3
> 1 14 15 13 7 12 3
On Thu, 20 Nov 2008, Sharma, Dhruv wrote:
Hi,
In a binary classification dataset is there a way to tune rpart to now
allow any false positives.
There may be no such tree, but if there is, ?rpart.control gives you the
means to achieve it.
I.e. is there some setting that will enforce no fa
Hi,
In a binary classification dataset is there a way to tune rpart to now
allow any false positives.
I.e. is there some setting that will enforce no false positives (prior
etc)? (and yet Without making the tree pick all 0 or 1 outcome)
Dhruv
[[alternative HTML version deleted]]
Actually "drm" as posted before fits a sigmoid curve (a generalized
logistic function with 4 parameters, see ?LL.4), so I didn't get the
point of your new question.
Dani Valverde schrieb:
Thank you all for your answers. If you look at the plot resulting from
my data, it seems that it is some
> txt <- "ID CA1 CA3 CA2 MA2 MA1 MA3
+ 1 14 15 13 7 12 3
+ 2 19 7 12 10 14 5
+ 3 21 12 19 6 8 9 "
> dta <- read.table(textConnection(txt), header=TRUE)
> dta
ID CA1 CA3 CA2 MA2 MA1 MA3
1 1 14 15 13 7 12 3
2 2 19 7 12 10 14 5
3 3 21 12 19 6 8 9
>
# use [ , ] with a
On 21/11/2008, at 9:19 AM, David Kaplan wrote:
Hi all,
In using the identify command, I get the following message
plot(hatvalues(scireg3))
abline(h=.0154,lty=2) # plots a reference line at (k + 1)/n
identify(1:1165, hatvalues(scireg3),row.names(sciach))
Error in xy.coords(x, y) : 'x' and '
2008/11/20 David Kaplan <[EMAIL PROTECTED]>:
> Hi all,
>
> In using the identify command, I get the following message
>
>> plot(hatvalues(scireg3))
>> abline(h=.0154,lty=2) # plots a reference line at (k + 1)/n
>> identify(1:1165, hatvalues(scireg3),row.names(sciach))
>
> Error in xy.coords(x, y) :
Hi All,
I am making the switch to R and uncertain which of the several packages for
mixed models is appropriate for my analysis. I am waiting for Pinheiro and
Bates' book to arrive via inter-library loan, but it will be a week or more
before it arrives.
I am trying to fit a generalized line
Hi,
If you wish to connect each point to the next with a different
linetype, I think your best bet is to use segments()
x <- stats::runif(12); y <- stats::rnorm(12)
i <- order(x,y); x <- x[i]; y <- y[i]
plot(x, y)
s <- seq(length(x)-1)
segments(x[s], y[s], x[s+1], y[s+1], lty=1:10)
If, how
nothing at all! i have incidentally posted to the wrong list.
apologies.
vQ
David Winsemius wrote:
>
>
> What does this have to do with R?
>
> -- David Winsemius
> Heritage Labs
> On Nov 20, 2008, at 2:45 PM, [EMAIL PROTECTED] wrote:
(...)
__
R-help@
FAQ 7.10
or use colClasses to define your input
On Thu, Nov 20, 2008 at 3:12 PM, thoeb <[EMAIL PROTECTED]> wrote:
>
> Hello, I have a problem with reading a csv-file. One colum of the inputfile
> consists of characters and numbers. After reading the csv-file I create a
> new dataframe by dividing
Hi !
I am new to R. Can somebody help me in reformatting huge output
files ,i.e, rearranging sets of columns in specific order.
For example: I have data for three compunds 1, 2 and 3
file1:
ID CA1 CA3 CA2 MA2 MA1 MA3
1 14 15 13 7 12 3
2 19 7 12 10 14 5
3 21 12 19 6 8 9
to
File 2:
I tried vif()...but I can only get it to run with lm and glm. Is there a way
to run vif() with lmer() or glmmPQL()? Or is there another way to check for
collinearity between parameters while taking a random effect into account?
> Date: Wed, 19 Nov 2008 23:26:53 -0300> From: [EMAIL PROTECTED]>
Hi all,
In using the identify command, I get the following message
> plot(hatvalues(scireg3))
> abline(h=.0154,lty=2) # plots a reference line at (k + 1)/n
> identify(1:1165, hatvalues(scireg3),row.names(sciach))
Error in xy.coords(x, y) : 'x' and 'y' lengths differ
which doesn't allow me to
I am confused by the behavior of the lines function when the lty argument is a
vector. ?lines indicates that lty is a valid parameter, but says nothing else
about it. ?plot.xy (which I think is what gets called) refers back to ?lines.
?plot.default says to see ?par. In ?par, about lty it say
Dear list,
I'm trying to get two lattice plots aligned on a page. They should
share a common x axis, hence the need for perfect alignment, but the
data is taken from unrelated, separate sources (it is therefore
inappropriate to combine them and use facetting to get an automatic
layout: t
Hello, I have a problem with reading a csv-file. One colum of the inputfile
consists of characters and numbers. After reading the csv-file I create a
new dataframe by dividing the values of that colum into more colums (then a
colum contains just characters or numbers) but the numbers are converted
On Thu, 20 Nov 2008, Wacek Kusnierczyk wrote:
[EMAIL PROTECTED] wrote:
So what exactly is your understanding of a canonical entity and
perhaps that'd clarify your point?
when i studied medicine, i used to think about 'canonical' anatomy
(well, we'd just speak of human anatomy, with no can
?colorRamp
On Thu, Nov 20, 2008 at 2:51 PM, Liu, Hao [CNTUS] <[EMAIL PROTECTED]> wrote:
> Dear All:
>
> I have a question on how to get more choices of color to show gene expression
> up or down regulation, like choosing the neutral color as yellow, any way to
> work with the color transition?
Dear All:
I have a question on how to get more choices of color to show gene expression
up or down regulation, like choosing the neutral color as yellow, any way to
work with the color transition?
I could not find any reference on those.
Thanks for your help
Hao
[[alternative HTML
Achim Zeileis wrote:
On Thu, 20 Nov 2008, David Kaplan wrote:
Hi all,
I'm looking for a program that will take the predicted probabilities
from a logistic regression using glm{stats}, dichotomize them
according to a threshold that I can control, and then use them to form
sensitivity, specif
On Thu, Nov 20, 2008 at 12:28 PM, Gavin Simpson <[EMAIL PROTECTED]> wrote:
> But Prof. Ripley has pointed out (off list) that
>
> ffg[rowSums(ffg) > 0, ]
I suggested much the same solution off-list (using apply rather than
rowSums, as I'm
apparently incapable of remembering the existence of the l
Thank you all for your answers. If you look at the plot resulting from
my data, it seems that it is some kind of sigmoid function, not only
polynomial. How could I fit it?
Best,
Dani
Daniel Valverde Saubí
Grup de Biologia Molecular de Llevats
Facultat de Veterinària de la Universitat Autònoma
On Thu, 20 Nov 2008, David Kaplan wrote:
Hi all,
I'm looking for a program that will take the predicted probabilities from a
logistic regression using glm{stats}, dichotomize them according to a
threshold that I can control, and then use them to form sensitivity,
specificity, false pos and f
Thanks.
Daniel
-
cuncta stricte discussurus
-
_
Von: Henrique Dallazuanna [mailto:[EMAIL PROTECTED]
Gesendet: Thursday, November 20, 2008 5:18 AM
An: Daniel Malter
Cc: [EMAIL PROTECTED]
Betreff: Re: [R] sub / gsub - extracting between
On Thu, 2008-11-20 at 17:08 +, Gavin Simpson wrote:
> On Thu, 2008-11-20 at 12:01 -0500, stephen sefick wrote:
> > ##I want to remove the rows where the row sums are zero and this is as
> > far as I have gotten
>
> Given your ffg,
>
> ## the which() call returns row indices for rows with rowS
I'm rather doubtful that you can improve on the uniform jittering
strategy
you originally considered. It would require intimate knowledge about
the non-uniformity of the density in the spacings between your
quantized version.
But if you really _knew_ the parent distribution
then something l
> I have some data measured with a coarsely-quantized clock. Let's say
> the real data are
>
> q<- sort(rexp(100,.5))
>
> The quantized form is floor(q), so a simple quantile plot of one
> against the other can be calculated using:
>
> plot(q,type="l"); points(floor(q),col="red")
>
Hi
[EMAIL PROTECTED] napsal dne 20.11.2008 16:15:44:
>
> I tried to do some cycles using the ?for? function, but then I had to
make
> multiple ?for?s and it didn?t work. One of the major problems is that I
> don?t know how to use the summation function on R.
I did not find function summation.
I don't know of an existing set of symbols for this, but if you can come up
with a set or descriptions of what you want, then the my.symbols function in
the TeachingDemos package can be used to place them on a plot.
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
The logspline package has tools for estimating a density function for interval
censored data (the old methods), you could use those to estimate the density of
your data, then compare that density to the theoretical density.
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain He
Hi
something like
ffg[!rowSums(ffg)==0,]
Petr Pikal
[EMAIL PROTECTED]
724008364, 581252140, 581252257
[EMAIL PROTECTED] napsal dne 20.11.2008 18:01:28:
> ##I want to remove the rows where the row sums are zero and this is as
> far as I have gotten
>
> ffg <- (structure(list(CD = c(0, 0, 0, 0
Hi all,
I'm looking for a program that will take the predicted probabilities
from a logistic regression using glm{stats}, dichotomize them according
to a threshold that I can control, and then use them to form
sensitivity, specificity, false pos and false neg rates.
Thanks in advance.
David
On Thu, 2008-11-20 at 12:01 -0500, stephen sefick wrote:
> ##I want to remove the rows where the row sums are zero and this is as
> far as I have gotten
Given your ffg,
## the which() call returns row indices for rows with rowSum > 0
ffg[which(rowSums(ffg) > 0, ]
does the trick
HTH
G
>
> ff
To be clear, the problem is not the return statement in your function, but
the extra argument, i, in your lapply statement:
lapply(1:4,fn)
works just fine with your original function. You need to read ?lapply more
carefully: fn receives the values of the first argument (1:4) in turn
automatical
lapply already passes the first arg to fn and by specifying the
i (which is undefined -- its only defined within fn) it would be
trying to to pass a second arg to fn yet fn takes only takes
one arg. Try these:
lapply(1:4, fn)
lapply(1:4, "^", 2)
On Thu, Nov 20, 2008 at 11:31 AM, megh <[EMAIL P
##I want to remove the rows where the row sums are zero and this is as
far as I have gotten
ffg <- (structure(list(CD = c(0, 0, 0, 0, 3.125, 0, 0, 0, 0, 1.6, 3.125,
0, 0, 6.25, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 3.125, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 1.6, 0, 0, 0, 0, 0,
0, 0, 0, 0
[EMAIL PROTECTED] wrote:
>
>
>>> A canonical human body will have canonical parts and those canonical
>>> parts will
>>> have canonical subparts and so on.
>>
>> ... and?
>>
>>> Can't think of anyone who would fit that
>>> description.
>>
>> is this considered an argument for that there cannot poss
on 11/20/2008 10:31 AM megh wrote:
> I have written following codes, with intention to get a list with values
> 1,2,9,16 :
>
> fn <- function(i) return(i^2)
> lapply(1:4, fn, i)
>
> However I got following error :
> Error in FUN(1:4[[1L]], ...) : unused argument(s) (1)
>
> Can anyone please tell
Hi,
you are feeding lapply "i" as an optional argument, which is passed to
fn() and causes an error. Just use lapply(1:4, fn), or better yet,
sapply,
> fn <- function(i) return(i^2)
> sapply(1:4, fn)
[1] 1 4 9 16
Hope this helps,
baptiste
On 20 Nov 2008, at 16:31, megh wrote:
I ha
rm(list=ls(all=TRUE))
??remove will get you there (I don't think 'clear' is what many people
would call this).
RON70 wrote:
Is there any function in R to clear values from all variables stored in
current and previous sessions? R site search do not give the desired result
at all, all results
I believe that the epicalc package has a "zap" function that might
accomplish that.
--
David Winsemius
On Nov 20, 2008, at 11:14 AM, RON70 wrote:
Is there any function in R to clear values from all variables stored
in
current and previous sessions? R site search do not give the desired
I have written following codes, with intention to get a list with values
1,2,9,16 :
fn <- function(i) return(i^2)
lapply(1:4, fn, i)
However I got following error :
Error in FUN(1:4[[1L]], ...) : unused argument(s) (1)
Can anyone please tell me what will be the correct code here?
Regards,
--
Another approach:
? jitter
plot(jitter(q, factor=1),type="l")
factor = 1 by default but can get increased so the spaces get filled
in to your satisfaction:
plot(q,type="l"); points( jitter(floor(q), factor=2) ,col="red")
plot(q,type="l"); points( jitter(floor(q), factor=3), col="red")
I su
On Thu, Nov 20, 2008 at 10:04 AM, Dieter Menne
<[EMAIL PROTECTED]> wrote:
> hadley wickham gmail.com> writes:
>
>> > library(plyr)
>> > dat = data.frame(SUBJECT_ID=sample(letters[1:5],100,TRUE),HR=rnorm(100))
>> > daply(dat,.(SUBJECT_ID),sd)
>> > ddply(dat,.(SUBJECT_ID),sd)
>>
>> Well that calcula
Is there any function in R to clear values from all variables stored in
current and previous sessions? R site search do not give the desired result
at all, all results are about to clear screen only not variables.
Thanks,
--
View this message in context:
http://www.nabble.com/How-to-Clear-all-v
Dear Brian, Mose, Peter and Stefan,
Thanks a lot for your replies - the issues are now clearer to me. (and
I apologize for not using the appropriate list).
Best wishes,
Emmanuel
2008/11/19 Peter Dalgaard <[EMAIL PROTECTED]>:
> Stefan Evert wrote:
>>
>> On 19 Nov 2008, at 07:56, Prof Brian Ri
Nicklas Pettersson stat.su.se> writes:
> I wonder if anyone knows how to generate a list of objects, e.g. ten
> vectors with names: vect1, vect2, ... , vect10.
> My own idea was to use something like:
>
> for (i in 1:10)
> print(paste("vect", i,"<-NULL",sep=""))
>
for (i in 1:10)
ev
hadley wickham gmail.com> writes:
> > library(plyr)
> > dat = data.frame(SUBJECT_ID=sample(letters[1:5],100,TRUE),HR=rnorm(100))
> > daply(dat,.(SUBJECT_ID),sd)
> > ddply(dat,.(SUBJECT_ID),sd)
>
> Well that calculates sd on the whole data frame. (Like sd(dat)).
Not really, it looks like the b
Thanks for your comment. I would typically follow this approach too,
but I'm wondering whether one could find a more sophisticated
solution. Ideally, I'd like to be able to select the text that is
annotating the figure. There are very few cases where I can see a real
need for raster text, t
On Thu, 20 Nov 2008, Erika Melissari wrote:
Hello all,
I have read the suggested manuals, but I do not manage to install the SPIA
package yet.
I have installed Rtools28.exe following the indications, I have put the
package tar.gz in c:\ directory and I have extracted it by using command
wind
I have some data measured with a coarsely-quantized clock. Let's say
the real data are
q<- sort(rexp(100,.5))
The quantized form is floor(q), so a simple quantile plot of one
against the other can be calculated using:
plot(q,type="l"); points(floor(q),col="red")
which of course sho
Doug Bates and I have exchanged ideas on the issue of singularities in nonlinear models a number of times over the years. Both perspectives are "right", and though I will characterize them as adversarial, they are really complementary. These views can be overly simplified as
- if there's a singula
On Thu, 2008-11-20 at 09:45 -0500, stephen sefick wrote:
> #for instance this
> ordisurf(bug.4, env.savannah[,"TSS"]+env.savannah[,"TIN.TP"])
>
> This is mod1?
Hi Stephen,
Yes, but it is mod1 expressed as
## the sum of TSS and TIN.TP
y <- env.savannah[,"TSS"]+env.savannah[,"TIN.TP"]
## sites s
I tried to do some cycles using the “for” function, but then I had to make
multiple “for”s and it didn’t work. One of the major problems is that I
don’t know how to use the summation function on R.
Best regards,
P.Branco
stephen sefick wrote:
>
> What have you tried so far?
>
> On Thu, Nov 20
Dear Gerit,
Here is a start using a data set which first column is numeric and the rest
are factors 'f1', 'f2',,'f1381' (I'm using only 3):
# Data set
x <- c(1,4,2,6,8,3,4,2,4,5,1,3)
y <- as.factor(c(2,2,1,1,1,2,2,1,1,2,1,2))
z <- as.factor(c(1,2,2,1,1,2,2,3,3,3,3,3))
mydata=data.frame(x,y,z)
On Thu, 2008-11-20 at 09:45 -0500, stephen sefick wrote:
> #for instance this
> ordisurf(bug.4, env.savannah[,"TSS"]+env.savannah[,"TIN.TP"])
>
Stephen,
According to ordisurf documentation, this is correct if 'bug.4' is an
ordination result, and the sum of those two env.savannah columns is the
si
Can someone recommend a package in R that will perform a two-sample
Kolmogorov–Smirnov test on left censored data? The package "surv2sample"
appears to offer such a test for right censored data and I guess that I can
use this package if I flip my data, but I figured I would first ask if there
was
Gerard M. Keogh justice.ie> writes:
>
>
> Hi everyone,
>
> newbee query!
>
> I've installed R 2.8.0 and tried to run this simple glm -
> x is no of cars in a given year, y is the number voted in an election
> that year while n is the population 18+:
I strongly suspect that you're co
Jim,
Thank you so much. There is a lot for me here to dig into, learn and
understand. But you have made my task so much easier by giving me sufficient
material to get started. Once again, thanks a lot.
/Ravi
- Original Message
From: jim holtman <[EMAIL PROTECTED]>
To: ravi <[EMAIL PROT
Wijffels, Jan thomascook.be> writes:
>
> Hi
>
> I'm quite new to optimization algorithms and I could use some advice or
> pointers. I'm using ?optim (method L-BFGS-B) to optimize a function over
> a 60-dimensional parameter space. The function itself takes about 1 to 6
> minutes to compute. It
What have you tried so far?
On Thu, Nov 20, 2008 at 7:40 AM, P.Branco <[EMAIL PROTECTED]> wrote:
>
> Dear all,
>
> I was wondering if you could help me with the following.
>
> I want to do calculate this equation:
>
> Ps(t)= (∑_(r=1)^N Xtr* 1/√drs)/(∑_(r=1)^N 1/√drs)
>
> Ps(t) – Probability that
Hello all,
I have read the suggested manuals, but I do not manage to install the SPIA
package yet.
I have installed Rtools28.exe following the indications, I have put the
package tar.gz in c:\ directory and I have extracted it by using command
window on c:\ folder and the following instruction
#for instance this
ordisurf(bug.4, env.savannah[,"TSS"]+env.savannah[,"TIN.TP"])
This is mod1? I am new to gam models, and will buy Simon's book when
I have the funds.
thanks for being patient
Stephen Sefick
(I will send you data off list if you wish with reproduvible code, but
I believe that it
Hi everyone,
newbee query!
I've installed R 2.8.0 and tried to run this simple glm -
x is no of cars in a given year, y is the number voted in an election
that year while n is the population 18+:
votes <- data.frame(x = c(0.62,0.77,0.71,0.74,0.77,0.86,1.13,1.44),
+
Hello,
Take a look at this course:
http://www.r4all.group.shef.ac.uk/index.html
I don't think they teach tools for working with the genome, but it
might be helpful anyway.
On Thu, Nov 20, 2008 at 11:16 AM, Peter Saffrey <[EMAIL PROTECTED]> wrote:
> (apologies if this is the wrong list)
>
> I'm
Look at summaryBy in the doBy package.
On Thu, Nov 20, 2008 at 9:16 AM, Gerit Offermann <[EMAIL PROTECTED]> wrote:
> Dear list,
>
> I reduced my data to the following:
>
> x <- c(1,4,2,6,8,3,4,2,4,5,1,3)
> y <- as.factor(c(2,2,1,1,1,2,2,1,1,2,1,2))
> z <- as.factor(c(1,2,2,1,1,2,2,3,3,3,3,3))
>
>
Dear list,
I reduced my data to the following:
x <- c(1,4,2,6,8,3,4,2,4,5,1,3)
y <- as.factor(c(2,2,1,1,1,2,2,1,1,2,1,2))
z <- as.factor(c(1,2,2,1,1,2,2,3,3,3,3,3))
I can produce the statistical summary just fine.
s1 <- tapply(x, y, summary)
d1 <- tapply(x, y, sd)
s2 <- tapply(x, z, summary)
d2
You can start by taking a look here:
https://help.ubuntu.com/8.04/serverguide/C/index.html
Daniel
Daniel Høyer Iversen
Brøsetveien 155, 14
7050 Trondheim
Mob.: 48 22 90 21
Privat: [EMAIL PROTECTED]
Skole: [EMAIL PROTECTED]
1 - 100 of 137 matches
Mail list logo
