Dear R users,
I have two time series
Calculate and plot cross correlation between two time series over nested time
periods. Each point in either time series is for a week (not exactly a calendar
week, but the first week in a calendar year always starts from Jan 1, and the
other weeks in the s
Hi Bob,
Given the other answers I may be off target, but perhaps this will help:
# create two "nights" worth of data
Times<-strptime(
paste(c("2015-05-13","2015-05-14"),paste(rep(c(18:23,0:6),2),":30:00",sep="")),
"%Y-%m-%d %H:%M:%S")
# telescope the two nights into repeated hours
Hours<-strptim
Dear List
I am trying to do an association network using some expression data I have, the
data is really huge: 300 samples and ~30,000 genes. I would like to apply a
gaussian graphical model to my data using the huge R package.
Here is the code I am using
> dim(data)
#[1] 317 32200
> huge.out
Hi Richard,
Thank you very much for your cooperation. I installed the packages "HH" and
run the following model:
Bud.aov<-aov(terms(Budburst~CO2*SoilTemp*Photoperiod+Greenhouse/(SoilTemp*Photoperiod),keep.order=TRUE),data=data)
> summary(Bud.aov)
Hi Richard,
Thank you very much for your cooperation. I installed the packages "HH" and
run the following model:
Bud.aov<-aov(terms(Budburst~CO2*SoilTemp*Photoperiod+Greenhouse/(SoilTemp*Photoperiod),keep.order=TRUE),data=data)
> summary(Bud.aov)
On May 13, 2015, at 2:01 PM, Vin Cheng wrote:
> Hi David,
>
> I tried both solutions you provided(lapply(dlist, function(x) rbind(x,x) ) &
> Map( rbind, smaller, smaller)) and they seem to transpose and reshape the
> data into a different structure. Whenever I added the str - I only get a
Hi David,
I tried both solutions you provided(lapply(dlist, function(x) rbind(x,x) ) &
Map( rbind, smaller, smaller)) and they seem to transpose and reshape the data
into a different structure. Whenever I added the str - I only get a NULL.
> You basically want a list of the same general str
Hmmm. Follow the instructions at
http://www.r-project.org/
After choosing a mirror, you will reach a page that says, in part:
partial quote
Source Code for all Platforms
Windows and Mac users most likely want to download the precompiled
binaries listed in the upper box, not the sou
This example might help:
tmp <- data.frame(x=1:10, y=rnorm(10))
foo <- lm(y~x, data=tmp)
plot(tmp$x, residuals(foo))
It appears that eval$bty_avg is not what you think it is.
-Don
--
Don MacQueen
Lawrence Livermore National Laboratory
7000 East Ave., L-627
Livermore, CA 94550
925-423-1062
Yes
> bar <- function(y = "12345"){
+ message(y)
+ }
> foo <- function(x = "a", y = "b"){
+ bar(y = y)
+ }
> bar()
12345
> bar(y = "abc")
abc
> foo()
b
> foo(y = "xyz")
xyz
ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek / Research Institute for Nature and
Forest
team Biometrie &
Hi,
On Wed, May 13, 2015 at 5:58 AM, venkadesan venky wrote:
> Hello Team,
>
> I have data like this
>
> have Small doubts on the following calculation
>
>
> For example,
>
> Employee sizeCamp 1camp 2 Camp 3
>
> 11 0 0
>
> 2
Hello,
I googled and looked on r-bloggers website for how to plot several REC
curves (rminer package) on the same plotting area but couldn't find
anything useful. I want to compare several regression (not logic regression
and not classification) predictive models analyzing the REC curves and the
e
Hello,
I am using shiny package to develop web application but got stuck while
creating multi-page web application. I googled my problem but did not get
proper help. I just found R is not supporting multi-page web application
but it was to old question.
*Is R supports multi-page web application d
Hello Team,
I have data like this
have Small doubts on the following calculation
For example,
Employee sizeCamp 1camp 2 Camp 3
11 0 0
20 0 1
3
I have a dataframe with the following structure:
'data.frame': 13095 obs. of 1433 variables:
$ my : Factor w/ 624 levels "19631","19632",..: 1 1 1 1 1 1 1 1 1 1 ...
$ s1 : num NA NA NA NA NA NA NA NA NA NA ...
Where my is a factor with the number of the month, s1,..,,s1426 vectors t
Your problem is that PatientID, FatherID, MotherID are factors. The authors of kinship2
(myself and Jason) simply never thought of someone doing this. Yes, that is an oversight.
We will correct it by adding some more checks and balances. For now, turn your id
variables into character or nume
Hi nusrat,
The ifelse function returns the number of values that result from the
logical expression in the first argument. If you use &&, you get one
logical value. If you use & you get logical values for the number of
conditionals that you specify. For example:
1:10 > 0 && 1:10 < 11:20
[1] TRUE
Hi Frederic,
It looks like you have a number of daily rainfall values and you want
to get the mean rainfall on days where there was some rain. What you
probably want is:
daily_rainfall<-rpois(50,2)
mean(daily_rainfall[daily_rainfall>0])
You may have to add na.rm=TRUE if there are NAs in your data
Hi Shivi82,
The error message suggests that "eval$bty_avg" is a function. What does:
str(eval)
say about its components?
Jim
On Wed, May 13, 2015 at 7:33 PM, Shivi82 wrote:
> HI All,
> I Am creating a residual plot for my linear model.
> the code I created is : plot(eval$bty_avg,residuals,yla
This problem is discussed in Winston Chang's R Graphics Handbook,
section 3.9. Adapting his code to this example:
test <- data.frame(variables = c("PE_35", "PE_49"),
value1=c(13,3),
value2=c(75,31),
value3=c(7,17),
v
HI All,
I Am creating a residual plot for my linear model.
the code I created is : plot(eval$bty_avg,residuals,ylab="residuals",
xlab="Score", main = "Residual Analysis")Here data set is eval. eval$bty_avg
is my response variable and residual is the var I have created using resid
function to stor
Dear Glenn,
Please keep the mailing list in cc.
In this case I would drop the "signature" like values and use correct
defaults.
TermStructure <- function(rates.data, method = "ns")
you can drop if(missing(method)) method = "ns" in that case.
Best regards,
ir. Thierry Onkelinx
Instituut voor n
22 matches
Mail list logo
