That's not too hard:
x11(width=7,height=5)
par(cex.axis=.8)
fin_month<-month
fin_month[fin_month>6]<-fin_month[fin_month>6]-12
fin_month<-fin_month+6
plot(fin_month~year,axes=FALSE,type="l", ylab="Month", xlab="Year",
main="Month of occurrence in year")
axis(1,at=seq(1981,2014,3))
Dear Jim,
Thank you very much for nice suggestion and figure there. But what I need
is the y axis start from let say 7 (July) and end at 6 (June).
In those sense, I can said clearly that the occurrence, minima in this
case, is fall during winter season, because Jan and Dec close to each other.
Hi Ani,
There are a number of ways to modify this sort of plot. Here is one:
x11(width=7,height=5)
par(cex.axis=.8)
plot(month~year,xaxt="n", type="l", ylab="Month", xlab="Year",
main="Month of occurrence in year")
axis(1,at=seq(1981,2014,3))
library(plotrix)
Dear R-Help,
I have 35 data that is month when the annual minima happened. So I want to
plot those data but the order of y axis is not from 1 to 12, but let say
start from 9,10,11,12,1,..8. The reason to do this is when 12 (Dec) meet 1
(Jan) in the following year the graph is not quite good (for
If a crosstabs on the two factors has any zero counts it might explain.
—
David
Sent from my iPhone
> On Nov 6, 2019, at 5:42 PM, Jim Lemon wrote:
>
> Hi Greg,
> I tried this:
>
> cp<-read.table(text="Birey Grup Time y
> 11 Cp1 0.7916386
> 11 Cp3 1.7463777
> 11 Cp7 1.2008390
Hi Greg,
I tried this:
cp<-read.table(text="Birey Grup Time y
11 Cp1 0.7916386
11 Cp3 1.7463777
11 Cp7 1.2008390
11 Cp14 0.6311380
11 Cp21 2.1563557
11 Cp28 1.2008390",
header=TRUE)
library(prettyR)
xtab(Grup~Time,cp)
Crosstabulation of Grup by Time
Time
Hello,
I did Mendelian randomization using this software:
https://cran.r-project.org/web/packages/MendelianRandomization/vignettes/Vignette_MR.pdf
library(MendelianRandomization)
f=read.table("246LDout272Biobank_Retina.txt", header=T)
> head(f)
rs exposure.beta exposure.se outcome.beta
This might be impossible to answer without all the data. You might wish to
contact the package maintainer for a question like this.
Cheers,
Bert
On Wed, Nov 6, 2019 at 1:58 PM greg holly wrote:
> I got the following error message after running t2way(y ~ Grup*Time, data
> = cp)
> Error in
I got the following error message after running t2way(y ~ Grup*Time, data
= cp)
Error in x[[grp[i]]] :
attempt to select less than one element in get1index
a part of the data is given below. Your help is highly appreciated.
Greg
> head(cp)
Birey Grup Time y
1 11 Cp1
9 matches
Mail list logo