search for: can1

...pose Map1 (which has black borders between Canadian provinces) onto Map2 (which is also a map of Canada)? Thanks a lot for your hints! Dimitri ### A. Reading in Canada data at the province and then at the county level: library(raster) getData('ISO3') # Canada's code is "CAN" can1<-getData('GADM', country="CAN", level=1) can2<-getData('GADM', country="CAN", level=2) ### B. Creating Map1- with borders between provinces: can1 at data[["groups"]]<-1 # want all provinces to have the same color, but ideally I'd like th...

...' ' 1 this is different than the output I get with SAS: Eigenvalue Difference Proportion Cumulative Ratio F Value Num DF Den DF Pr > F 1 0.1068 1.0000 1.0000 0.90349416 31.88 2 597 <.0001 I am also wondering how to plot the can1*can1 like it is done in SAS. proc plot; plot can1*can1=species; format species spechar.; title2 'Plot of Constits_vs_cassettes'; run; Thanks [[alternative HTML version deleted]]

...'m stumped. > # fit randomized block model, .~Block+Gp, where Gp = Contour:Depth > # Soils data is in car package > library(car) > soils.mod1 <- lm(cbind(pH,N,Dens,P,Ca,Mg,K,Na,Conduc) ~ Block + Gp, + data=Soils) > # Do the canonical discriminant analysis for Gp: > can1 <-candisc(soils.mod1, term="Gp") > > str(can1$scores) 'data.frame': 48 obs. of 11 variables: $ Block: Factor w/ 4 levels "1","2","3","4": 1 2 3 4 1 2 3 4 1 2 ... $ Gp : Factor w/ 12 levels "D0","D1",&q...

...6 5.47 2 2.58 4.46 2 2.16 6.22 2 3.27 3.52 I put these in a table and then a linear model >newdata <- lm(cbind(three, five) ~ species, data=rawdata) and then do a candisc on them >candata<-candisc(newdata) Here are my scores; >candata$scores species Can1 1 1 -2.3769280 2 1 -2.7049437 3 1 -3.4748309 4 1 -0.9599825 5 2 4.2293774 6 2 2.6052193 7 2 2.6820884 and here are my coefficients > candata$coeffs.raw Can1 three -5.185380 five -2.160237 > candata$coeffs.std Can1 three -2.5...

I am trying to create a frequency distribution and I am a bit confused. Here are the commands I have entered: > data <- read.csv(file="40609_sortedfinal.csv",head=TRUE,sep=",") > NumberOfActionsByStatus = data$STATUS > NumberOfActionsByUser = data$ETS_LOGIN > NumberOfBidOffer = data$BID_OFFER > NumberOfActionsByUser.freq = table(NumberOfActionsByUser) >