[R] Shaded area
Dear R users, I have an xlsx file (attached to this mail) that shows the values of a "der" series observed on a daily basis from January 1, 2017 to January 25, 2017. This series is strictly positive during two periods: from January 8, 2017 to January 11, 2017 and from January 16, 2017 to January 20, 2017. I would like to plot the series with two shaded areas corresponding to the positivity of the series. Specifically, I would like to draw 4 vertical lines intersecting the x-axis in the 4 dates mentioned above and shade the two areas of positivity. Thanks for your help. __ 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] Logistic regression for large data
Dear R users,
I have a database called Base.csv (attached to this email) which
contains 13 columns and 8257 rows and whose the first 8 columns are dummy
variables which take 1 or 0. The problem is when I wrote the following
instructions to do a logistic regression , R runs for hours and hours
without giving an output:
Base=read.csv("C:\\Users\\HP\\Desktop\\New\\Base.csv",header=FALSE,sep=";")
fit_1=glm(Base[,2]~Base[,1]+Base[,10]+Base[,11]+Base[,12]+Base[,13],family=binomial(link="logit"))
Apparently, there is not enough memory to have the requested output. Is
there any other function for logistic regression that handle large data and
return output in reasonable time.
Many thanks
Kind regards
George
__
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] Mapping of groups of countries
Thanks a lot Jim. I really appreciate your help. Happy new year 2022.
George
On Wed, Dec 29, 2021 at 10:10 PM Jim Lemon wrote:
> Hi George,
> This example is obviously a mess, but it shows the basic elements. The
> labels may be placed manually using something like boxed.labels
> (plotrix) with lines connecting the labels to the smaller countries.
> If you want the labels within the countries you will have to use a big
> graphics device.
>
> library(maps)
> map("world")
> countries=c("China","Bangladesh","Egypt","Indonesia","South
> Korea","Nigeria",
> "Brazil","Pakistan","Philippines","Vietnam","South Africa",
> "India","Mexico","Turkey")
> countrycol<-c(rep("red",6),rep("blue",4),"green",rep("orange",3))
> for(i in 1:length(countries)) {
> mapinf<-map(region=countries[i],col=countrycol[i],add=TRUE)
> xpos<-sum(mapinf$range[1:2])/2
> ypos<-sum(mapinf$range[3:4])/2
> text(xpos,ypos,countries[i])
> }
>
>
> Jim
>
> On Thu, Dec 30, 2021 at 7:34 AM george brida
> wrote:
> >
> > Dear R users,
> >
> > I would like to map with R the following countries (with their names on
> the
> > map) as follows: the first group in red composed by China, Bangladesh,
> > Eypt, Indonesia, South Korea, Nigeria. The second group in blue composed
> by
> > Brazil, Pakistan, Philippines, Vietnam.
> >
> > The third group in green solely composed by South Africa.
> >
> > Finally, the fourth group in orange composed by India, Mexico and Turkey.
> >
> > I'd also like a legend that assigns the different colors to the various
> > group.
> > I tried several times with library maps but it seems that I failed to
> find
> > the real setting of xlim and ylim. Can you please help me, I will be very
> > grateful.
> >
> > Best regards
> > George
> >
> > [[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.
>
[[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.
[R] Mapping of groups of countries
Dear R users, I would like to map with R the following countries (with their names on the map) as follows: the first group in red composed by China, Bangladesh, Eypt, Indonesia, South Korea, Nigeria. The second group in blue composed by Brazil, Pakistan, Philippines, Vietnam. The third group in green solely composed by South Africa. Finally, the fourth group in orange composed by India, Mexico and Turkey. I'd also like a legend that assigns the different colors to the various group. I tried several times with library maps but it seems that I failed to find the real setting of xlim and ylim. Can you please help me, I will be very grateful. Best regards George [[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] Mapping of countries
Thanks a lot Jim..
Thank you Ani and Rasmus.
Best
On Wed, Apr 1, 2020 at 2:54 AM ani jaya wrote:
> Maybe simply add:
>
>
> points(station$Lon, station$Lat, col="red", pch=16, label="Your Country")
> text(station$Lon, station$Lat,"Your Country",
> col="black", pos=3, cex=1)
>
> station$Lon and Lat in your coordinate position.
>
> Regards,
> Ani
>
>
> On Tue, Mar 31, 2020 at 10:58 PM george brida
> wrote:
> >
> > Dear Jim,
> >
> > Is it possible to add also a title to this map?
> >
> > Many thanks
> >
> >
> >
> > On Tue, Mar 31, 2020 at 2:29 PM george brida
> wrote:
> >
> > > Dear Jim,
> > >
> > > Thank you very much. I obtained now the required map. I would like to
> know
> > > how to add the names of the countries.
> > >
> > > Best
> > > George
> > >
> > > On Tue, Mar 31, 2020 at 10:10 AM Jim Lemon
> wrote:
> > >
> > >> Hi George,
> > >> Try this:
> > >>
> > >> library(maps)
> > >> map("world",xlim=c(34.353,60.369),ylim=c(16.7,32.193),
> > >> regions="Saudi Arabia",col="yellow",fill=TRUE)
> > >> map("world",regions="Bahrain",col="yellow",fill=TRUE,add=TRUE)
> > >> map("world",regions="Kuwait",col="lightblue",fill=TRUE,add=TRUE)
> > >> map("world",regions="Qatar",col="yellow",fill=TRUE,add=TRUE)
> > >> map("world",regions="United Arab
> > >> Emirates",col="lightblue",fill=TRUE,add=TRUE)
> > >> map("world",regions="Oman",col="lightgreen",fill=TRUE,add=TRUE)
> > >>
> > >> Jim
> > >>
> > >> On Tue, Mar 31, 2020 at 12:39 PM george brida >
> > >> wrote:
> > >> >
> > >> > Dear R users,
> > >> >
> > >> > i would like to plot the maps of the Gulf Cooperation Council (GCC)
> > >> > countries (KSA, Qatar, Bahrain, Kuwait, UAE and Oman) with these
> > >> > constraints: i/ KSA , Qatar and Bahrain have the same face color ,
> ii/
> > >> > Kuweit and UAE with the same face color and iii/Oman with another
> face
> > >> > color. Is there any code in R doing this task.
> > >> >
> > >> > Many thanks.
> > >> >
> > >> > [[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.
> > >>
> > >
> >
> > [[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.
>
[[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] Mapping of countries
Dear Jim,
Is it possible to add also a title to this map?
Many thanks
On Tue, Mar 31, 2020 at 2:29 PM george brida wrote:
> Dear Jim,
>
> Thank you very much. I obtained now the required map. I would like to know
> how to add the names of the countries.
>
> Best
> George
>
> On Tue, Mar 31, 2020 at 10:10 AM Jim Lemon wrote:
>
>> Hi George,
>> Try this:
>>
>> library(maps)
>> map("world",xlim=c(34.353,60.369),ylim=c(16.7,32.193),
>> regions="Saudi Arabia",col="yellow",fill=TRUE)
>> map("world",regions="Bahrain",col="yellow",fill=TRUE,add=TRUE)
>> map("world",regions="Kuwait",col="lightblue",fill=TRUE,add=TRUE)
>> map("world",regions="Qatar",col="yellow",fill=TRUE,add=TRUE)
>> map("world",regions="United Arab
>> Emirates",col="lightblue",fill=TRUE,add=TRUE)
>> map("world",regions="Oman",col="lightgreen",fill=TRUE,add=TRUE)
>>
>> Jim
>>
>> On Tue, Mar 31, 2020 at 12:39 PM george brida
>> wrote:
>> >
>> > Dear R users,
>> >
>> > i would like to plot the maps of the Gulf Cooperation Council (GCC)
>> > countries (KSA, Qatar, Bahrain, Kuwait, UAE and Oman) with these
>> > constraints: i/ KSA , Qatar and Bahrain have the same face color , ii/
>> > Kuweit and UAE with the same face color and iii/Oman with another face
>> > color. Is there any code in R doing this task.
>> >
>> > Many thanks.
>> >
>> > [[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.
>>
>
[[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] Mapping of countries
Dear Gunter, It is noted. Thanks Best On Tue, Mar 31, 2020 at 6:04 AM Bert Gunter wrote: > Probably better posted to the r-sig-geo list, where you are more > likely to find the relevant expertise. > > Perhaps see also https://cran.r-project.org/web/views/Spatial.html, > depending on what you plan to do with your maps. > > Bert Gunter > > "The trouble with having an open mind is that people keep coming along > and sticking things into it." > -- Opus (aka Berkeley Breathed in his "Bloom County" comic strip ) > > On Mon, Mar 30, 2020 at 8:03 PM Erin Hodgess > wrote: > > > > Hello George! > > > > Do you mean to have a map of the world with these countries filled in, or > > to have just those countries on the map, please? > > > > Thanks, > > Erin > > > > Erin Hodgess, PhD > > mailto: erinm.hodg...@gmail.com > > > > > > On Mon, Mar 30, 2020 at 7:38 PM george brida > wrote: > > > > > Dear R users, > > > > > > i would like to plot the maps of the Gulf Cooperation Council (GCC) > > > countries (KSA, Qatar, Bahrain, Kuwait, UAE and Oman) with these > > > constraints: i/ KSA , Qatar and Bahrain have the same face color , ii/ > > > Kuweit and UAE with the same face color and iii/Oman with another face > > > color. Is there any code in R doing this task. > > > > > > Many thanks. > > > > > > [[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. > > > > > > > [[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. > [[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] Mapping of countries
Dear Jim,
Thank you very much. I obtained now the required map. I would like to know
how to add the names of the countries.
Best
George
On Tue, Mar 31, 2020 at 10:10 AM Jim Lemon wrote:
> Hi George,
> Try this:
>
> library(maps)
> map("world",xlim=c(34.353,60.369),ylim=c(16.7,32.193),
> regions="Saudi Arabia",col="yellow",fill=TRUE)
> map("world",regions="Bahrain",col="yellow",fill=TRUE,add=TRUE)
> map("world",regions="Kuwait",col="lightblue",fill=TRUE,add=TRUE)
> map("world",regions="Qatar",col="yellow",fill=TRUE,add=TRUE)
> map("world",regions="United Arab
> Emirates",col="lightblue",fill=TRUE,add=TRUE)
> map("world",regions="Oman",col="lightgreen",fill=TRUE,add=TRUE)
>
> Jim
>
> On Tue, Mar 31, 2020 at 12:39 PM george brida
> wrote:
> >
> > Dear R users,
> >
> > i would like to plot the maps of the Gulf Cooperation Council (GCC)
> > countries (KSA, Qatar, Bahrain, Kuwait, UAE and Oman) with these
> > constraints: i/ KSA , Qatar and Bahrain have the same face color , ii/
> > Kuweit and UAE with the same face color and iii/Oman with another face
> > color. Is there any code in R doing this task.
> >
> > Many thanks.
> >
> > [[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.
>
[[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] Mapping of countries
Hi Erin, Thanks for the reply. I would like to have just those countries on the map. Best George On Tue, Mar 31, 2020 at 5:02 AM Erin Hodgess wrote: > Hello George! > > Do you mean to have a map of the world with these countries filled in, or > to have just those countries on the map, please? > > Thanks, > Erin > > Erin Hodgess, PhD > mailto: erinm.hodg...@gmail.com > > > On Mon, Mar 30, 2020 at 7:38 PM george brida > wrote: > >> Dear R users, >> >> i would like to plot the maps of the Gulf Cooperation Council (GCC) >> countries (KSA, Qatar, Bahrain, Kuwait, UAE and Oman) with these >> constraints: i/ KSA , Qatar and Bahrain have the same face color , ii/ >> Kuweit and UAE with the same face color and iii/Oman with another face >> color. Is there any code in R doing this task. >> >> Many thanks. >> >> [[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. >> > [[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.
[R] Mapping of countries
Dear R users, i would like to plot the maps of the Gulf Cooperation Council (GCC) countries (KSA, Qatar, Bahrain, Kuwait, UAE and Oman) with these constraints: i/ KSA , Qatar and Bahrain have the same face color , ii/ Kuweit and UAE with the same face color and iii/Oman with another face color. Is there any code in R doing this task. Many thanks. [[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.
[R] diff doesnot work
Dear R users,
I have a txt file entitled coc composed by one column of numeric values
without header and having 50 rows. This file is under the following path:
C:\\Users\\intel\\Documents\\TR
I have written the following lines:
xcx=read.table("C:\\Users\\intel\\Documents\\TR\\coc.txt",header=F)
xc=log(xcx)
I have obtained the series xc but when I wrote
x=diff(xc)
I obtained the following message:
data frame with 0 columns and 50 rows
I don't know where is the problem. Can you help me please
[[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] Matrix multiplications
Thank you very much Peter. On Sat, May 21, 2016 at 9:18 PM, peter dalgaard <pda...@gmail.com> wrote: > > > On 21 May 2016, at 21:00 , george brida <george.br...@gmail.com> wrote: > > > > Dear R users: > > > > I have written the following lines : > > > >> x=c(10,11,12,13,14,17,15,16,10,11,41,25,26,14,12,14,15,20,14,22) > >> x=matrix(x,ncol=2) > >> a=matrix(1,nrow(x),1) > >> X=cbind(a,x) > >> y=c(12.00, 11.00, 13.00, 12.50, 14.00, 18.50, 15.00, 12.50, 13.75, > 15.00) > > > >> b=solve(t(X)%*% X)%*% t(X)%*% y > > > > when I wrote the following line > >> (t(y-X %*% b)%*%(y-X %*% b)/(length(y)-ncol(X)))*solve(t(X)%*% X) > > I have obtained an error message, I don't know why namely (t(y-X %*% > > b)%*%(y-X %*% b)/(length(y)-ncol(X))) is a scalar: > >> (t(y-X %*% b)%*%(y-X %*% b)/(length(y)-ncol(X))) > > [,1] > > [1,] 3.620354 > > > > > > Can you please help me. > > Mistake is that (t(y-X %*% b)%*%(y-X %*% b)/(length(y)-ncol(X))) is NOT a > scalar, but a 1 x 1 matrix. > > This works: > > as.vector((t(y-X %*% b)%*%(y-X %*% b))/(length(y)-ncol(X)))*solve(t(X)%*% > X) > > as does recognizing the first term as a sum of squares: > > sum((y-X %*% b)^2)/(length(y)-ncol(X))*solve(t(X)%*% X) > > (And, for the Illuminati, a Kronecker product works too: > > (t(y-X %*% b)%*%(y-X %*% b)/(length(y)-ncol(X))) %x% solve(t(X)%*% X) > > This could be useful for multivariate y.) > > > > > Thank you > > > > George > > > > [[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. > > -- > Peter Dalgaard, Professor, > Center for Statistics, Copenhagen Business School > Solbjerg Plads 3, 2000 Frederiksberg, Denmark > Phone: (+45)38153501 > Office: A 4.23 > Email: pd@cbs.dk Priv: pda...@gmail.com > > > > > > > > > > [[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] Matrix multiplications
Roy, Yes, t(y-X %*% b) is the transpose of y-X %*% b. In principle the product of t(y-X %*% b) *(y-X %*% b) is a scalar, I don't know why I have an error message after the following line: (t(y-X %*% b)%*%(y-X %*% b)/(length(y)-ncol(X)))*solve(t(X)%*% X) Thanks On Sat, May 21, 2016 at 9:10 PM, Roy Mendelssohn - NOAA Federal < roy.mendelss...@noaa.gov> wrote: > > str(t(y-X %*% b)) > num [1, 1:10] 0.595 -1.7538 -0.0498 -1.651 -0.6328 ... > > str((y-X %*% b)) > num [1:10, 1] 0.595 -1.7538 -0.0498 -1.651 -0.6328 … > > -Roy > > > > On May 21, 2016, at 12:00 PM, george brida <george.br...@gmail.com> > wrote: > > > > Dear R users: > > > > I have written the following lines : > > > >> x=c(10,11,12,13,14,17,15,16,10,11,41,25,26,14,12,14,15,20,14,22) > >> x=matrix(x,ncol=2) > >> a=matrix(1,nrow(x),1) > >> X=cbind(a,x) > >> y=c(12.00, 11.00, 13.00, 12.50, 14.00, 18.50, 15.00, 12.50, 13.75, > 15.00) > > > >> b=solve(t(X)%*% X)%*% t(X)%*% y > > > > when I wrote the following line > >> (t(y-X %*% b)%*%(y-X %*% b)/(length(y)-ncol(X)))*solve(t(X)%*% X) > > I have obtained an error message, I don't know why namely (t(y-X %*% > > b)%*%(y-X %*% b)/(length(y)-ncol(X))) is a scalar: > >> (t(y-X %*% b)%*%(y-X %*% b)/(length(y)-ncol(X))) > > [,1] > > [1,] 3.620354 > > > > > > Can you please help me. > > > > Thank you > > > > George > > > > [[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. > > ** > "The contents of this message do not reflect any position of the U.S. > Government or NOAA." > ** > Roy Mendelssohn > Supervisory Operations Research Analyst > NOAA/NMFS > Environmental Research Division > Southwest Fisheries Science Center > ***Note new address and phone*** > 110 Shaffer Road > Santa Cruz, CA 95060 > Phone: (831)-420-3666 > Fax: (831) 420-3980 > e-mail: roy.mendelss...@noaa.gov www: http://www.pfeg.noaa.gov/ > > "Old age and treachery will overcome youth and skill." > "From those who have been given much, much will be expected" > "the arc of the moral universe is long, but it bends toward justice" -MLK > Jr. > > [[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.
[R] Matrix multiplications
Dear R users: I have written the following lines : >x=c(10,11,12,13,14,17,15,16,10,11,41,25,26,14,12,14,15,20,14,22) > x=matrix(x,ncol=2) > a=matrix(1,nrow(x),1) > X=cbind(a,x) >y=c(12.00, 11.00, 13.00, 12.50, 14.00, 18.50, 15.00, 12.50, 13.75, 15.00) >b=solve(t(X)%*% X)%*% t(X)%*% y when I wrote the following line >(t(y-X %*% b)%*%(y-X %*% b)/(length(y)-ncol(X)))*solve(t(X)%*% X) I have obtained an error message, I don't know why namely (t(y-X %*% b)%*%(y-X %*% b)/(length(y)-ncol(X))) is a scalar: > (t(y-X %*% b)%*%(y-X %*% b)/(length(y)-ncol(X))) [,1] [1,] 3.620354 Can you please help me. Thank you George [[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.
[R] Problem with plm package
Dear friends, I have a csv data entitled GFS in the following path: c:\TEG This data is attached with this mail. When I wrote the following lines: y=read.csv(c:\\TEG\\GFS.csv, header=FALSE, sep=;) purtest(y, pmax = 4, exo = intercept, test = hadri) I obtained the following message: Error in Math.factor(x) : cumsum not meaningful for factors In addition: There were 50 or more warnings (use warnings() to see the first 50) where's the problem? can you please help me Thank you in advance George __ R-help@r-project.org mailing list 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] Problem with plm package
Hi Arun you can't open the file? On Wed, Jan 29, 2014 at 9:23 PM, arun smartpink...@yahoo.com wrote: Hi, Your attachment didn't came through. A.K. On Wednesday, January 29, 2014 12:03 PM, george brida george.br...@gmail.com wrote: Dear friends, I have a csv data entitled GFS in the following path: c:\TEG This data is attached with this mail. When I wrote the following lines: y=read.csv(c:\\TEG\\GFS.csv, header=FALSE, sep=;) purtest(y, pmax = 4, exo = intercept, test = hadri) I obtained the following message: Error in Math.factor(x) : cumsum not meaningful for factors In addition: There were 50 or more warnings (use warnings() to see the first 50) where's the problem? can you please help me Thank you in advance George __ R-help@r-project.org mailing list 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 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] Problem with plm package
Ah ok. the problem is somewhat great, I can't rewrite it in txt file. On Wed, Jan 29, 2014 at 10:48 PM, george brida george.br...@gmail.comwrote: Ah ok. the problem is somewhat great, I can't rewrite it in txt file. On Wed, Jan 29, 2014 at 9:52 PM, arun smartpink...@yahoo.com wrote: HI George, Some of the attachments get stripped off by the mailing list. Check this link: http://www.r-project.org/mail.html May be .txt files will be okay. Arun On Wednesday, January 29, 2014 3:43 PM, george brida george.br...@gmail.com wrote: Hi Arun you can't open the file? On Wed, Jan 29, 2014 at 9:23 PM, arun smartpink...@yahoo.com wrote: Hi, Your attachment didn't came through. A.K. On Wednesday, January 29, 2014 12:03 PM, george brida george.br...@gmail.com wrote: Dear friends, I have a csv data entitled GFS in the following path: c:\TEG This data is attached with this mail. When I wrote the following lines: y=read.csv(c:\\TEG\\GFS.csv, header=FALSE, sep=;) purtest(y, pmax = 4, exo = intercept, test = hadri) I obtained the following message: Error in Math.factor(x) : cumsum not meaningful for factors In addition: There were 50 or more warnings (use warnings() to see the first 50) where's the problem? can you please help me Thank you in advance George __ R-help@r-project.org mailing list 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 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 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.
