similar to: Decision boundaries for lda function?

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,

Hello I have serious problem with zhone z-plex 10 the device sometimes get red alarm on some channel and then the device suddenly shutdown. I try to unplugged the power chord and plug again, but the device still didn't up. I have to wait for long time to get this device up again. How can i fix this problem, so i can get this device immediately? Eris Riswanto

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 R experts. I performed an Linear Discriminant Analysis (lda) and now I want to plot the first two axes (LDA1 and LDA2). Well MASS package have the plot.lda and pairs.lda to do that. But, they don’t let me personalize them, once they don’t accept the type=’n’ plot. So I start looking at the lda object properties, trying to found out my groups (n=4) axis coordinates, curiously I couldn’t. There

Hello all, I''ve been running Xen 2.0x for some few months and have written a simple top(1)-like tool - xenConsole - which continously displays all running domains, along with their CPU-consumption, load averages, percentual CPU-/memory-usage, domain uptimes and so forth. A screenshot probably explains it better than me: http://www.gothix.dk/xenConsole/xcon-0.10b1-ok.PNG It requires

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

# lda.r > # > # Author: Amsha Nahid, Jairus Bowne, Gerard Murray > # Purpose: Perform Linear Discriminant Analysis (LDA) > # > # Input: Data matrix as specified in Data-matrix-format.pdf > # Output: LDA plot > # > # Notes: Missing values (if any) are replaced by the half of the lowest > # value in the entire data matrix. >

hello, i already include the error in blue color word. i hope it can help you to understand my question. if not burden you, please give me a guide how to correct the error or maybe you can correct the coding cause error. thank you. > # lda.r > # > # Author: Amsha Nahid, Jairus Bowne, Gerard Murray > # Purpose: Perform Linear Discriminant Analysis (LDA) > # > #

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

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

HELLO, anybody... could you help me to check the below coding for volcano. what is the mistake? what the plot could not display? # volcano_plot.r # # Author: Amsha Nahid, Jairus Bowne, Gerard Murray # Purpose: Produces a volcano plot # # Input: Data matrix as specified in Data-matrix-format.pdf # Output: Plots log2(fold change) vs log10(t-test P-value) # #

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

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 <-

Dear R users, I want to do discriminant analysis on my data. I have successfully followed the discriminant analysis in V & R on the iris data: > ir <- rbind (iris3[,,1],iris3[,,2],iris3[,,3]) > ir.species <- c(rep("s",50),rep("c",50),rep("v",50)) > a <- lda(log(ir),ir.species) > a$svd^2/sum(a$svd^2) [1] 0.996498601 0.003501399 > a.x <-

I found the following in eqscplot, library MASS: When given just one point, it fails: > eqscplot(x=1, y=1) Error in plot.window(xlim, ylim, log, asp, ...) : need finite xlim values After inserting browser() just before the last line in eqscplot() (in which plot() gets called), I found the following: > eqscplot(x=1, y=1) Called from: eqscplot(x = 1, y = 1) Browse[1]> xlim [1] NaN

This is not a bug. There is no way one can scale the plot axes equally when no indication has been given of the scaling expected for either! It has always seemd to be that plot(x=1, y=1) shoul dnot work: it has not supplied enough information! On Thu, 8 Nov 2001 pflugshaupt@geobot.umnw.ethz.ch wrote: > I found the following in eqscplot, library MASS: > > > When given just one

Dear all I was writing some code that needed functions from packages 'MASS' and 'CircStats', and I received a warning saying that the function 'eqscplot' in one of the packages was masked by another 'eqscplot' from the other package (in fact 'eqscplot' from 'CircStats' seems a short version of 'eqscplot' from 'MASS'). This masking of

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

> # >> # Generate the figures (on screen) >> # >> # Image generation - function definition >> 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 + >