similar to: Trouble with data dim

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

Dear r-help, I am trying to build a new function (to process rpart objects) that will output matrix that has a row for each node and a column for each feature. With each entry in the table is a numerical property at that node for that feature (e.g. surrogate split agreement, improvement). My current trouble is that the only clue to the identity of the feature is stored as the *name* of the

Dear All, I was trying to follow some of the examples in Venables and Ripley "Modern applied ... with S-plus" I have downloaded a copy of the iris data set and loaded into R. : however I cannot use the apply command (from p47): > apply (iris, 2 ,mean) Error in sum(..., na.rm = na.rm) : invalid "mode" of argument > apply (iris, c(2) ,mean) Error in sum(..., na.rm =

Perhaps this has been mentioned and I wasn't paying attention. When I try to import a postscript file (which is EPS according to the R help on postscript) my Framemaker filter chokes. If I simply edit the first line of the file to put %!PS-Adobe-3.0 EPSF-3.0 rather than %!PS-Adobe-3.0 it works just fine. It sure would be simpler if the postscipt driver put this in the file. Paul Gilbert

Dear r-help, Is there any available description of the components of lqs objects found in the package "lqs"? > names(slts) [1] "crit" "sing" "coefficients" "bestone" [5] "fitted.values" "residuals" "scale" "terms" [9] "call"

Courtesy follow up back to the unix-dev list (mainly for those like me who searched the archives) Thanks to the very swift response of Darren, a patched auth.c did the trick Patch enclosed below (apply after the p20) I'm guessing he'll either update p20 or issue p21 soon. Many thanks Andrew --- auth.c.orig 2003-06-25 23:14:16.000000000 +1000 +++ auth.c 2003-06-25

Dear All, I am trying to use ode solver "rk4" to solve an ODE system, however, it keeps saying: Error in eval(expr, envir, enclos) : object "dIN" not found. The sample codes are enclosed as follows, please help me. Thank you very much! rm(list=ls()) library(odesolve) # The ODE system ode <- function(t,x,p){ with(as.list(c(x,p)),{

Dear All, I am trying to use ode solver "rk4" to solve an ODE system, however, it keeps saying: Error in eval(expr, envir, enclos) : object "dIN" not found. The sample codes are enclosed as follows, please help me. Thank you very much! rm(list=ls()) library(odesolve) # The ODE system ode <- function(t,x,p){ with(as.list(c(x,p)),{

All - I have a large data frame that looks like ID p1 p2 p3...p20 Lat1 Lat2 Lat3...Lat20 Long1 Long2 Long3...Long20 1 0 0 1 0 NA NA 29.xx NA NA NA -89.xx NA 2 1 0 0 1 27.xx NA NA 29.00 -88.00 NA NA -89.xx 3 0 0 0 0 NA

Buenas, Se me hace raro que la función no admita un parámetro para eso. Una opción es hacerlo en dos pasos: P1: Reemplaza p10, p20, ... por algo como q10, q20, ... P2: Haz el reemplazo que quieres. Igualmente puedes hacer el reemplazo que haces y en un segundo paso reemplazar p10_integra_datos por p10_cuales_rep, ... Un saludo, Manuel. El jue, 09-03-2023 a las 10:20 -0300, juan manuel dias

I have written a few different summary functions. I want to calculate the statistics by groups and I am having trouble getting the output as a dataframe. I have attached one example with a small dataset that calculates summary stats and percentiles, I have others that calculate upper confidence limits etc. I would like the output to be converted to a dataframe with one of the columns as the

Hi, I'm having a problem using sample() within a function. Basically I get an error reading: Error in sample(v, 1, prob = h) : non-positive probability Can anyone advise me as to the possible origin of this error? Here is my code #Discretised Gillespie algorithm function (From SMfSB, D.J. Wilkinson) gillespied=function (N, T=100, dt=1, ...) { tt=0 n=T%/%dt x=N$M

Hola, Algo que me está pasando con esta función es lo siguiente: A "p1" lo tengo que reemplazar por "p1_integra_datos" y a "p10" por "p10_cuales_rep", pero como p10 contiene p1, y así pasa con otras preguntas (p2 y p20, por ejemplo), me los reemplaza en ambas, cuándo lo que busco es otra cosa. reglas$condicion_final <-

I've been struggling to get a loop to work. I want to create a new variable in each loop with data from some function. for example: # part of the names for the variables to be created Frags <- c("F04", "F05", "F07", "F09", "F11", "F13", "F14", "F17", "F18", "F19", "P20",

We have FreeBSD-server with dovecot installed on it as IMAP-server. My user and password database is a text file with plaintext passwords. Clients connect to imap-server via TLS protocol and plaintext password. All works fine. But I want to configure ability to authorize with a client certificates. I have generated a client certificate and imported it to email-client. Also I have configured

John R Pierce wrote: > On 5/25/2016 11:44 AM, Kelly Lesperance wrote: >> The HBA is an HP H220. > > OH. its a very good idea to verify the driver is at the same revision > level as the firmware. not 100% sure how you do this under CentOS, my > H220 system is running FreeBSD, and is at revision P20, both firmware > and driver. HP's firmware, at least what I

On 11/05/2016 06:40 AM, Steven Tardy wrote: > What firmware is/are on the cards? In the dmesg output I included: mpt3sas0: LSISAS3008: FWVersion(10.00.00.00), ChipRevision(0x02), BiosVersion(08.25.00.00) > The driver version difference between C6 and C7 could interact badly with old firmware. I suspected that was the case, and tried running elrepo's "kmod-mpt3sas" from

Dear List, This is just a programming problem which i cannot seem to figure out. I am trying to get a set of power from a test (say, kolmogorov smirnov) out of a distribution (say, G-K distribution) as follows. I am trying to reduce to pain of writing the whole set of data points (p# below) using "for" loop. However, I seem to have some problem in it as the output "M" does not

I have added "auth_debug_password=yes" to "10-logging.conf" and restarted dovecot.?But I do not see any information about the password in the logs. Does this mean that the thunderbird does not send the password??Although it asks for the password and I enter one. New log: dovecot: master: Warning: Killed with signal 15 (by pid=19769 uid=0 code=kill) dovecot: master: Dovecot

Hi useRs - I was wondering if anyone out there can tell me where to find R-code to do mixes of tables and graphics. I am thinking of something similar to this: http://yost.com/information-design/powerpoint-corrupts/ or like the excel routines people are demonstrating: http://infosthetics.com/archives/2006/08/excel_in_cell_graphing.html My aim is to provide small graphics to illustrate