IdentifiantMot de passe
Loading...
Mot de passe oublié ?Je m'inscris ! (gratuit)
Navigation

Inscrivez-vous gratuitement
pour pouvoir participer, suivre les réponses en temps réel, voter pour les messages, poser vos propres questions et recevoir la newsletter

R Discussion :

Question Rstudio : PLS regression


Sujet :

R

Mode arborescent

Message précédent Message précédent   Message suivant Message suivant
  1. #1
    Membre habitué
    Homme Profil pro
    Bioingénieur
    Inscrit en
    Décembre 2015
    Messages
    10
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Âge : 32
    Localisation : Belgique

    Informations professionnelles :
    Activité : Bioingénieur

    Informations forums :
    Inscription : Décembre 2015
    Messages : 10
    Par défaut Question Rstudio : PLS regression
    Bonjour à tous,

    Je suis débutant en Rstudio et j'aimerai faire une Partial Least Square regression. Cependant, j'ai un problème avec cette fonction. J'ai fait l'exemple proposé par R avec les données de "gasoline" et tout marchait. Mais, une fois que je veux faire les mêmes manipulations avec mes données, cela ne fonctionne plus.

    J'ai donc plusieurs questions concernant cette fonction.

    QUESTION 1

    Mes données que vous trouverez en pièce jointes sont composés de 4 variables avec 2200 observations. J'aimerai construire un modèle permettant de prédire la "nitrogenConcentration" en fonction de la "reflectance" des différentes "lettreFiltres" ET de savoir quels "filtres" sont les plus pertinents dans ce modèle.

    Mon code R est le suivant :

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12
    13
    14
    15
    16
    17
    18
    19
    20
    21
    22
    23
    24
    25
    26
    27
    28
     
    # 1. importation donnée
    # -------------------------
    data1 <- read.table("tableReflec.txt", header = TRUE)
     
    # 2. unstack data
    # -----------------
    #unstack nitrogen content
    unstack_Nc <- unstack(data1,nitrogenContent ~ lettreFiltre)
    Nc <- unstack_Nc[,1]
    remove(unstack_Nc)
    #unstack reflectance
    unstack_reflec <- unstack(data1, reflectance ~ lettreFiltre)
    myName <- paste("reflectance",names(unstack_reflec),sep = ".")
    names(unstack_reflec)<-myName
    remove(myName)
     
    # 3. Données comprenant Nc et la reflectance pour les différents filtres
    # ------------------------------------------------------------------------------
    #ré-organisation des données
    dataReflec <- data.frame(Nc)
    dataReflec <-cbind(dataReflec, unstack_reflec)
     
    # 4. PLS
    # -------
    # tools -> install package -> PLS
    library(pls)
    model1 <-plsr(Nc ~ reflectance, ncomp = 10, data = dataReflec, validation = "LOO")
    Cependant, quand je lance la dernière ligne, j'ai une erreur qui est la suivante :

    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
     
     Error in eval(expr, envir, enclos) : object 'reflectance' not found
    Je ne comprends pas pourquoi cette fonction fonctionne avec "gasoline" et pas avec mes data. Quelqu'un peut-il m'aider ? Qu'est ce que je dois faire pour que cela fonctionne ?


    QUESTION 2

    En faisant des recherches sur internet pour mon problème, je suis tombé sur une autre fonction qui est la plsreg1 qu'il est possible d'obtenir avec le package "plsdepot". Quelqu'un peut m'expliquer la différence entre cette fonction et la fonction plsr ? De plus, j'ai l'impression qu'il y a peu d'information sur la plsreg1. Je n'ai pas trouvé d'exemple utilisant cette fonction sur le net. Si quelqu'un a des infos, je suis preneur !


    QUESTION 3

    Je cherche à connaitre les filtres les plus pertinents à partir de la PLS. Quelle serait la fonction me permettant de savoir cela ?


    Voila, je sais que cela fait beaucoup de question et je vous remercie d'avance pour le temps que vous consacrerez pour y répondre.
    Fichiers attachés Fichiers attachés

+ Répondre à la discussion
Cette discussion est résolue.

Discussions similaires

  1. aide library(pls) question sur data(yarn) plz
    Par redando dans le forum R
    Réponses: 3
    Dernier message: 20/09/2016, 11h28
  2. Réponses: 3
    Dernier message: 31/08/2016, 10h58
  3. Réponses: 6
    Dernier message: 17/06/2013, 15h23
  4. Question de faisabilité
    Par lisarasu dans le forum CORBA
    Réponses: 3
    Dernier message: 14/05/2002, 11h26
  5. [HyperFile] 2 questions de débutant
    Par khan dans le forum HyperFileSQL
    Réponses: 2
    Dernier message: 29/04/2002, 23h18

Partager

Partager
  • Envoyer la discussion sur Viadeo
  • Envoyer la discussion sur Twitter
  • Envoyer la discussion sur Google
  • Envoyer la discussion sur Facebook
  • Envoyer la discussion sur Digg
  • Envoyer la discussion sur Delicious
  • Envoyer la discussion sur MySpace
  • Envoyer la discussion sur Yahoo