similar to: Prediction with two fixed-effects - large number of IDs

I have no direct experience with such horrific models, but your formula is a mess and Google suggests the biglm package with ffdf. Specifically, you should convert your discrete variables to factors before you build the model, particularly since you want to use predict after the fact, for which you will need a new data set with the exact same levels in the factors. Also, your use of I() is

Hello Jim, Thanks again. I am getting the two-headed arrow but I cannot seem to get the coordinates right for the arrow to appear beneath the map. These coordinates puts the arrow on the left hand side. Thanks again! Sincerely, Milu On Tue, Apr 12, 2016 at 1:15 PM, Jim Lemon <drjimlemon at gmail.com> wrote: > Hi Milu, > There is a two-headed arrow on the image you sent, and it

I have a big data set that includes character variables of many different values. I'm trying to use ff to read the data and then use biglm.big.matrix to build linear models. However, since big.matrix will convert all character vectors to factors and the character labels will be lost. I decided to create a lookup table outside of R for my character columns and use numbers to represent different

Dear David, The device was the issue. The quartz() device works fine but pdf() does not. Now I just need to figure out the limits for map for Europe. Thanks for all your help and patience. Sincerely, Milu On Sun, Apr 10, 2016 at 7:10 PM, David Winsemius <dwinsemius at comcast.net> wrote: > > > On Apr 10, 2016, at 4:12 AM, Miluji Sb <milujisb at gmail.com> wrote: > >

Hi Milu, My fault here. As I don't have the data to make the map and try out my suggestions I mixed up the x and y coordinates. Try this: par(xpd=TRUE) arrows(-19.75966,53,33.60000,53,code=3) par(xpd=FALSE) Jim On Tue, Apr 12, 2016 at 10:11 PM, Miluji Sb <milujisb at gmail.com> wrote: > Hello Jim, > > Thanks again. I am getting the two-headed arrow but I cannot seem to get

Hello David, This is exactly what I want but I still can't get the arrows. R and R studio is updated. Thanks again! Sincerely, Milu On Sat, Apr 9, 2016 at 10:29 PM, David Winsemius <dwinsemius at comcast.net> wrote: > > > On Apr 9, 2016, at 1:27 PM, David Winsemius <dwinsemius at comcast.net> > wrote: > > > > > >> On Apr 9, 2016, at 11:18 AM,

This is what I was looking for. Thank you everyone! Sincerely, Milu <https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail> Mail priva di virus. www.avast.com <https://www.avast.com/sig-email?utm_medium=email&utm_source=link&utm_campaign=sig-email&utm_content=webmail>

> On Apr 10, 2016, at 1:45 PM, Miluji Sb <milujisb at gmail.com> wrote: > > Dear David, > > The device was the issue. The quartz() device works fine but pdf() does not. Now I just need to figure out the limits for map for Europe. Thanks for all your help and patience. After plotting a map of Europe with base graphics the coordinates of the lower-left and upper-right

Dear Jim, Thanks again! I do want the arrows at the bottom (beneath the map). This is what I am doing: # Draw the map eps_europe <- mapCountryData(n, nameColumnToPlot="eps_score", mapTitle="EPS Score - Europe",colourPalette=colourPalette, catMethod="fixedWidth", missingCountryCol = "white", mapRegion="Europe", addLegend=FALSE) # ISO3 codes on

Ista, et. al: efficiency? (Note: I needed to correct my previous post: do.call() is required for pmax() over the data frame) > x <- data.frame(matrix(runif(12e6), ncol=12)) > system.time(r1 <- do.call(pmax,x)) user system elapsed 0.049 0.000 0.049 > identical(r1,r2) [1] FALSE > system.time(r2 <- apply(x,1,max)) user system elapsed 2.162 0.045 2.207 ##

> On Apr 10, 2016, at 4:12 AM, Miluji Sb <milujisb at gmail.com> wrote: > > Hello David, > > This is exactly what I want but I still can't get the arrows. R and R studio is updated. Thanks again! I didn't try it in Rstudio until just now (and I don't remember that you ever mentioned RStudio as a possible issue.) The plotting I see in the default graphics Rstudio

or in base R : ?xtabs ?? as in: xtabs(~previous_location + current_location,data=x) (You can convert the 0s to NA's if you like) Cheers, Bert 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 Tue, May 8, 2018 at 9:21 AM, Huzefa

> On Apr 9, 2016, at 11:18 AM, David Winsemius <dwinsemius at comcast.net> wrote: > > >> On Apr 9, 2016, at 10:46 AM, Miluji Sb <milujisb at gmail.com> wrote: >> >> Dear David, >> >> Thank you for your answer. Sorry for the embarrassing mistake. However, even with when I generate a map for the whole world using: >> >> eps <-

xtabs does this automatically if your cross classifying variables are factors with levels all the cities (sorted, if you like): > x <- sample(letters[1:5],8, rep=TRUE) > y <- sample(letters[1:5],8,rep=TRUE) > xtabs(~ x + y) y x c d e a 1 0 0 b 0 0 1 c 1 0 0 d 1 1 1 e 1 1 0 > lvls <- sort(union(x,y)) > x <- factor(x, levels = lvls) > y <-

On Tue, Feb 20, 2018 at 11:58 AM, Bert Gunter <bgunter.4567 at gmail.com> wrote: > Ista, et. al: efficiency? > (Note: I needed to correct my previous post: do.call() is required for > pmax() over the data frame) > > > x <- data.frame(matrix(runif(12e6), ncol=12)) > > > system.time(r1 <- do.call(pmax,x)) > user system elapsed > 0.049 0.000

Dear all. Is it possible to convert.identify iso3 country names to the seven continent names? # Asia, Africa, Antarctica, Australia, Europe, South America, and North America, I have tried the following: ### region <- merge(countryExData,df,by.x='ISO3V10',by.y='iso3') where df is the name of my dataset with iso3 the identification variable but there seems to be a a lot of

Dear David, Thank you for your answer. Sorry for the embarrassing mistake. However, even with when I generate a map for the whole world using: eps <- mapCountryData(n, nameColumnToPlot="eps_score", mapTitle="EPS Score",colourPalette=colourPalette, catMethod="fixedWidth", missingCountryCol = "white", addLegend=FALSE) And then use:

Dear all, I am trying to use 'bam' to run the following generalized additive model: ### m <- bam(result ~ factor(city) + factor(year) + lnpopulation + s(lnincome_pc) + ,data=full_df,na.action=na.omit,family=ziP(theta = NULL, link = "identity",b=0)) But getting the following error: ### Error in bam(result : extended families not supported by bam The documentation for

Hi Milu, There is a two-headed arrow on the image you sent, and it seems to be where you specified. Did you want it beneath the map, as: par(xpd=TRUE) arrows(-22,54.75,-22,74,code=3) par(xpd=FALSE) Jim On Tue, Apr 12, 2016 at 7:58 PM, Miluji Sb <milujisb at gmail.com> wrote: > Dear Jim, > > Thanks again! I do want the arrows at the bottom (beneath the map). This is > what I am

Hi Milu, I just realized that by "the bottom of the map" you may mean "beneath the map", in which case you should use: par(xpd=TRUE) arrows(...) par(xpd=FALSE) Jim On Mon, Apr 11, 2016 at 11:50 PM, Miluji Sb <milujisb at gmail.com> wrote: > Dear David, > > Thank you very much for your replies! I didn't know about par('usr'). > > I get