search for: dimens

..., # plot margins (b,l,t,r) + las=1) # horizontal labels + plot(matrix, # data to plot + cex=cex_val, # font size + dimen=2 # dimensions to plot + ) + title(main=title) # title of plot + } > # Plot LDA scores with sample names > pic_onscr(lda_result,"Linear Discriminant Analysis") Error in plot(matrix, cex = cex_val, dimen = 2) : error in evaluating the argument 'x...

..., # plot margins (b,l,t,r) + las=1) # horizontal labels + plot(matrix, # data to plot + cex=cex_val, # font size + dimen=2 # dimensions to plot + ) + title(main=title) # title of plot + } > # Plot LDA scores with sample names > pic_onscr(lda_result,"Linear Discriminant Analysis") Error in plot(matrix, cex = cex_val, dimen = 2) : error in evaluating the argument 'x...

Hi there, I am a master student in Denmark, and apply R to analyze the drug data. I use the function 'lda' to obtain a linear discriminant object, and then wish to use 'plot.lda' to do the plot on this object like below. drug.lda <- lda(Inhibition ~ NET_CHARGE + PKA_1 + MW + MLOGP, data = drug.class, method = "moment") plot.lda(drug.lda,

This is a simple problem, but a mystery to me. I'm trying to grab $Family "Scelionidae" from one dataframe and put it into another dataframe occupied with NA in $Family. The result is a "158" ends up there instead of Scelionidae. Simply put fam$Family[1] <- least$Family[1] If I have made a mistake here, can somebody point it out. I've included the simple

...gt; pic_onscr<-function(matrix, title, cex_val=1) > + {x11() + par(mgp=c(5,2,0), # axis margins + # (title, > labels, line) + mar=c(7,4,4,2), # plot margins (b,l,t,r) + > las=1) # horizontal labels + plot(matrix, # data to plot + > cex=cex_val, # font size + dimen=2 # dimensions to plot + > ) + title(main=title) # title of plot + } >> # Plot LDA scores with sample names >> pic_onscr(lda_result,"Linear Discriminant Analysis") > Error in plot(matrix, cex = cex_val, dimen = 2) : error in > evaluating the argument 'x...

Your data frame fam contains factors. Turn it into character strings using fam$Family = as.character(fam$Family) and try again. It may be helpful if you read up on R's factors, see ?factor. HTH, Peter On Tue, Nov 21, 2017 at 2:14 PM, Glen Forister <gforister at gmail.com> wrote: > This is a simple problem, but a mystery to me. > I'm trying to grab $Family

Hi, I have done different discriminant function analysis of multivariat data. With the CV=True option I was not able to perform the predict() call. What do I have to do? Or is there no possibility at all? You also need the predicted values to produce a plot of the analysis, as far as I know. Here my code: pcor.lda2<-lda(pcor~habarea+hcom+isol+flowcov+herbh+inclin+windprot+shrubcov+baregr,

hi, When I try to use the function coordProj {mclust} " coordProj(diabetes[,-1],dimens=c(2,3),what="uncertainy",uncertainty=diabetesModel$uncertainty,parameters=diabetesModel$parameters) " to identify uncertainty, it errors and send this warning message: " Warning message: In coordProj(diabetes[, -1], dimens = c(2, 3), what = "uncertainy", : what...

Well, ?factor does not say anything about this behaviour (assigning numeric code instead of level of factor). And actually if you do assignment for whole vector the result is different (vector in data frame is changed to factor). > temp2$fff[1]<-vec[1] > head(temp2,2) pokus minuty fff 1 T42 240 3 2 T42 300 <NA> > temp2$fff<-vec > head(temp2,2) pokus

Hi !! I am trying to analyze some of my data using linear discriminant analysis. I worked out the following example code in Venables and Ripley It does not seem to be happy with it. ============================ library(MASS) library(stats) data(iris3) ir<-rbind(iris3[,,1],iris3[,,2],iris3[,,3]) ir.species<-factor(c(rep("s",50),rep("c",50),rep("v",50)))

Hi, I am using the lda function from the MASS library. I would to find the decision boundaries of each class and subsequently plot them. I wonder if anybody can offer any help on this topic? Below I applied the lda function on a small dataset of mine. Any help will be much appreciated. > library(MASS) >

...(134.155.91.46) Merging of two data frames by variables with different names (using by.x and by.y) is no longer possible (it was possible in version 1.4.1). Call: mvoteint.classd.matadd <- merge( unique(mvoteint.classd), unique(subset(dimens.mpp, select=c(nation,year,mecode, econlrdiff, eeconlrdiff, welfdiff, labmiddiff,...

I've writen this function (imagenrgb) to display a (m,n,3) array as a RGB image with ngris^3 colors and,optionally, stretching. If option ver=F, it does not display but saves a pseudocolor version of the image as a list (so that subsequent displays are faster). I'd appreciate feedback and improvements and hope that it's useful for others. Example of use: > dim(imatest) [1] 100

Hello, I am using the parameterized EM- Algorithm in clustering context. I use generated datasets, containing 5000 points in 20 dimensions with 5 clusters. The following poses a problem for me: emEst<-me(modelName="EEV",data=data[,1:numDimensions],z=unmap(x)) emEst<-me(modelName="VEV",data=data[,1:numDimensions],z=unmap(x)) coordProj(main="Model EEV",data,dimens=c(1,2),type=&qu...

Hi, Thanks Sarah. Unfortunately I did not get a step further. My question, perhaps a bit clearer, is how to display the case control status (or any other arbitrary point label) after clustering in a plot: With a bit of pseudo code, where dataset is a data.frame, parameters are those column names where we find numerical values (no NAs) and nclasses is the desired number of classes. fit <-

...quality=100, # image quality (percent) # bg="white", # background colour # res=300, # image resolution (dpi) # units="in", width=8.3, height=5.8) # image dimensions (inches) # par(mgp=c(5,2,0), # axis margins # # (title, labels, line) # mar=c(6,6,4,2), # plot margins (b,l,t,r) # las=1 # horizon...

Hi, I'm using plot.lda (R1.4.0 on linux) and, if dimen > 2, I get the error: Error in text.default(x, y, txt, cex = cex, font = font) : X11 font at size 16 could not be loaded The traceback indicates that the problem comes from pairs(). I've tried several cex and font values, always with the same error. The same command works fine with pdf() If I'm missing fonts, how

Dear all When plotting the results of lda and predict.lda, I expect the axis dimensions to remain essentially the same - predict.lda projects new observations onto linear discriminants. However, I am getting different ranges of values along the x-axis when producing plots following each: plot(lda) yields values in the range approximately -3 to +4, while plot(predict.lda) yields va...

Hi all, I am trying to plot the results of a discriminant analysis done with lda(MASS) but my groups appear in two different plots (in the same graphics device) and I want to combine them in one plot. My code looks like: BirdTrain.lda <- lda(Bdisperser~., data=BirdTrain.mx) predict(BirdTrain.lda) plot(BirdTrain.lda) I have two types of Bdisperser, so I only get one linear discriminant

Dear R I'm trying to plot the lda means onto a 2 D plot of discriminant scores. Preferably I'd like these to be in a larger font compared to the discriminant scores. I tried skull.mean.pred <- predict(skulls.lda, as.data.frame(skulls.lda$means), dimen=2) from which I got skull.mean.pred $class [1] 1 2 3 4 5 Levels: 1 2 3 4 5 $posterior 1 2 3 4