Re: [R] HISTOGRAM
First, a histogram would not be appropriate (your data appear to be categorical - a histogram is for continuous numeric vales) - you would need a bar plot. You should make two vectors (one for the category names and the other for the frequencies) and use the barplot function. On Fri, Nov 9, 2018 at 1:46 PM Medic wrote: > What would be the correct code (simplest version) (without gplot()) > for histogram (with 7 bars), which would include 7 names of bars under > the X-axis. The data are: > > name number > ds6277 > lk 24375 > ax46049 > dd70656 > az216544 > df 220620 > gh641827 > > (I'm attaching mydata.r, making with dput.) > > My attempt is: > > options(scipen=999) > with (mydata, hist(number)) > > P.S. I can't understand how the column "name" to include in a code > __ > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see > > https://na01.safelinks.protection.outlook.com/?url=https%3A%2F%2Fstat.ethz.ch%2Fmailman%2Flistinfo%2Fr-help&data=02%7C01%7Crab45%40pitt.edu%7C55f2571738b246f9dc1e08d646739b27%7C9ef9f489e0a04eeb87cc3a526112fd0d%7C1%7C0%7C636773859721875419&sdata=haHuxll2%2FO0Dci0fSpd0evjfi0MTmLi0JoghvHxlz3o%3D&reserved=0 > PLEASE do read the posting guide > https://na01.safelinks.protection.outlook.com/?url=http%3A%2F%2Fwww.R-project.org%2Fposting-guide.html&data=02%7C01%7Crab45%40pitt.edu%7C55f2571738b246f9dc1e08d646739b27%7C9ef9f489e0a04eeb87cc3a526112fd0d%7C1%7C0%7C636773859721875419&sdata=aL6arSyKx4m9LBdOMhdhSpsdJmGCHubfdI%2Fns4Rgytw%3D&reserved=0 > and provide commented, minimal, self-contained, reproducible code. > [[alternative HTML version deleted]] __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Re: [R] Read
Hello, I've just tested Jeff's solution, it works but the second code line should be dsh <- sh[ , -length( sh ) ] (dsh doesn't exist yet.) Hope this helps, Rui Barradas Às 02:46 de 10/11/2018, Jeff Newmiller escreveu: Your file has 5 commas in the first data row, but only 4 in the header. R interprets this to mean your first column is intended to be row names (has no corresponding column label) rather than data. (Row names are "outside" the data frame... use str(dsh) to get a better picture.) Basically, your file does not conform to consistent practices for csv files of having the same number of commas in every row. If at all possible I would eliminate the extra comma. If you have many of these broken files, you might need to read the data in pieces... e.g. dsh <- read.csv( "dat.csv", header=FALSE, skip=1 ) dsh <- dsh[ , -length( dsh ) ] dshh <- read.csv( "dat.csv", header=TRUE, nrow=1) names( dsh ) <- names( dshh ) On Fri, 9 Nov 2018, Val wrote: HI all, I am trying to read a csv file, but have a problem in the row names. After reading, the name of the first column is now "row.names" and all other column names are shifted to the right. The value of the last column become all NAs( as an extra column). My sample data looks like as follow, filename = dat.csv The first row has a missing value at column 3 and 5. The last row has a missing value at column 1 and 5 x1,x2,x3,x4,x5 12,13,,14,, 22,23,24,25,26 ,33,34,34, To read the file I used this dsh<-read.csv(file="dat.csv",sep=",",row.names=NULL,fill=TRUE,header=TRUE,comment.char = "", quote = "", stringsAsFactors = FALSE) The output from the above is dsh row.names x1 x2 x3 x4 x5 1 12 13 NA 14 NA NA 2 22 23 24 25 26 NA 3 33 34 34 NA NA The name of teh frist column is row,banes and all values of last columns is NAs However, the desired output should be x1 x2 x3 x4 x5 12 13 NA 14 NA 22 23 24 25 26 NA 33 34 34 NA How can I fix this? Thank you in advance __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. --- Jeff Newmiller The . . Go Live... DCN: Basics: ##.#. ##.#. Live Go... Live: OO#.. Dead: OO#.. Playing Research Engineer (Solar/Batteries O.O#. #.O#. with /Software/Embedded Controllers) .OO#. .OO#. rocks...1k __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Re: [R] Read
Your file has 5 commas in the first data row, but only 4 in the header. R interprets this to mean your first column is intended to be row names (has no corresponding column label) rather than data. (Row names are "outside" the data frame... use str(dsh) to get a better picture.) Basically, your file does not conform to consistent practices for csv files of having the same number of commas in every row. If at all possible I would eliminate the extra comma. If you have many of these broken files, you might need to read the data in pieces... e.g. dsh <- read.csv( "dat.csv", header=FALSE, skip=1 ) dsh <- dsh[ , -length( dsh ) ] dshh <- read.csv( "dat.csv", header=TRUE, nrow=1) names( dsh ) <- names( dshh ) On Fri, 9 Nov 2018, Val wrote: HI all, I am trying to read a csv file, but have a problem in the row names. After reading, the name of the first column is now "row.names" and all other column names are shifted to the right. The value of the last column become all NAs( as an extra column). My sample data looks like as follow, filename = dat.csv The first row has a missing value at column 3 and 5. The last row has a missing value at column 1 and 5 x1,x2,x3,x4,x5 12,13,,14,, 22,23,24,25,26 ,33,34,34, To read the file I used this dsh<-read.csv(file="dat.csv",sep=",",row.names=NULL,fill=TRUE,header=TRUE,comment.char = "", quote = "", stringsAsFactors = FALSE) The output from the above is dsh row.names x1 x2 x3 x4 x5 112 13 NA 14 NA NA 222 23 24 25 26 NA 3 33 34 34 NA NA The name of teh frist column is row,banes and all values of last columns is NAs However, the desired output should be x1 x2 x3 x4 x5 12 13 NA 14 NA 22 23 24 25 26 NA 33 34 34 NA How can I fix this? Thank you in advance __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. --- Jeff NewmillerThe . . Go Live... DCN:Basics: ##.#. ##.#. Live Go... Live: OO#.. Dead: OO#.. Playing Research Engineer (Solar/BatteriesO.O#. #.O#. with /Software/Embedded Controllers) .OO#. .OO#. rocks...1k __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
[R] Read
HI all, I am trying to read a csv file, but have a problem in the row names. After reading, the name of the first column is now "row.names" and all other column names are shifted to the right. The value of the last column become all NAs( as an extra column). My sample data looks like as follow, filename = dat.csv The first row has a missing value at column 3 and 5. The last row has a missing value at column 1 and 5 x1,x2,x3,x4,x5 12,13,,14,, 22,23,24,25,26 ,33,34,34, To read the file I used this dsh<-read.csv(file="dat.csv",sep=",",row.names=NULL,fill=TRUE,header=TRUE,comment.char = "", quote = "", stringsAsFactors = FALSE) The output from the above is dsh row.names x1 x2 x3 x4 x5 112 13 NA 14 NA NA 222 23 24 25 26 NA 3 33 34 34 NA NA The name of teh frist column is row,banes and all values of last columns is NAs However, the desired output should be x1 x2 x3 x4 x5 12 13 NA 14 NA 22 23 24 25 26 NA 33 34 34 NA How can I fix this? Thank you in advance __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Re: [R] list contingency tables
Yes, exactly (heh, heh).
?fisher.test
is probably what is wanted.
For arbitrary rxc tables with fixed marginals, this is a difficult problem.
Mehta's efficient network algorithm to solve it can be found by a web
search on "algorithm for Fisher exact test."
-- Bert
On Fri, Nov 9, 2018 at 12:36 PM David Winsemius
wrote:
> Seems like you are trying to recreate the calculations needed to perform
> an exact test. Why not look at the code for that or even easier, just
> use the function.
>
> --
>
> David.
>
> On 11/8/18 8:05 PM, li li wrote:
> > Hi all,
> >I am trying to list all the 4 by 2 tables with some fixed margins.
> >For example, consider 4 by 2 tables with row margins 1,2,2,1 and
> > column margins 3,3. I was able to do it using the code below. However,
> > as seen below, I had to first count the total number of tables with
> > the specific row margins and column margins in order to create space
> > to store the tables.
> > Is there a way to skip the step of counting the number of tables?
> >Also, wanted to avoid for loops as much as possible since it can be
> > extremely slow and inefficient.
> > Thanks so much in advance for you insight and help.
> > Hanna
> >
> >
> >
> >> library(gtools)
> >> A <- permutations(n=4,r=2,v=0:3, repeats.allowed=TRUE)
> >> B <- apply(A, 1, sum)
> >> rmg <- c(1,2,2,1)
> >> cmg <- c(3,3)
> >> m1 <- t(A[which(B==1),])
> >> m2 <- t(A[which(B==2),])
> >> m3 <- t(A[which(B==2),])
> >>
> >> ##count number of tables with row margins 1,2,2,1 and column margins
> 3,3.
> >> num <- 0
> >> for (i in 1:ncol(m1)){
> > + for (j in 1:ncol(m2)){
> > + for (k in 1:ncol(m3)){
> > + M <- t(cbind(m1[,i], m2[,j], m3[,k]))
> > + M1 <- rbind(M, cmg-apply(M,2,sum))
> > + num <- num+(sum(M1[4,] < 0) == 0)
> > + }}}
> >>
> >> #create space to store the tables
> >> C <- array(NA, dim=c(4,2,num))
> >>
> >> # list all the tables with fixed margins
> >> num <- 0
> >> for (i in 1:ncol(m1)){
> > + for (j in 1:ncol(m2)){
> > + for (k in 1:ncol(m3)){
> > + M <- t(cbind(m1[,i], m2[,j], m3[,k]))
> > + M1 <- rbind(M,cmg-apply(M,2,sum))
> > + if (sum(M1[4,] < 0) == 0) {
> > + num <- num+1
> > +C[,,num] <- M1
> > + }
> > + }}}
> >> C
> > , , 1
> >
> > [,1] [,2]
> > [1,]01
> > [2,]02
> > [3,]20
> > [4,]10
> >
> > , , 2
> >
> > [,1] [,2]
> > [1,]01
> > [2,]11
> > [3,]11
> > [4,]10
> >
> > , , 3
> >
> > [,1] [,2]
> > [1,]01
> > [2,]11
> > [3,]20
> > [4,]01
> >
> > , , 4
> >
> > [,1] [,2]
> > [1,]01
> > [2,]20
> > [3,]02
> > [4,]10
> >
> > , , 5
> >
> > [,1] [,2]
> > [1,]01
> > [2,]20
> > [3,]11
> > [4,]01
> >
> > , , 6
> >
> > [,1] [,2]
> > [1,]10
> > [2,]02
> > [3,]11
> > [4,]10
> >
> > , , 7
> >
> > [,1] [,2]
> > [1,]10
> > [2,]02
> > [3,]20
> > [4,]01
> >
> > , , 8
> >
> > [,1] [,2]
> > [1,]10
> > [2,]11
> > [3,]02
> > [4,]10
> >
> > , , 9
> >
> > [,1] [,2]
> > [1,]10
> > [2,]11
> > [3,]11
> > [4,]01
> >
> > , , 10
> >
> > [,1] [,2]
> > [1,]10
> > [2,]20
> > [3,]02
> > [4,]01
> >
> > __
> > R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
> > https://stat.ethz.ch/mailman/listinfo/r-help
> > PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> > and provide commented, minimal, self-contained, reproducible code.
>
> __
> R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide
> http://www.R-project.org/posting-guide.html
> and provide commented, minimal, self-contained, reproducible code.
>
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
Re: [R] list contingency tables
Seems like you are trying to recreate the calculations needed to perform
an exact test. Why not look at the code for that or even easier, just
use the function.
--
David.
On 11/8/18 8:05 PM, li li wrote:
Hi all,
I am trying to list all the 4 by 2 tables with some fixed margins.
For example, consider 4 by 2 tables with row margins 1,2,2,1 and
column margins 3,3. I was able to do it using the code below. However,
as seen below, I had to first count the total number of tables with
the specific row margins and column margins in order to create space
to store the tables.
Is there a way to skip the step of counting the number of tables?
Also, wanted to avoid for loops as much as possible since it can be
extremely slow and inefficient.
Thanks so much in advance for you insight and help.
Hanna
library(gtools)
A <- permutations(n=4,r=2,v=0:3, repeats.allowed=TRUE)
B <- apply(A, 1, sum)
rmg <- c(1,2,2,1)
cmg <- c(3,3)
m1 <- t(A[which(B==1),])
m2 <- t(A[which(B==2),])
m3 <- t(A[which(B==2),])
##count number of tables with row margins 1,2,2,1 and column margins 3,3.
num <- 0
for (i in 1:ncol(m1)){
+ for (j in 1:ncol(m2)){
+ for (k in 1:ncol(m3)){
+ M <- t(cbind(m1[,i], m2[,j], m3[,k]))
+ M1 <- rbind(M, cmg-apply(M,2,sum))
+ num <- num+(sum(M1[4,] < 0) == 0)
+ }}}
#create space to store the tables
C <- array(NA, dim=c(4,2,num))
# list all the tables with fixed margins
num <- 0
for (i in 1:ncol(m1)){
+ for (j in 1:ncol(m2)){
+ for (k in 1:ncol(m3)){
+ M <- t(cbind(m1[,i], m2[,j], m3[,k]))
+ M1 <- rbind(M,cmg-apply(M,2,sum))
+ if (sum(M1[4,] < 0) == 0) {
+ num <- num+1
+C[,,num] <- M1
+ }
+ }}}
C
, , 1
[,1] [,2]
[1,]01
[2,]02
[3,]20
[4,]10
, , 2
[,1] [,2]
[1,]01
[2,]11
[3,]11
[4,]10
, , 3
[,1] [,2]
[1,]01
[2,]11
[3,]20
[4,]01
, , 4
[,1] [,2]
[1,]01
[2,]20
[3,]02
[4,]10
, , 5
[,1] [,2]
[1,]01
[2,]20
[3,]11
[4,]01
, , 6
[,1] [,2]
[1,]10
[2,]02
[3,]11
[4,]10
, , 7
[,1] [,2]
[1,]10
[2,]02
[3,]20
[4,]01
, , 8
[,1] [,2]
[1,]10
[2,]11
[3,]02
[4,]10
, , 9
[,1] [,2]
[1,]10
[2,]11
[3,]11
[4,]01
, , 10
[,1] [,2]
[1,]10
[2,]20
[3,]02
[4,]01
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
Re: [R] list contingency tables
Don't give up on for loops entirely... some of the largest time savings in
optimizing loops are achieved by managing memory effectively. [1]
[1] https://www.r-bloggers.com/r-tip-use-vectormode-list-to-pre-allocate-lists
On November 8, 2018 8:05:39 PM PST, li li wrote:
>Hi all,
> I am trying to list all the 4 by 2 tables with some fixed margins.
> For example, consider 4 by 2 tables with row margins 1,2,2,1 and
>column margins 3,3. I was able to do it using the code below. However,
>as seen below, I had to first count the total number of tables with
>the specific row margins and column margins in order to create space
>to store the tables.
>Is there a way to skip the step of counting the number of tables?
> Also, wanted to avoid for loops as much as possible since it can be
>extremely slow and inefficient.
> Thanks so much in advance for you insight and help.
> Hanna
>
>
>
>> library(gtools)
>> A <- permutations(n=4,r=2,v=0:3, repeats.allowed=TRUE)
>> B <- apply(A, 1, sum)
>> rmg <- c(1,2,2,1)
>> cmg <- c(3,3)
>> m1 <- t(A[which(B==1),])
>> m2 <- t(A[which(B==2),])
>> m3 <- t(A[which(B==2),])
>>
>> ##count number of tables with row margins 1,2,2,1 and column margins
>3,3.
>> num <- 0
>> for (i in 1:ncol(m1)){
>+ for (j in 1:ncol(m2)){
>+ for (k in 1:ncol(m3)){
>+ M <- t(cbind(m1[,i], m2[,j], m3[,k]))
>+ M1 <- rbind(M, cmg-apply(M,2,sum))
>+ num <- num+(sum(M1[4,] < 0) == 0)
>+ }}}
>>
>>
>> #create space to store the tables
>> C <- array(NA, dim=c(4,2,num))
>>
>> # list all the tables with fixed margins
>> num <- 0
>> for (i in 1:ncol(m1)){
>+ for (j in 1:ncol(m2)){
>+ for (k in 1:ncol(m3)){
>+ M <- t(cbind(m1[,i], m2[,j], m3[,k]))
>+ M1 <- rbind(M,cmg-apply(M,2,sum))
>+ if (sum(M1[4,] < 0) == 0) {
>+ num <- num+1
>+C[,,num] <- M1
>+ }
>+ }}}
>>
>> C
>, , 1
>
> [,1] [,2]
>[1,]01
>[2,]02
>[3,]20
>[4,]10
>
>, , 2
>
> [,1] [,2]
>[1,]01
>[2,]11
>[3,]11
>[4,]10
>
>, , 3
>
> [,1] [,2]
>[1,]01
>[2,]11
>[3,]20
>[4,]01
>
>, , 4
>
> [,1] [,2]
>[1,]01
>[2,]20
>[3,]02
>[4,]10
>
>, , 5
>
> [,1] [,2]
>[1,]01
>[2,]20
>[3,]11
>[4,]01
>
>, , 6
>
> [,1] [,2]
>[1,]10
>[2,]02
>[3,]11
>[4,]10
>
>, , 7
>
> [,1] [,2]
>[1,]10
>[2,]02
>[3,]20
>[4,]01
>
>, , 8
>
> [,1] [,2]
>[1,]10
>[2,]11
>[3,]02
>[4,]10
>
>, , 9
>
> [,1] [,2]
>[1,]10
>[2,]11
>[3,]11
>[4,]01
>
>, , 10
>
> [,1] [,2]
>[1,]10
>[2,]20
>[3,]02
>[4,]01
>
>__
>R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
>https://stat.ethz.ch/mailman/listinfo/r-help
>PLEASE do read the posting guide
>http://www.R-project.org/posting-guide.html
>and provide commented, minimal, self-contained, reproducible code.
--
Sent from my phone. Please excuse my brevity.
__
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
Re: [R] HISTOGRAM
Hello, You probably want a bar plot, not a histogram. old.sci <- options(scipen=999) with(mydata, barplot(number, space = 0, names.arg = name, beside = TRUE)) options(scipen = old.sci) #- mydata <- read.table(text = " name number ds6277 lk 24375 ax46049 dd70656 az216544 df 220620 gh641827 ", header = TRUE) mydata Hope this helps, Rui Barradas Às 18:45 de 09/11/2018, Medic escreveu: What would be the correct code (simplest version) (without gplot()) for histogram (with 7 bars), which would include 7 names of bars under the X-axis. The data are: name number ds6277 lk 24375 ax46049 dd70656 az216544 df 220620 gh641827 (I'm attaching mydata.r, making with dput.) My attempt is: options(scipen=999) with (mydata, hist(number)) P.S. I can't understand how the column "name" to include in a code __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
[R] HISTOGRAM
What would be the correct code (simplest version) (without gplot()) for histogram (with 7 bars), which would include 7 names of bars under the X-axis. The data are: name number ds6277 lk 24375 ax46049 dd70656 az216544 df 220620 gh641827 (I'm attaching mydata.r, making with dput.) My attempt is: options(scipen=999) with (mydata, hist(number)) P.S. I can't understand how the column "name" to include in a code __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
Re: [R] randomForrest-imputation
Hi, The missRanger package performs predictive mean matching which should generate positive values if the non-missing values are positive. Regards, Denes On 11/09/2018 01:11 PM, Rebecca Bingert wrote: Hi! How can I generate only positive data with randomForrest-imputation? I'm working with laboratory values which are always positive. Can anybody help out? Thanks! (P.S.: I needed to send this request again because there were problems by delivering it) __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code. __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
[R] randomForrest-imputation
Hi! How can I generate only positive data with randomForrest-imputation? I'm working with laboratory values which are always positive. Can anybody help out? Thanks! (P.S.: I needed to send this request again because there were problems by delivering it) __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
[R] [R-pkgs] New Package pgsc
Dear all, I am pleased to announce a new CRAN package, "pgsc" version 1.0.0 < https://cran.r-project.org/package=pgsc>. This package implements the extension of the synthetic control method to allow for continuous and time-varying treatments described in Powell (2017) . Functions for both estimation and general hypothesis testing are provided. The vignette presents an extended example where panel regression suffers from omitted variables bias -- even with time and unit fixed effects -- and shows how the generalized synthetic control method produces unbiased estimates of treatment effects. The vignette can be found at < https://cran.r-project.org/web/packages/pgsc/vignettes/pgsc_vignette.pdf>. Suggestions and comments much appreciated, Philip Barrett Email: pobarrett at gmail dot com Github: https://github.com/philipbarrett/pgsc [[alternative HTML version deleted]] ___ R-packages mailing list r-packa...@r-project.org https://stat.ethz.ch/mailman/listinfo/r-packages __ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
