I'm translating a matlab routine to R and I need some equivalent to deconv():
Description: deconv()
[q,r] = deconv(v,u) deconvolves vector u out of vector v, using long
division. The quotient is returned in vector q and the remainder in
vector r such that v = conv(u,q)+r .
If u and v are vectors
See package polynom (and 'S Programming', for which it is a running
example).
On Wed, 26 Mar 2008, Ferran Carrascosa wrote:
> I'm translating a matlab routine to R and I need some equivalent to deconv():
>
> Description: deconv()
> [q,r] = deconv(v,u) deconvolves vector u out of vector v, using
I am trying to obtain adjusted means and standard errors for a three way
ANOVA
I have three effects, two continuous; fire frequency and annual
precipitation, and one categorical; soil type in an unbalanced design.
I am testing the effect of annual precipition (AP), soil type (ST), and fire
hadley wickham wrote:
> Please read this first:
> http://www.perceptualedge.com/articles/visual_business_intelligence/dual-scaled_axes.pdf
>
> It's a reasoned discussion of why it's a bad idea and proposes some
> alternative methods.
>
> Another good article is:
> K. W. Haemer. Double scales are
Dear mailing list,
I'm still quite a newbie in the statistical analysis of
genotype/allele data, resp. more generally in the analysis of
categorical variables. Moreover, I'm currently totally confused by the
many R packages available to do such analysis.
Here is my case: I've got a list of genes,
Dear all,
I was wonderin how to copy data from an Excel sheet
to the R prompt.
For a single culumn it works
to select the data in Excel,
copy to the clipboard with ctrl-c
and to use in R:
---
x<-scan()
[paste data with ctrl-v (paste)]
[ctrl-z, to end input]
---
But with multiple columns E
Use
read.delim("clipboard")
Uwe Ligges
G. Draisma wrote:
> Dear all,
> I was wonderin how to copy data from an Excel sheet
> to the R prompt.
>
> For a single culumn it works
> to select the data in Excel,
> copy to the clipboard with ctrl-c
> and to use in R:
> ---
> x<-scan()
> [paste d
On Tue, Mar 25, 2008 at 05:29:38PM -0500, hadley wickham wrote:
> Please read this first:
> http://www.perceptualedge.com/articles/visual_business_intelligence/dual-scaled_axes.pdf
Thanks for this pointer, interesting read. As an additional alternative
to dual scales, le'ts not forget about scatte
Easiest way to do it is to try it out and time it. Here is a case
where I generated two sets of data with 120,000 characters each (just
random numbers converted to character strings) and then asked for the
intersection of them. Came up with 3 matched in about 0.2 seconds.
That would seem fastest
Apologies! I though that the Orange dataset comes with R, but it is in
fact in the package "datasets".
So here's another "Orange2" dataset for the example:
Tree_v = rep(c(1:5),each = 5)
age_v = rep(seq(1,25, by = 5),5) + 10*runif(25)
circumference_v <- rep(seq(21,45, by = 5), 5)*Tree_v + 25*runif
If you want to read in the numeric values from EXCEL and the use them,
here is a function I use to go the other way -- take a POSIXct value
and convert it to what EXCEL wants. So you can just write the inverse
of this function:
unix2EXCEL <- function (time) time/86400 + 25569
On Tue, Mar 25, 20
Is this what you want? You have to add "na.rm=TRUE" to the calculations:
> week<-c(28,28,28,28,28,28,28,29,29,29,29,29,29,29,30,30,30,30,30,30,30)
> td<-c(0.015,0.012,NA,0.015,NA,0.014,0.014,0.013,0.013,0.013,0.013,NA,0.010,0.013,0.015,0.011,0.014,0.014,0.014,0.014,0.016)
> pd<-c(0,0,NA,80,NA,45,
Tribo,
It looks like geom_line() accepts only 4 linetypes and you asked for 5.
library(ggplot2)
Tree_v <- rep(c(1:5),each = 5)
age_v <- rep(seq(1,25, by = 5),5) + 10*runif(25)
circumference_v <- rep(seq(21,45, by = 5), 5)*Tree_v + 25*runif(25)
#This will work
Orange2 <- data.frame(Tree = as.fact
Hello,
I have a gene expression matrix with columns being samples and rows
being genes. I would like to display the expression values for each
gene. I have two groups which I colour differently. The aim is to see
if there is any difference between the two groups consistently across genes.
So t
try this:
dat <- matrix(rnorm(20), 10, 2)
dimnames(dat) <- list(1:10, c("sample1", "sample2"))
matplot(dat, pch = 16, xaxt = "n", col = 1:2,
ylab = "Expression Value")
axis(1, at = 1:10, labels = paste("Gene", 1:10), cex.axis = 0.7)
I hope it helps.
Best,
Dimitris
Dimitris Rizopoulos
Check this out, this might help
http://arxiv.org/abs/0801.2934
Ben Spycher
PhD student
Division of Clinical Epidemiology & Biostatistics
Institute of Social and Preventive Medicine
University of Berne
Finkenhubelweg 11 - CH-3012 Berne, Switzerland
Tel +41 31 631 35 07 / Fax +41 31
Hi,
I have some scripts creating GNUplot graphs that I would like to
convert to R. It would be a great help if there would be a
reasonable converter that takes over plot commands and initial
settings and put these into R commands (I'm an R beginner and
thus this kind of Kickstart would be helpful
Hi,
In Mathematica, if we have a 3X3 matrix A, we can get its dimensions by:
Dimensions[A] = {3,3}
TensorRank[A] = 2
In R, we have dim() for Dimensions[]. Is there any direct function similar
to TensorRank in R?
Thank you
[[alternative HTML version deleted]]
___
Hello everone,
I am a computer science researcher in the UK - and have a quick question
regarding importing data into R and more specifically latentnet.
I am trying to import data such as a text file containing: (This is based on
the Sampson data set)
---
Network attributes:
verti
On Wednesday 26 March 2008 01:31:46 pm Andreas Tille wrote:
AT> reasonable converter that takes over plot commands and initial
AT> settings and put these into R commands
Not that I know. But I doubt that someone would make this effort since the
capabilities of gnuplot and R are quite different.
Hello.
I have some further problems with modelling an
optimization problem in R:
How can I model some optimization problem in R with a
linear objective function with subject to some
nonlinear constraints?
I would like to use "optim" or "constrOptim", maybe
with respect to methods like "Simulated
Hello, my name is Alfonso, I'm trying to apply bootstrap to a process to
calculate the confidence intervals for a final value. I've revised all whole
process and any mistake has been encuntered in the script. But when I try to
apply boostrap the begining error message appears:
Error in
length(dim())?
Amit Soni wrote:
> Hi,
>
> In Mathematica, if we have a 3X3 matrix A, we can get its dimensions by:
> Dimensions[A] = {3,3}
> TensorRank[A] = 2
>
> In R, we have dim() for Dimensions[]. Is there any direct function similar
> to TensorRank in R?
>
> Thank you
>
> [[alternative
Dear all, I have two matrices:
matrix 1: height
23121 2125 25 25 2121
0 21 216767 33 22725
0 1 92 190 7171 192 067
0 00 8671 0 067
0 00 0 0 0 0 67
0
John and Rich,
Thankyou for the suggestion of I've gone over my data and found that
one of the text values had an extra space in it. Error fixed. It's good
to know what the error means, now if I encounter it again I don't need
to ask about it here.
Thanks again!
-Max
It happens that John Fox
yet another possibility...
if you look at 'help(package=xlsReadWrite)' you'll find a function
'dateTimeToStr' by which you can convert your date-numbers from column 3
into a string. Something along (untested):
myRead <- function( fn ) {
dat <- read.xls( file = fn, dateTimeAs = "numeric" )
dat
Hi,
Can someone tell me how to control sample size (n) in bootstrap function
boot in R. Can we give some option like we give for # of repeated
samples(R=say 100).
Will appreciate any help.
thanks
__
R-help@r-project.org mailing
Dear Max,
So Rich's guess that you had a factor level with only one observation (and
consequently an empty cell in the two-way table of means) proved correct.
I'll modify the plotMeans() function so that it prints a more informative
message under this circumstance.
Regards,
John
---
R 2.6.2 gets hung up when I use the "Save". It saves the R script
document but I must use "Force Quit" (not responding) and reboot after
saving.
I tried updating R and all its packages, ran the Mac Disk Utilities
program.
I hope that the following while help determine the problem.
Any su
Dear all,
I am now doing the logistic regression using R. (glm, family=binomial).
Besides the standardize summary statistics generated from R, I am also
interested in some more informations concerning the model fitting / prediction
etc; Particularly I am interested in "pseudo R squar" and "C s
Dear all,
I am now doing the logistic regression using R. (glm, family=binomial).
Besides the standardize summary statistics generated from R, I am also
interested in some more informations concerning the model fitting /
prediction etc; Particularly I am interested in "pseudo R squar" and "C
st
Greetings:
If you use LinkedIn please join our R Group.
***http://www.linkedin.com/e/gis/77616/0AEFE3574537*
Regards,
Ajit
--
_
Ajit Gemunu de Silva
Oakland CA 94619
skype: ajit_de_silva
[[alternative HTML version deleted]]
__
Dear Burak,
Since two of the explanatory variables are quantitative, it is unusual to
call this a three-way ANOVA (as opposed to a dummy-variable regression or
analysis of covariance). I'd also think about fitting the model with lm()
rather than aov(), so that you can more easily see the regressio
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] On Behalf Of Zaihra T
> Sent: Wednesday, March 26, 2008 7:57 AM
> To: Jan T. Kim; R-help@r-project.org
> Subject: [R] sample size in bootstrap(boot)
>
>
>Hi,
>
>Can someone tell me how to control sample s
Dear all, I have a problem with a loop, if anyone has any knowledge on
these things I would appreciate some comments. The code below is
designed to allow me to extract the top record of the data frame, and
them remove rows from the data frame which have an index close to the
extracted top record.
Basically you are moving the data up and then incrementing to the next
row. Here is an example; assume that you are at the 2nd entry:
1
2 <== here
3
4
Now your loop index is 2 and you remove the current data ('2') and are
left with:
1
3 <== index of 2 points here
4
Now you increment the inde
Hi Dan,
Thanks for response yes i do know that bootstrap samples generated by
function boot are of the same size as original dataset but somewhere in the
R-help threads i saw a suggestion that one can control sample size (n) by
using the following command(plz see below) but my pro
I have obtained from transfer functions, the state space matrices for the
following state space model.
x* = Ax + Bu
y = Cx + Du
I have A, B, C, and D, now I would like to take the exogenous inputs and
simulate the data using the state space model. I know there is a simulate
function in the
It all depends. If your looks like your sample, then you will
probably want to read the data with 'readLine' and bring in the entire
line without trying to parse it. You would then go through the data
with 'grep' to find the lines of interest and then you can use
'strsplit' or regular expressions
Hello there,
I'm trying to do lots of plots in one for-loop. But somehow ggplot does
not evaluate arguments as expected. Here is an example:
library(lattice)
library(ggplot2)
pl <- list()
pl2 <- list()
cDat <- as.data.frame(cbind(x1=0:100,x2=0:10,x3=1:20))
for(obs in c("x1", "x2")) {
pl[[obs]]
Hi,
As far as I can tell, The R function unix.time calculates elapsed CPU
time. Is there a command within R to measure actual elapsed time (I think
this is sometimes referred to as wall time)?
For example, the time command from GNU time calculates the actual elapsed
time, as far as I can tell
Hi,
Can anyone point me to documentation of how R objects such as vectors,
lists, matrices and data frames are laid out in memory? For example, which
of these are laid out as a 1-d array, 2-d array and so forth?
The R Internals document seems like the natural place to look, but I don't
see it
A new package 'epiR' is available on CRAN. Package description as follows:
Package: epiR
Version: 0.9-3
Date: 2008-03-24
Title: Functions for analysing epidemiological data
Author: Mark Stevenson <[EMAIL PROTECTED]> with contributions
from Telmo Nunes, Javier Sanchez, and Ron Thornton.
Ma
On Wed, Mar 26, 2008 at 1:17 PM, Andreas Klein <[EMAIL PROTECTED]> wrote:
> I have some further problems with modelling an
> optimization problem in R:
>
> How can I model some optimization problem in R with a
> linear objective function with subject to some
> nonlinear constraints?
> I would
Hi,
I would like to use odfWeave to output some contingency tables (the
output of "table") into OOo. I know I can do this in LaTex (using
"latex" in the Hmisc package), but I was wondering if it is possible in
OdfWeave. My documentation to odfTable says inputs can only be vector,
matrix or dat
Hi to the list,
I am trying to find a way to painlessly move structured data back and
forth between R and Matlab (also Octave). For this purpose I found the
R.matlab package great help. I wish to use a Matlab -v6 MAT file as an
intermediary format, because it is well read by both Matlab and
Octave
Hi R Users:
This package provides two generic functions for performing Markov
chain sampling in a naive way for a user-defined target distribution,
which involves only continuous variables. The function "gibbs_met"
performs Gibbs sampling with each 1-dimensional distribution sampled
with Metropoli
see ?system.time
Mark W. Kimpel MD ** Neuroinformatics ** Dept. of Psychiatry
Indiana University School of Medicine
15032 Hunter Court, Westfield, IN 46074
(317) 490-5129 Work, & Mobile & VoiceMail
(317) 204-4202 Home (no voice mail please)
mwkimpelgmailcom
**
Dear helpeR's,
I am estimating a hurdle model with the hurdle function of the pscl package.
(negative binomial for the count data and logit for the hurdle).
I am interested in robust estimates. In stata I can do this with robust
cluster(group). Since R is faster by roughly factor 20 in the mode
I believe that work is in progress to make odfTable work on tables as
well, in the meantime if you use "unclass" on a table you get back a
matrix that works with odfTable.
Hope this helps,
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
[EMAIL PROTECTED]
(801) 40
Where did you get that idea from? Not from the help page:
'system.time' calls the function 'proc.time', evaluates 'expr',
and then calls 'proc.time' once more, returning the difference
between the two 'proc.time' calls.
'unix.time' is an alias of 'system.time', for compatibi
Try this. If the format varies you may need to modify it from below. Replace
textConnection(rawLines) with the name of your file. This extracts lines with
a = and rereads them using read.table with = as the separator. It also
gets lines with a number that has spaces on both sides, deletes everyt
Hi R,
I have a dataframe with
dim(test)
[1] 435150 4
class(test)
[1] "data.frame"
In the third column every time a number with "-" appears I need to
replace with previous value
I am using the following code
s=which(substr(as.character(test[,3]),1,1)=="-")
for(i in 1:length(s)) test[s[i],3] =
Hi Thierry
Thanks for replying. I tried your code, but it spit an error on me:
> Orange2 <- data.frame(Tree = as.factor(Tree_v), age = age_v,
+ circumference = circumference_v, Lines =
factor(as.numeric(Orange2$Tree) + %% 4))
Error: unexpected SPECIAL in:
"Orange2 <- data.frame(Tree = as.factor(T
On Wed, Mar 26, 2008 at 11:03 AM, Sebastian Weber
<[EMAIL PROTECTED]> wrote:
> Hello there,
>
> I'm trying to do lots of plots in one for-loop. But somehow ggplot does
> not evaluate arguments as expected. Here is an example:
>
> library(lattice)
> library(ggplot2)
> pl <- list()
> pl2 <- lis
On 3/26/2008 12:20 PM, Faheem Mitha wrote:
> Hi,
>
> Can anyone point me to documentation of how R objects such as vectors,
> lists, matrices and data frames are laid out in memory? For example, which
> of these are laid out as a 1-d array, 2-d array and so forth?
Vectors, lists and matrices ar
I may have misread your original posting, I thought you were using the boot
package but I don't see any function bootstrap. I also don't see any
documentation for the sampler parameter. What package are you using? If it is
possible, I am sure I will be corrected.
Sorry, I can't be of more h
On Wed, 26 Mar 2008, Prof Brian Ripley wrote:
> Where did you get that idea from? Not from the help page:
No. From direct observation. Ie. comparing what the function actually
returned vs actual elapsed time. In one case I observed, these two are
dramatically different. It is possible I made
I am trying to generate a set of random numbers that fulfill the following
constraints:
X1 + X2 + X3 + X4 = 1
aX1 + bX2 + cX3 + dX4 = n
where a, b, c, d, and n are known.
Any function to do this?
Thanks,
-Ala'
[[alternative HTML version deleted]]
_
I am trying to generate a set of random numbers that fulfill the
following constraints:
X1 + X2 + X3 + X4 = 1
aX1 + bX2 + cX3 + dX4 = n
where a, b, c, d, and n are known.
Any function to do this?
Thanks,
-Ala'
__
R-help@r-project.org mailing list
ht
Hi all,
I recently switched to R version 2.6.2 (from 2.5.0 I was using) on Windows XP.
I have a fair amount of previous data saved in the form of "approxfun" objects
created by version 2.5.0, with code basically like this:
> x <- rlnorm(100,10,2)
> f <- approxfun(1:100,x)
> f(5)
[1] 313811.3
> sa
I have been trying to run this analysis and finally want to ask for help.
I am trying to run a paired t-test in which I have 19 subjects who have
reaction times for 5 different conditions. The data frame has 114 rows
total (6 rows for each subject representing six average reaction times by
condit
Apologies for the previous post - sent by mistake. Here is the full post
below:
I have been trying to run this analysis and finally want to ask for help.
I am trying to run a paired t-test in which I have 19 subjects who have
reaction times for 5 different conditions. The data frame has 114 row
Hey everyone,
For my research, I needed a peak detector algorithm that can spot peaks in
an intensity graph. The excel
http://www.nabble.com/file/p16313309/sample%2Bdata.xls sample data that
I've attached contains columns of such intensity graphs.
I came found an algorithm developed by the
At 01:13 PM 3/26/2008, Ala' Jaouni wrote:
>I am trying to generate a set of random numbers that fulfill the following
>constraints:
>
>X1 + X2 + X3 + X4 = 1
>
>aX1 + bX2 + cX3 + dX4 = n
>
>where a, b, c, d, and n are known.
>
>Any function to do this?
1. Generate random variates for X1, X2, based
Hi,
I failed to mention that the X values have to be positive and between 0 and 1.
e.g.
0.1812*X1 + 0.1871*X2 + 0.1847*X3 + 0.2745*X4 + 0.1304*X5 = 0.2
so one possible combination of X values can be:
0.319, 0.201, 0.084, 0.26, 0.136
another possible combination:
0.151, 0.253, 0.197, 0.256, 0.14
You have 4 random variables that satisfy 2 linear constraints, so you
are trying to generate a point in a (4-2) = 2 dimensional linear
(affine, in fact) subspace of R^4.
If you don't have any further requirement for the distribution of the
random points you want to generate, there are infinitely
Hi Chaser --
Ask on the Bioconductor mailing list
(http://www.bioconductor.org/docs/mailList.html) to connect with other
Bioconductor users likely to have experience with the package.
Use an informative subject line, so that people able to help can see
that the help they can provide is needed.
On Wed, Mar 26, 2008 at 7:27 PM, Ala' Jaouni <[EMAIL PROTECTED]> wrote:
> I failed to mention that the X values have to be positive and between 0 and
> 1.
Use Robert's method, and to do his step 1, use runif (?runif) to get
random numbers from the uniform distribution between 0 and 1.
Paul
___
On 26-Mar-08 20:13:50, Robert A LaBudde wrote:
> At 01:13 PM 3/26/2008, Ala' Jaouni wrote:
>>I am trying to generate a set of random numbers that fulfill
>>the following constraints:
>>
>>X1 + X2 + X3 + X4 = 1
>>
>>aX1 + bX2 + cX3 + dX4 = n
>>
>>where a, b, c, d, and n are known.
>>
>>Any function
Hello,
I want to draw 3D cylinder objects in R.
Given is the length and the diameter of the cylinder.
Has anybody an example?
Thank you very much!
Best regards
--
--
Dr. Hans-Joachim Klemmt
Forstoberrat
Organisationsprogrammierer IHK
Bayerische Landesanstalt für Wald
A friend of mine can't send emails to the
R-list from his work and he had a question that
he asked me to send because I don't know the answer
to his question. I did suggest the cran packages
list and the search function but I think he wasn't
successful ?
He needs to do a quick prototype in R with
Hi R,
I have a list
> class(pp2)
[1] "list"
> length(pp2)
[1] 1244
It is in the below format
RIC Trade.Date Close.Price Currency.Code Convertion.Rate New.Price
ABCD.SZ 2008/02/29 15.30CNY 0.1408 2.154240
ABCD.SZ 2008/01/31 15.27CNY 0.1392
http://cran.r-project.org/contrib/extra/dcom/00ReadMe.html
is the easiest way.
On Wed, 26 Mar 2008, [EMAIL PROTECTED] wrote:
> A friend of mine can't send emails to the
> R-list from his work and he had a question that
> he asked me to send because I don't know the answer
> to his question. I di
Not exactly clear what the data in the list looks like. Is it a
dataframe of the format you listed? Do you want just a single range
for the dates? Do you want to create one large dataframe and then
partition it by date to find the range? You can do
"do.call(rbind,pp2)" to create the one large d
Can you give a little more detail?
Do you want the nice 3d shading effect? Or is an ellipse, 2 lines and another
ellipse good enough? Does the viewing angle matter (flatness of ellipses). Do
the cylinders need to be placed at certain locations? Or are you trying to do a
barplot with 3d cylind
To add more clarity to my question
My data pp2 is a list
(pp2[[1]])
RIC Trade.Date Close.Price Currency.Code Convertion.Rate New.Price
ABCD.SZ 2008/02/29 15.30 CNY 0.1408
2.154240
ABCD.SZ 2008/01/31 15.27 CNY 0.1392
2.040048
ABCD.SZ
I think something like this should work. I will give you the range
for each date across all the data:
x <- do.call(rbind, pp2)
tapply(x$New.Price, x$Trade.Date, range)
On 3/26/08, Ravi S. Shankar <[EMAIL PROTECTED]> wrote:
> To add more clarity to my question
>
> My data pp2 is a list
>
> (pp2
X1,X2,X3,X4 should have independent distributions. They should be
between 0 and 1 and all add up to 1. Is this still possible with
Robert's method?
Thanks
On Wed, Mar 26, 2008 at 12:52 PM, Ted Harding
<[EMAIL PROTECTED]> wrote:
> On 26-Mar-08 20:13:50, Robert A LaBudde wrote:
> > At 01:13 PM 3/2
Dear list,
Is there any way of making barplots as a Trellis graphic?
Currently I do something like:
class <- unique(mydata[,1])
par(mfrow=c(3,3))
for(i in class){
sel <- mydata[mydata[,1]==i,]
barplot(sel)
title(i)
}
where class would take values between 1 and 9
and the first column of myd
X1,X2,X3,X4 should have independent distributions. They should be
between 0 and 1 and all add up to 1. Is this still possible with
Robert's method?
Thanks
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read t
Ala' Jaouni wrote:
>
> I am trying to generate a set of random numbers that fulfill the
> following constraints:
>
> X1 + X2 + X3 + X4 = 1
>
> aX1 + bX2 + cX3 + dX4 = n
>
> where a, b, c, d, and n are known.
>
> Any function to do this?
>
You must give more information.
How are those numbers
Hi,
I need to compute an array from a matrix and an array:
A <- array(1:20,c(2,2,5))
B <- matrix(1:10,5)
And I would like the result to be an array consisting of the following:
rbind(A[1,,1]*B[1,],
A[2,,1]*B[1,])
rbind(A[1,,2]*B[2,],
A[2,,2]*B[2,])
rbind(A[1,,3]*B[2,],
A[2,,3]*B[2,])
etc.
He
On 26-Mar-08 21:26:59, Ala' Jaouni wrote:
> X1,X2,X3,X4 should have independent distributions. They should be
> between 0 and 1 and all add up to 1. Is this still possible with
> Robert's method?
>
> Thanks
I don't think so. A whileago you wrote
"The numbers should be uniformly distributed" (but
In the code of my previous message
the barplot should be:
barplot(apply(sel[,-1],2,mean))
instead of
barplot(sel)
Sorry for the confusion.
Agus
Mensaje original
Asunto: barplot as Trellis graphic
Fecha: Wed, 26 Mar 2008 22:24:04 +0100
De: Agustin Lobo <[EMAIL PROTECTED]>
Respon
Hi,
I've been searching CRAN and the web for a recursive multivariate
filter with time-varying coefficients.
What I mean is the following:
I have a series of square matrices A_t
an initial value vector y_0
and I need to compute
y_t =A_t%*%y_t-1
As these y_t may diverge quickly and/or lead t
I did the following
DF<-do.call(rbind, pp2)
DF1=na.omit(DF)
DF1[,2]=as.Date(DF1[,2])
str(DF)
'data.frame': 18660 obs. of 6 variables:
I tried the following code
temp<-lapply(split(DF1,DF1$Trade.Date), function(.df) {
+
data.frame(DATE=.df$Trade.Date,RANGE=max(.df$New.Price)-min(.df$New.Pric
e
Hi
When I copy a block of code to R console, the code would be displayed
line by line very slowly.
I am using R version 2.6.2(2008-2-8). Does anybody know what is the
reason for the slow display? It was not like this when I used earlier
versions.
Huilin
__
Dear all,
I an running R on a Windows 2000 machine (1.5Gb RAM) and am trying to load the
lme4 package, however, whenever I attempt to load the library "lme4" I get the
following error message
I have installed lme4 and Matrix locally from the zip file with no problems.
R version 2.6.2 (2008-0
If you have lots of memory there is an obvious strategy:
d12 <- prod(dim(A)[1:2])
A <- A * array(rep(B, each = d12), dim = dim(A))
I don't really see much wrong with the obvious for() loop, though:
for(b in 1:length(B)) A[,,b] <- A[,,b] * B[b]
Bill Venables
CSIRO Laboratories
PO Box 120, Cleve
Mac OSX 10.4.10 and R 2.6.2
Whenever I resize a quartz window the graph disapears. I have to size
the window and then reprint the graph. Should I report this to the
R-developer's list? Are there any work arounds?
thanks
Stephen Sefick
--
Let's not spend our time and resources thinking about
OOPS! A mistake below. I should have written:
This raises a general question: Does anyone know of
an R function to sample uniformly in the interior
of a general (k-r)-dimensional simplex embedded in
k dimensions, with (k-r+1) given vertices?
On 26-Mar-08 22:06:54, Ted Harding wrote:
> On
This statement:
temp<-lapply(split(DF1,DF1$Trade.Date), function(.df) {
+
data.frame(DATE=.df$Trade.Date,RANGE=max(.df$New.Price)-min(.df$New.Pric
e))
+ })
has many results for DATE and one for RANGE; that is the reason you
are getting multiple copies. Instead you need to write:
temp<-lapply(sp
Hi,
I am creating a model with loess but predict complains about missing
variables. This doesn't help with single variable formulas or with lm with
any formula. Looks like a bug to me
# doesn't work
> l = loess(y ~ log(x) + log(z) , data.frame(z = 1:100, x=1:100 +
rnorm(100), y=1:100))
> predict(
>From: "Wang, Xiaojing" <[EMAIL PROTECTED]>
>Date: 2008/03/26 Wed AM 10:08:01 CDT
>To: "'r-help@r-project.org'"
>Subject: [R] pseudo R square and/or C statistic in R logistic regression
The analogue to the regression R rsquared in a GLM
is the deviance. I don't know what C is but hopefully someo
Ala' Jaouni gmail.com> writes:
>
> X1,X2,X3,X4 should have independent distributions. They should be
> between 0 and 1 and all add up to 1. Is this still possible with
> Robert's method?
>
NO.
If they add to 1 they are not independent.
As Ted remarked, the constraints define two simplexes an
Check out the 'rgl' package for creating 3D objects.
On Wed, Mar 26, 2008 at 3:02 PM, Hans-Joachim Klemmt
<[EMAIL PROTECTED]> wrote:
> Hello,
>
> I want to draw 3D cylinder objects in R.
>
> Given is the length and the diameter of the cylinder.
>
> Has anybody an example?
>
> Thank you very much!
Is there an on-line source that explains index vectors? I used 'strsplit' to
generate a lengthy list of pairs of character vectors. I only want the first
of each pair. Brackets, double brackets, and '$' operator don't work
[[alternative HTML version deleted]]
_
Is this what you were asking for?
> x <- c('a b', 'c d', 'e f')
> y <- strsplit(x, ' ')
> y
[[1]]
[1] "a" "b"
[[2]]
[1] "c" "d"
[[3]]
[1] "e" "f"
> sapply(y, '[[', 1)
[1] "a" "c" "e"
>
>
On Wed, Mar 26, 2008 at 8:55 PM, <[EMAIL PROTECTED]> wrote:
> Is there an on-line source that explains in
I find R books to be helpful for beginners. There is a fairly
comprehensive list on the R web site, www.r-project.org.
Is this what you had in mind though? Without example code, it's hard
for me to tell exactly your situation.
test <- list(first = c("string1a", "string2a"), second = c("string1
1 - 100 of 110 matches
Mail list logo
