similar to: NAN problem

Hi useR's, Has anyone used the "Lars" package in R before? If so, is there any tutorial (not manual) or worked out example online for this R package that one can go through to figure out how one can use this package with lasso regression? I appreciate any help I can get in this direction. Sincerely, Vishal -- *Vishal Thapar, Ph.D.* *Scientific informatics Analyst Cold Spring

Hi useR's, I am resending this request since I got no response for my last post and I am new to the list so pardon me if I am violating the protocol. I am trying to use the "Kernlab" package for training and prediction using SVM's. I am getting the following error when I am trying to use the predict function: > predictSvm = predict(modelforSVM, testSeq); Error in

Hi All, This is more of a general statistics question that I am not so sure about so I need some thoughts from you experts out there. I have Chip-Seq read count data for 2 different treatments (T1 and T2) and my final aim is to calculate differential expression between them. I also have Chip-Seq read count data from a "control / wild type (C)" experiment. I have 3 replicates (2

Hi, I apologize for posting this here, I am also trying to post this on machine learning emailing lists. I have a set (18K) of sequences (22 nt long) and I have their counts at 4 different stages. The difference in counts from one stage to the next represents how well the sequence performed in the transition. The total counts remain about the same in each stage. So if a 1 sequence loses some

Hi, I'm writing code to generate a plot, in which I draw a series of rectangles. So my code is of the form plot.new() plot.window( ... ) draw rectangle draw rectangle ... Is there a way for me to insert a PNG or PDF graphic at a specific position in the plot (ideally in plot coordinates)? I realize that this might probably be better done in a separate image editor, but if it could

Dear useR's I have a small basic problem which I am hoping to get some help with. I have a data frame, testSeq_df, with 1 row and 500 columns. Each column is a character (a,c,g or t). I want this sequence to have 4 factors (a,c,g,t). When I try the following: for(i in 1:500){ if (length(levels(testSeq_df[,i]))==1) levels(testSeq_df[,i]) <-

Hello, I have a problem that I think can be solved in R but I'm not sure how to tie things together. I have a digital image of a crystal growth in 2 dimensions. And my aim is to calculate the fractal dimension of the crystal. I was planning to use the box counting method. So I need to read in the image in R (for which I can use the pixmap or rimage packages) and then draw a grid over at a

Hi, I am writing to seek some guidance regarding using Lasso regression with the R package LARS. I have introductory statistics background but I am trying to learn more. Right now I am trying to duplicate the results in a paper for shRNA prediction "An accurate and interpretable model for siRNA efficacy prediction, Jean-Philippe Vert et. al, Bioinformatics" for a Bioinformatics project

Hello , My optmatch package is loaded on R and otherwise running fine. I get an error after lcds successfully completes a logistic regression and I then try to obtain a propensity score: pdist <- pscore.dist(lcds) Error: could not find function "pscore.dist" Does anyone know if pscore.dist is part of the optmatch package, or is it contained in another package that I need to

Hello friends, I am trying to test NAS I/O performance over a network and trying to see the numbers for write and read speed. I have successfully configured and ran a lot of tests. However the numbers have not increased, I have reached a bottleneck. I tried playing around with the smb.conf file including all possible variations like, socket options = TCP_NODELAY SO_RCVBUF=(tried values like

I have a dataframe say: date price_g price_s 0.34 0.56 0.36 0.76 . . . . . . and so on. say, 1000 rows. Is it possible to add two columns to this dataframe, by computing say diff(log(price_g) and diff(log(price_s)) ? The elements in the first row of these columns cannot be computed, but

*** BLURB HERE *** Bean Huo (5): fs/buffer: clean up block_commit_write fs/buffer.c: convert block_commit_write to return void ext4: No need to check return value of block_commit_write() fs/ocfs2: No need to check return value of block_commit_write() udf: No need to check return value of block_commit_write() fs/buffer.c | 24 +++++++-----------------

Hello, In a way similar to speexclient voip example, I want to have a voip application where 3 or more users can interact. This creates the following problems "1. How to playback 2 or more audio streams such that they appear overlayed (parallel)". The audio streams are speex encoded. I have few approaches in mind. 1. Create multiple decoder states and multiple ALSA audio playback

Mich besch?ftig folgende Fragestellung. Ich kenne die Verteilung (lognormal) zus?tzlich weiss ich das 99%, das 90% und das 1% Quantil. Gibt es in R eine M?glichkeit die Lognormalverteilung zu finden, das heisst den korrespondierenden logmean und logsd? Vielen Dank f?r ihre Hilfe Gruss Yvonne

Hello R Users, I am doing a Latin Hypercube type simulation. I have found the improvedLHS function and have used it to generate a bunch of properly distributed uniform probabilities. Now I am using functions like qlnorm to transform that into the appropriately lognormal or triangularly distributed parameters for my modes. However I have a parameter which I believe is bimodally distributed,

Hi Matt, All, Thanks for looking into this and your suggestions. I compiled the program with -O3 optimization level (clang -O3) on X86_64 target. After your suggestion, i ran the program with iteration of 10000 runs and found that average runtimes are (attaching data collected as well) Forward loop traverse : 2.006 milli seconds Reverse loop traverse : 1.531 milli seconds Yes, i agree that

date <- as.POSIXlt(Sys.time()) #present date for (i in 1:difftime(as.POSIXlt(Sys.Date()),"2007-10-01")) if (date$wday != 0 & date$wday != 6) {print(date);assign("date", (date-86400))} else (assign("date", (date-86400))) I am trying to print dates from present day to a day in the past, but omitting weekends. I am not doing something right, but can't

Hi Vishal, Ruijie, Thanks for your interest in the project. To get started, the first task would be to merge the 5 patches on top of trunk llvm. The list of patches are listed in the project description: http://llvm.org/OpenProjects.html#llvm_mergesim Please create an account in llvm phabricator (reviews.llvm.org) if you haven't already, and put your patches there. Let me know if you have

I am trying to download a bunch of files from a server, for which I am using download.file( ) within a for loop. The script is working fine except until download.file hits a URL which has no file, at which point it exits. I want to change this behavior to simple log the failure and maintain state within the for loop and iterate to next. I read about try / tryCatch but am having trouble

Can someone throw light on the following problem I am having with RODBC? There's an Excel file I am trying to read from, it has one sheet named 'nameclass'. Thanks in anticipation. Vishal Belsare > library(RODBC) > con = odbcConnectExcel(file.choose()) > tbls <- sqlTables(con) > tbls TABLE_CAT TABLE_SCHEM TABLE_NAME TABLE_TYPE REMARKS 1