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
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
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 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
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 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
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
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
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
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,
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
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
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 &
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
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
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
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 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)
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 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 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
22 matches
Mail list logo
